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\ Computer is like your 
M ony 8 


BHukBau 
IIQY DRAWERS. 


drawers for numbers 
with decima\ points. 
AC hese are called 


REAL or floating pot 
numbers. TAAL 
and it has teeny ting ators 


en late er numbers, 


decimal points or fractional parts, 
ntegers ave used for Gountin g 
and things like that. 


“Lo talk about the contents of a 


big drawer, give the \ 


drawer i fe) 
N bolia =E 8 - 
gable YS. 


Cant it H0G A }i 
because that isn’t nice and because it has more 
than S[X characters. You can make vp Names 
Using any Combination of \etters or numbers but 
Ye Olde Gomputer restricts you to names with Six 


Gharacters or less, just so you'll know it’s boss. 
“rhe first character inavariable name must 
be a letter, however. Variable names can be made 
Up to suit your faney (if you wear a 


fancy), 
One thing more: REAT variable names can't 
Start with 
ral) 


{,J. if 36, Me oN 


TW ake up nice Mnemonic names Starting with 


Lv, WL, 5 or N for the cute 
little integer drawers. Vol Oy a 


N\Q are o.k. as integer 


Variable names if your mnemonica has a 


tin ear, bout MAXVAL rT ASX 


might have more ASSOGIATIVE valve to 
somebody trying to understand your program, 
Who knows ° The Somebody someday could be you! 
ee m iS a good , mnemonic two character 
YAN variable name, ts Conese and 
descriptive -not flowery, 
(Awright. So it is flowery.) 


IN] one since WI doesn’t start with one of 


the Special Letters, it wovld be a good Rea/ 
Variable name for a location in which te Keep 


3. 1415 9265 BE BI 73 © Memnamnin a mei mene 


a i a, 


Whad AYO mean? 


Actually, momma's drawers aren't big 
enovagh to store an | NN Ss = number 
of significant digits. LENT & 
iT here’s only room for about seven Significant 
digits, 50 you'd need to be happy with 
3.141599 Be that as itmay, IP Jis still 
a good name for it and it 1S plenty good 
enovgh for Government Work where errors 
of plus or minus 5 Billion ave acceptable. 


JH FPA is a egal Inecse Vasunne 


NAME, since it has 


“= { toG characters 
<> STARTS with an I, and 


=P the other characters are 


alphabetic or numeric 
lbut not recherché exotica 
suchas dollar signs $ $$ or 


ampersands Be & & .¢ So that’s what 
those things are: 


Say, that’s the 


strangest computer 
Tve ever seen! What 


brevity, you can use a as 
a variable name, if you wish, 

Choosing a smale character for brevit 
brings to mind Hamlet. LH owever, the FORTRAN 
compiler would consider HAMLET @ six 
character real name. or are typical 
one character real variable names and or 
[Ml are possible one character integer names. 


But I digress. 
A\\ your digressions a YO} 


are giving me terrible 
in digression. 


Vn certain situations IHT! P2 might have 
mnemonic valve, but it would have dubious 
merit ina Kinetoelastostaticmagnetohydro- 
quvasiviscopolumeric computer program. 


He dov’t say ! 
a nese are called VARIA PLE names 


because tne Name sticks with the drawer but 
the Gontents of the drawer May change. 


ZY vow vo FILL 
LS OUR PRAWERS 
Cy > Wh 


To put a number ina drawer have the Name 
of the drawer appear on the ere of an 
equal sign, hus, 


PI=3.1415 9 


stores the yeo\| or floating point constant 
value 3.14159 in the variable location named FI. 


Similarly, 


BIGX = 2). 


stores the real constant 25, (it has a 
decimal point, so it’s real, didya notice 7?) 
in the location named B\GX. 


This is called an 


2° Assignment Statement 


arhe quantity on the Right of the equal 
sign Re places the contents of the draver, 


_ a 
The old contents areflost forevermore, 
When a variable appears on the [rant of an 


WN ttipt till jeg 
equal Sign man assignment statement, its 
Valve is copied from the storage location 
with that name but the contents of the location 


remain unchan ged. 


=A Ff A . 
a 

a 

Ze Give thema 


| 


Simple Example 


(Los deer, Se 
them a Simple 


Example / 


A\\ right, Momma. 
(\\ They’ re your drawers | 


N = (N+3)*MOM 
is an Integer A ssignment Statement, Don't 
think of the equal Sign as being an equal sign 
but think of it as being a gozzinta sign. Ot 
MEANS the computer should first Figure out 
the value of what’s on the right. 


the stuff on the right has been boiled down 
toa single number. lHlere’s where the 
gozzinta Sign Swings into action. Once the 
computer has oa single number on the right, 
— that valve gozzinta the 


location onthe left. Dn this case, the location 
on the left is the drawer named “N” The 
old value of N 1s thrown away and the new 
value is stored, Thus, the value of N changes 
during the gozzinta part of this operation, 
The valve of MOM 1S used but remains 
unchanged. 

This definition of gozzinta differs 
from the rigorous mathematica) definition 
with which you are undoubted|y familiar. 
that is, four gozzinta seven not quite twice, 

ZNs another example, Suppose a Variable 
named THAT contained the number 10.5 
and a variable named THI contained-8.9, 
If the computer then came to the assignment 
statement THAT = THIS it is clear 
that THIS ond THAT would both then 
contain - 6.0, Thus , this THAT takes on the 
valve Formerly in that HIS, So THIS and 
THAT end vp with the same contents but 
those contents of that THIS are undisturbed, 1.e. 


itis these contents of that THAT that get 
replaced, since that THAT is what ig onthe 
left of the equal sign. This is, of course 
clear to anyone whe Knows one hand Cate 


other 


Do you Know one 
hand from the other ? 


se OA 
WEA > An equa sign in 
ForRTRAN 
isn’t at all like 
an equal sign in 


Al\gebral 


(ln FORTRAN it means 
pire: Figure ovt the valve of what. 
La 


on the Right Side &+ 


Ther Store that valve as the new 
contents of the location 
named on the Left. 


\ 


its you try to assign Oo floating point 
number to an integer location, there isn’t 
enough space, so the fractional part gets 
scraped off and lost on the floor, 


4 


hus, this Sequence of statements first 
gives ml the value 3, ae then gives aT 
the valve 3 > and Finally gives Pi the 
value 3, 0000cee- 

a cee ae 


Don’t M1 % integers and reals on the right 
’ 

hand side of an equal sign or you'll be 

making a big mistake.! 

TA pologies to Alka Seltzer” 


* This joke goes back before your time. Ask your folks to 
explain it, if they remember the ear |y days of television, 


BESIEES BARRY 
GOLD WATER, WYAT IS 
FOUND TO DHE RIGHT? 


Whe stuff on the right of an equal sign is 
an Arithmetic Ex pression which boils 
down toa Single musrmber that 
can be stored in the location indicated on 
the left. A\ typical assignment statement 
might be 

SUM =PEHTOTAL*3.95 - OLDX 


Arithmetia Expression 

ZAM variables in an arithmetic 
expression (must) have been 
assigned valves earlier in the program, 
hus, the expression can be evaluated 
to yield o Single Number which is 
then stored as the contents of the 
location named “ SUM: 


Jn FORTRAN, valid arithmetic operators 
which can be used in ways obvious to even the least 
literate perci pient of this cogent, efficacious 
tome for the preschooler are: 


+ SS Which means ADD 


—_ FEE Which means SUBTRACT 
ZB Which means MULTIPLY’ 
a TEE Which means DIVIDE 


ZX short Parenthetical Note is in order: (((Paren- 
theses can be used to group terms and to make 
your intentions (which we hope are honorable) 
clearer. )(Be sure every left parenthesis has a 
Corresponding right parenthesis or youll be 
Subjected to unnecessary Humiliation and 
Ridicule, ) The FORTRAN compiler So Stickler 
for suntactical precision, )) 
(Ess of noted 
CF 
A\ less obvious arithmetic operator is two asterisks 
in succession, which means “raise to the power, “Thus, 
A** 50 


means “raise X to the 50% power,” 


*% 1S Sometimes called the Nathan Hale operator 
since he had but one asterisk. 


Just three pages ago, Dr. KX (as he is known in 
Washington) told you not to mix integers and 
reals inan arithmetic expression. What’s called 
a Mixed Mode” expression. On gome computers 
it’s legal, but it can often cause more trouble 
than a “Code in da Node.” So, 1f you shouldn't 
mix integers and reals in an expression, what 
gives with Kx*X5O meaning e e 


© ATART of PAINFUL EXPONENTIATION | 
3 EYPLANATION 


: (14 you hate it- skip it!) 
(If your teacher likes it- flunk! ) 
CS 


Exponentiation 1S@ Special case.When 
possible, try to use Integer Exponents, 
ALPHA **1 

isn’t considered a mixed-mode expression. 

Z It simply means multiply ALPHA by itself 
I times, That’s a cinch for the Computer, 

ou can also raise a. real toaveal power, 

as in ALPHA ** BETA 

Z This isabit-more hassle to figuve out Whe 
poor computer needs tofind the log of 

: ALPHA, multiply it by BETA, and then 


un-logify the whole mess as in ePeAAnsitin 


My learned colleague, 
the famous mathematician 
C heby shev fvom St Yetersburg 
(his name is sometimes spelled 
Tschebytschew or Tehebycheff but 
his friends call him “Pafnuti;’) 
feels the term “unlogify "is 
only acceptable to pure mathe- 
maticians. Appli ed scientists 
shovid use the term "to dislog.” 


Piven basssass 


Last but not least, things equally high on the 
pecking order get evaluated from Left to R; hy 
When in doubt, throw in some parentheses | 
and be sure, @nlu use good quality parentheses 
with nice round sides. 


Arhus, amazinaly enovgh, most things on the computer 
work the way you'd expect them to. 


L$ you wrote the FORTRAN expression 
A+B/O** 1.5 *DELTAtE/(F+G) 


it would be evaluated, according to the hierarchy, 


aS A+ BAvELTA 2. E 
CG F+G 


On the other hand, for had you written 


(A+(B/C)**(1.5 XDELTA)#E)/(F+@) 
it would be figured as 


_ (By PENTA 
F+G 


Bn a demented sort Away, this all makes sense, 
doesn’t it ? One thing you might not expect, how- 
ever, is that INTEGER DIVISION chops 
off and throws away anu Fractional remainder 
Thus 10/3 comes ovt 3 in FORTRAN, while 
10./3. Comes out 0.$ 3.333333,,. Similarly 
-372 comes ovt as -1 while -3./72. would 
come ovt -1.5 ln the olde daus, 36/24/36 
would come out Marilun Monroe, but now it 
comes out Ghavvinist. In algebra, Something / 0 
comes out infinity but onthe computer it comes 


ovt ERROR MESSAGE 


Se <p" 
(OC jo 
) 
CT | \ =a 


LA library of common|y used functions is stored 
inthe computerso yov don't need to develop 
your own programs to find Sines, absolute 
valves, Square Yoots, and similar everyday 
trivia, For instance, to find the Square Root 
of a Variable XxX, You simply write 


SQRT (X) 
The sine of a + ot Might be written 


SIN (P1/4.+ ALPHA) 


A botanist, writing (oN FORTRAN program to Study 
daffodils might write 
SQRT (DAFDIL) 


OF course, Some daffodils don’t have Sqvare roots. 
If you plant a negative DAFDIL in the SQRT 


function, you re apt to find your whole program 


potted! Ve 


SQRT (DAFDIL) KO 2 
<< 
PRIIORA, oe 
he quantity in the parentheses of these 
library functions is an arithmetic EXPvession. 
(In this case it’s a very simple expression, consisting 
of just the one variable name.) The expression is 
first evaluated. Whe VALUE of that 
eae ed 
expression is the ARGUMENT of the function 
and is PACE y, off to a secret location in 
PERBLACK F ORE ST where a colony of 
enslaved hobbits calculate the value of the 
function. The (valve) of the function comes 


back by return mail to replace the whole 
Kit and Kaboodle, 
Thus, if you had pA = 4. 


1 


=e ee 


. be 
Y= aS 
you would end vp with x= 4 EF Y= 5. 


ek 


Remind me- 
What was he 
trying to Say? 
wy 


Thus,an Expression 
consists of a combination of 


Liritometic Operators: 
Functions & Parentheses. 


CHI 
@dd ly enough, NO valves get assigned to 
variables by extraterrestrial forces, Whey get 
there because Domebody put them there, 
hopefully by intent. Fo avoid recycling 
somebody else's garbage, be SUAE to 
assign specific values to variables foefore) 
you try to use these variables in an expression, 
Variables can’t appear on the right side of an 
equa| Sign unless they were assigned valves 
earlier in the program! 


PROGRAM ?? 
What’s aA, dah, 
program ° 

don’t Know anything 
abovuta pro gram! 


IA masterful piece of understatement, 
my naive little polyhistor! At this point 
you could write an occasional, disjointed 
statement in Gompvtereze, but you couldn't 
actually solve a problem. 

Before you try to solve a problem vsing 
a computer, it’s only fair that I tell you 
a few things about what computers are like 

once Knew a fellow who called upa 
computer datin service, Whey tried tofix 
him up with an IBM machine. Needless to 
say, it didn’t work out, He was witty, hand- 
some, a dashing bon vivant raconteur — 
terribly debonair, Whe computer was dumb, 
unimaginative, totally lac king In creativity. 
He spoke fluent Francais— the machine only 
understood FORTRAN. Phe poor computer 
couldn’t follow his exuberant dialogue — 
it could comprehend only the simplest of 


grammatical constructions, 

(OF; course, it didw’t work out. For one 
thing, the computer was a lousy cha cha 
dancer, For another thing, there was a 
difference of religion. The computer 
thought it was God and he was upset about 
how they’d raise the children. 

o be fair, we should recognize that 
computers have many good points as well. 
rue, they are basically stupid. lowever, 
they have Great memories, Whey are 
incredibly fast, they are terrific t pists, 
and they make good money, last but not least, 
Zey vill follow inztruKshuns to zee letter! 


say a,dahya 
bureav can do 


all that ? 


oa 
How will you 
solve my problem, 
O Gomputer ? 


[ee 


How would you 
solve it by hand, 


\_ Murgatroyd? 
ee 


) hate to shatter your illusions, but using a 
computer to solve your problem doesn’t relieve you 
of the responsibility to think. Au contraire, 
mon frere! 

Because computers are so singularly literal- 
minded, you must be PAINSTAKINGLY Precise 
in telling them just whet you want them te do, 


The computer will do the RE work, 


Carrying ovt caleviations 4 Yle» -gedlu, 
thousands of times if necessary. Unfortunately, 
you need to fiqure out the (Sequence) of 
operations and you need to tell the computer 


how to make its decisions and how to figure ovt 
what te do next. 


Sulia Child, the French Chef, 
would call thata Kecipe for 
solving the problem, She 
would probably suggest a 
nice Ghablis to serve 
with the answer. 

Gom puterniks call it 


on ALGORITHM. 


A\\gorithms are a little 
hike Logarithms inthat neither one comes from 
the Babylonian word, “NOSHERAT “which 
Mesopotamian scholars have translated to 
mean "Vestal Virgins,” Ih] owever, that is where 
their similarity ends, 

Logarithms 5 AS You well Know, area 
Sin ing Qroup at M.L.'T. 

n Algorithm, onthe other hand 

Cnot that one- the Clean hand!) is a 
procedure for solving a problem, A\n important 
feature of any alaorithm is that it QUITS 
after a finite length of time. hike they Say, 
time is money ~ especially on a computer! $$$ 

Dn other words, an algorithm Is A Sequence of 
instvuctions which, _ 
comesto an end! St 


One can graphically picture an algorithm 
drawing a SAYFLOWCHARTAR.A 
flowchart isa sort of pictorial outline which 
Summarizes what happens when, and where 
to go to, should the calculations turn out one 
way or the other. The flowchart shows the 


GRILWD scheme of things, 


For a small program you don't necessarily 
need to actually putone down on paper, but 
doing so helps You to organize Your thoughts 
and be sure You havent overlooked some- 
thing vital / 

[LiKe most of US, You have proba) 
spent years ina dis organized fashion 
picking lint from Your belly button. See 
how much more efficiently this can be done 
when you systematize the process ina logica| 
fashion by drawing a Flowchart ? 


— 
a 


Look in Yile 
a 
No Pick it 
Ovt 
I 


‘Weave an 


Bar muff 


Toddle off fora 


lass of Gerito| or 
Other Liquid Refresher 


Using an expanded version of this flow- 
chart, in just a few menths 1 was able to 
produce a Shag Rug, two Tweed Jackets, 
three L'ies, and a Tea Cosy. 

Generally, in flowchart notation, the 
symbol (J) means ‘ Read cards,” the torn paper 
Sign | stands for “Write out on the printer” 
the rectangular blocks L_J are used for 


calcvlation segments, and the diamonds 


stand for decisions where there are forks 
in the road, ex = 

Don’t clutter up your Flowchart by trying 
to put every step of FORTRAN in the boxes, 
“Phe boxes are there to give a. CONCISE SUMmMary 
of the lo iC NECessary to solve the problem, 
Usually, the boxes don’t contain FORTRAN at all, 
but contain equations or English descriptions 
ot what's happeni ng. (In Yugoslavia they contain 
Serbo-Croatian descriptions, but the principle is 


the Same. ) Calculate Z using the 
>| DIMWIT Function, then 


We Vee 


When You Flowchart a problem, 
You included coy the logic however. Solving 
things by hand, Yov often take a lot for granted 
and make sim ple decisions unconscious |y in Your 
head. Whe computer won't take it upon itself to 
do anything, (lf you tell it to perform acalelation it 
will do so, but dom t stand around waiting for the 
answer. UF YoU want to Know the answer, you had 
better tell it youwant the answer! 

@nce your flowchart js complete and makes 


be sure 


sense when you go through it by hand, you can then 


implement the flowchart by means of an 


equivalent sequence of FORTRAN statements, 
to wit, a PROGRAM! 


How does it all get m 


the GOMPUTIN’ 
MACHINE ?? 


R a 56178 Grane aca woe wo0 71 12:73 74079 60 
ow 12 0 Yio Ys Soe le ee 
Row II \ a a ne — Ty 


[Each FORTRAN statement oets punched on 
an unstapled, unfolded, and unmutilated 
entity of hopefully uncorrugated cardboard, 


or your convenience, this material is 
generally provided in a form pretrimmed 


to the 


FaShionable &F Accugt omed 


size and shape. io wit? 


Figure 1 \ A Pep 
wr: AU/ 
Some left-handed deviates os HARVARD use 
cords of the form shown in Figure 2: 


Figure 2 


Figures three through Six show other commonly 
seen variants’ 


pos | DPD ee 


Kig.3 Fight = Figo Fig 
Left-handed  Cardas Right-handed Shredded 


card in card in Card upside 
backwards ees backwards oie 


Since the computer only reads the holes and 
pays no attention to the corner notch, for 35% 
Dr. Kavfman will let uou use either Fig.d cards or 
Fig.c cards, Just promise to Keep them right 
side up and not to tell your Mommy oy the 
Internal Revenue Service about the thirty-five cents, 


Trhereare MIGHT Y (80) (Count them/) 
columns ona punch card. Bach num bey, letter, 
or other symbol takes vp one column on the card 
and is represented by a coded combination of 
holes punched in the column. Where are |2 
rows that go across each card. For reasons 
Known only to Thomas Walsan: or God, the 
TOP row is Knownas the /2 edge and the bottom 
row is Known as the q edge. Y hos, one sometimes 
sees instructions, suchas the following alleged 
excer pt from Thomas Watson’s will, which 
specifies that when he dies, he is to be buried 
| edge face down,” 

ne) 


Statements get punched in Gol umns 7- 72 
of a card. Each statement goes ona separate card, 
and these cards stack up in sequence, IF a state- 
ment is too long for one card, it can be continued 
on following cards by punching a nonzero mark 
in column SIX on each of the extra cards. A 
letter G punched in column 4 of acard means 
“this card contains aGiomment for the bene fit of 
any humans who might want to read /t. “She 
computer lists Comments but ignores them when 
it thinks, (Each comment card must have a Gin 


$ 
Not Sherlock Holmes’ friend .ThomasWatson founded IBM. 


column 1, Since comments aren't statements, the 
other columns can be used any way you liKe, even 


column 6!) Gi; THEY MAKE Your PROGRAM CLEARER 
GATCH ON TO USING LOTS OF GOMMENTS 


Column 1 


The computer also Ignores) columns 73-80, Ya) 
you can write anything you want there. Go\umns 
cars ave used for Statement Numbers, which 
Will soon be cast like pearls before your very 


Cues, 


Lv typical FORTRAN statement might be 


| THETA=(2.* COS (PHI+PI/3,))/4.3 


(A NON-ZERO and NON- BLANK mark in 
co \umn SIX MEANS “Gontinved from 
preceding card.) 

Follow along and see if we agree on what this 
statement means. This wovld be evaluated b 
taking the Valve of the Variab e €PI} (PI is 
a vaviable-today you filled it with 3.14 but who 
Knows what you might put in tomorrow, maybe Pizza), 


dividing that number by the Constant 3, 
adding the value of PHI to the result, taking 
the Gosine of that number, doubling ics 
dividing that by 4.3, then multiplying the 
value of X by Duets raising that number 
to the ~2.5 power, and (qasp!) subtracting 
thatresult from the firs part of the answer. 
This Single Fina] Number would then 


Replace what ever garba e was currently stored 


as the contents of THETA. Right?? 


ae ae an, 
#8) (you snovld ave caidl 
(3 
S> 


al 


sree Lo Income laxSchedules Q,X,¢Z for other 
clarifying examples, 


Upon care fu\ scrutiny 

It 1s obvious to an 
intelligent observer that 
you have made a mistake 


There iS a serious bug in that Statement, 
You can te\l because there ave Two le ft 
parentheses but Three right parentheses! 


Just becavse you 
Say its FORTRAN 
doesn't make itso. 


Were there another left parenthesis, after 

the minus Sign, Os In 

ALI-( BOB eX) ¥-2,5 
then the statement wovld be evalvatedas 
described ear\ier. 

Suppose you had accidentally left off 
the right parenthesis after the ne as well 
as the left before the minus sign. Whe 
Statement would then read : 


: PTHETA=(2.%008 ( PHI+PI/3.))/4.3 
11:-5. 63 *X**- 2.5 


In this case, the computer would have 
raised the X to the -2.5 power, then 
multiplied that number by 5.63, and 
finally subtracted that result from 


2 Cosine (Phi + di) 
4.3 


p> that FORTRAN 
he's writing? 


(o.Tes Algebra, 
| you idiot! 
—_ 


Some mistakes, liKe two left parentheses and 
three right parentheses, the computer will 
catch and call to your attention ina helpful, if 
somewhat brusque and derogatory way. Other 
mistakes might still leave grammatically correct 


FORTRAN and really send you up the creek! 


Ww eatih JK 


ow To Ger Ofr Sue Sraaieyr 
\n> NaARRow Paty: 
ry 


Normally, statements are executed sequentially, 
one after the other. If Yyov want to jump from where 
you are in the pro ram to OO different spot, insert 


an Un GonDITIoNAL BRANGY in the program 


in the form | ( 


GO TO 50 
30 STUFF = TURK-EY 


it | 


50 TWAS =BRILIG + SLITHY 
ane 


T his Causes the computer to skip over 


the STUFF immediately Following the “GO TO“ 
and to pick vp execution with the statement 
numbered DO”, Any one to five digit integer 
punched in the first five columns of a statement 
cav\ be used as a statement number or * label," 
SO long as it doesn’t duplicate the number 
of some other statement, Further, my little 
chickadee, you can number any executable 
Statement your little heart © desires, (You 
can’t number comment cards. ) 


STATEMENT NUMBERS mere|y serve to identify 
places in the program. ahey need (net) Qe 
be sequentia\ or appear in any particular 
order, nh] owever, to avoid Floccilation' on the 
part of the readev (the Human Keader - not 
the Gard Reader) it is helpful if statement 
numbers appear in ascending numbers by tens, 
Then, afterthevghts can be neatly inserted 
withovt botching the orderly ascending 
arrangement and there’s less chance for 
confusion or duplicate numbers and theve’s less 
chance for confusion or duplicate numbers numbers, 


Only number those statements Yov plan to refer 


To branah, or 
Not to branch? 


Thatig the 
Question, 


i Blocci lation (flok”sa-la’shan) The deli rious 
search for Imaginary objects, 


COMPLI ONAL. 
BRANCHES 

allow You to branch to different 
Spots in the program, depending 
on the ovtcome of certain cal- 
culations. Yn other words, the 
program can decide what to 
do next, based ona test of 


its own ‘eeulian 


om (Aner 


nae ee allows Yov to branch 
to 1 of B locations depending 


on fos an expression is 


CD>Negative , 


nh ay 
Pulte, 


- 
g how vs an 


ARITHMETIC 


UF 
\f you will / 


Here's the form of a ty Pical Avithmetic If: 
IFCEXLAX) 20, 10, 30 


Between the parentheses is an Expression 
This expression will make You go someplace 
immediately! Df the expression is 


Negative ) 9° here 


Zero, 9o here 


Positive, go here 


In general, any sort of arithmetic expression 
‘can go between the parentheses, and any three 
statement numbers can go where I have 20, 
10,30. Just remember: the first statement 
number tells where to branch if the expression 
is (heqative); the second, if it is(Zerc); andthe 
third, if jt is Pesitive, Orhe same statement 
Number can be ve peated , if You want. For 


example, TECSUM4TOTAL) 100, 15, 100 


would cause a branch to 100, unless the SUM 
plus TOTAL came to zero, Then the program 
wovld branch to statement |5, 

[In a real-life Situation, you might see 
a confusing program like this sequence of 
excerpts: 


0 SPI Rd = BOARD / WALK + PARK XPLAGE ‘) 
IF (GRAFT ~SQRT(SPIRO)) 30,75; 60 


GOTO 50 Expression 
é GoLECT = 200% DOLARS 
SO PENSUN =850000.*xDOLARS 


t 
30 dal sasviay-aties 


fF the Expression has a yoitvelaie the 
program jumps to statement ¥ 60 

If the Expression has a Zero Value, the 
program passes the GO and reaches #75. 

Uf the valve of the Expression is Veg ks 
however, the program jumps to IMLS, it 
does not pass through the GO nor does it hit 
statement #757, 


Gouldja giveus 
dah..., another 
more tekni Kul 
xam pull ? 


Gertainly, my alodominous friend .Wou 
might say 
IF (P1-S@RT(VAR3)) 50, 25,25 
meaning : Jump to statement...” Ps \ 


—_ —P —_ —_ atl \ 
——=— - 


50 if tne square rootef VARS / 
1S bigger than PL... but te’ ) 


c5 “ig tne Sollsire eoules VARS : 
SB is equal to FY... but (Also) 
25 if PI is bigger than VVaAR3 
See how easy arithmetic ifs are. There's no 


reason not to learn all about them. (You just did) 
But there are ee Ways to play the same 


tune, SREP 


ee Za 
AThe [Logical Gu 


Dainaweuld be to Skip Ya 
all this DRIVEL and a 


turn to the -& 
CARTOON 8! Ce 
SZ . 
/ 


\m I glad Yov said drivel arhat reminds 
me. One can also write branch statements 
which are followed if a LOGICAL 
EXPRESSION is TRUE and \gnored if it is 
FALSE. The computer works on the principle: 

"Yo thine own calevlations be TRVE, 

And it must follow, as the night the dau, 

T hov canst not then be FALSE to any man, 
“Vulgate Edition ot Hamlet 


Be thatas it ma , and it seems a dubious 
hypothesis at best, one cov\da Sau in FORTRAN: 


| IF CHW.EQ. TRIVIA. AND. 
1/TEXT.EG, DULL) G0 T0 210 


210) (DATE= SEXY XX¥NTH 


Gontrary to Common belief, this is NOT 
an Unconditional Branch but IS O 


Logical Is 


ae 
Dt means” UF Both 


the real variable "HW equals 
the real variable TRIVIA’ 


the real variable “EEX 1 wae 
the same valve as the variable“ DULL 


THEN sKip tothe DATE 
OTHERWISE, plow on to the statement 
following the Logical Lf, 


Orhe stuff in parentheses in the Logica | \F 
IF ( ~~~) ~~ 


isa Logi calE xpression which cn either be 


or “C Logical Expressions 


Unside a Logical Expression, pairs of 
Avithmetic E x pressions are compared with 
one another by means of RELATIONAL 
OPERATORS. (lRemember-arithmetic 
expressions when evaluated, reduce to a single 
number, Whus , Relational Operators ave 


operators that sooneror \ater end up comparing 
Lwo numbers 1) 


Dah... that’s a 
nice simple - 

minded way to 
think of it! 


CrPhese RELATIONAL. OPERATORS 


= Oo. (onl) 
NE. (Not Equal) 
Gis (Greater Than) 
.GE. (Greater Than Or E qual To) 
LT. (less Than) 
LE. (Less Than Or Equal To) 


TDG>Wart Gore) 


Notice that the decimal points are Part of 


the operators themselves. That's how 
the computer Ky ows Aa Ls iS an operator 
meaning Less Than “instead of a variable 


meani ng ‘L eroy’s Tabylation.” 
Df you say, (nes Leroy? 
Re G TX 


YoU are asking the musical question 


eo ae aad aie 


fhe whole expression Y.GT.X” has the 
Logical Value of either being True or False. 


HES YOv Say 5 
ALPHA.Fa38.2 227 
SY e® 8 e 
hogtcal Expression 

and the answer is “Yup!” then the whole 

. ny Vv 
expression 1S trve. If the answer is Nope, then 
the expression is false. It’s aS simple as that. 


What's the 
difference between 
ALPHA .EQ. 38.2 
and 
ALPHA = 38.2 


All, THE DIFFERENGE istie WoRLD! 


On layman's terms, the answer to your excellent 
technica\ question 1s “A pproximatel the same 
difference as there is between an Elephant 
and a Prune.” Saying 


ALPHA = 38.2 


means “take the real constant 38.2 and 
store it in the location named ALPHA. 
DF Some thing is alread there, throw it 
away and put 38,2 in its place,” 
Saying 
ALPHA.EQ.382 


means ‘look at the number now stored in the 
location named ALPHA ond Gompare it 
with the number 38.2, Df they are equal 
SAY ‘Tyve” and if they differ, report back 
with “False.” P. S.- Whatever you do, DON’T 
change ALPHA — just look at jt!” 


ea 


eS, 


The power of Logical Expressions is further 
multi plied loy the 


sian .AND, 
P on OR. 
pérators NOT. 


and by the two LOGICAL GONSTANTS 


Be Py, A ae nei ( if, 


LOGICAL AND says: 


“T€ the stuff tothe Left 
of the CAND.) is True... 


oo 
baaagenaat 


AXTN LL 


the stuff to the Right 
‘sALSO True, <e ? 


then that whole piece of logical rn is. TRUE, 
in valve, Otherwise, its value is ,FALSE, 


[EF xam ple! 
OUM.EQ. TOTAL .AND TOTAL.NE, PI *2, 
NN ed VV ——_—_—~ LY 
This whole thing is Trve only if beth sides are Trve, 


You can’t paddle 
your canoe withovta 


LOGICAL OR f 


NSA ACK 


Logi cal Or means “ {Ff the term to the Left of 
the .OR, is True and the termto the Right 
is Also Trve, or else, if Bither term is 
individually Trve, then that whole chunk 
of logical expression is . LRUE. in value. 
Otherwise it js FALSE, .” 

Example: 


AGT 9.2.0R.Y.NE.Z 


If this Vis Teve OR thisis True, then the 
whole works is TRUE. in logi ca\ value. 


are both Binary operators. Whey each take 
a L ogi cal Rxpression (i.e, Something that is 


either Vrve or False ) on their Left and 
another Legical Expression on their Right. 
Bach expression indivi dvally is €ithey Trve 
Or False and the Binary operator then 
logica}| combines the two legica| values te 


end up with a single final answer of Trve 
or False, 


e Or, 5 not like that, 
NOT. is OO Unary operator which wee lA 


Y 


Suppose you had 
xX = 4,3 
Then the expression 


XxX 16) HO 
wovld have the valve of .FALSE. » Since 
X isw’t less than 3, ithe ex pression 


NOT.CX 4-79.50) 
woulda have the valve of TRUE. therefore. 
(IK the ex pression in the parentheses had been 
True, then the overall expression would be 
FALSE, in logical valve , 

S ome arammar ans land even an occasiona\ 
grand parian ) find .NOT.s confusing. B 
didd ling around with the Relational Operators 
nan expression, you can usvally write it in 
a cleaner form Withovt Using .NOT.s. Do 
whatever will make your program Cleavest, 
So you won’t Screw up. Don’t avoid skipping 
the use of NOT. S when You find not 
penys their negative aspects dou)\ 
contradictory to theiv unhelpfulness. 
(Only vs NOT wher you need to.) Some 
people aet so screwed UP they write programs 
with Square a OC. S. Yust don’t tell De Morgan, 
who you shovlan t confuse with Henry Morg an not 


that you would knew who he is, since he 
went off the air when you were still in 
swaddaling clothes. 

(zzee 


iY. 


RAE aVWhilt, racket Leg 
Expressions, one might Sau for example ° 
(OAL @ OR.Z.LT. (X+¥Y)). AND, (B-37,)*% X 
1} EQ SQRT (Z-X) 

SPhis titubant, peri phrastic though 
pers picuous orotunda and \uculent 


logodaedaly would be . TRUE. in value if 


Vex 


t This has gottabe trve and ene er 
gotta be trvé for the whole thing tobe true, 


SoA logica | If 
is just another 
way to makea 
branch, right? 


Wrong . Wrovev hog' cal branches 
barely use hrf the power of Logical Tfs/ 
Woof what we ave abovt to Cover Pex, You will 
be able to emBavkK on many diferent varieties 
of Legi ca\ fs. Al Pe, the key to trem all is the 
propev vse of logical EXpress] ons. ie a ) 

rhe genera) form for these dog gone Lf 
Statements is‘ 

TY (Logical E xpression) Domething 
At the valve of the expression is .VRUE. , the 
Domething gets done. Ifthe expression is ,P:ALSE, 


then the Domething gets iqnered and the 
Computer chugs on to the next statement in 


the program, 
Whe (Something can be almest any FORTRAN 


statement . OU CAN Sa 


TF (X.GL.Y-ALPHA) GO To 400 


o IF(PI.LT.BETA,AWD.MNEK) X=Y-PI 
(Only if this expressionis TRUE., does i Eo eee 
get executed.) 
Ina Loaical ‘ca the Something Gan’t be another 
If statement, however! 


[But You Could Say 
IF (JJEQ.500) READ (8,10) XY. 7 


Maybe You 
' can KEAD X,Y, Z, 
but Nobody ever 


taught ME to 
READ orWritE ! 


Unconfutabl 5 the veracity of the pubescent 
damsel?s trenchant alleaation necessitates 
expedient rectification of the omission. 

Gonstants whose values might vary each 


time a program is used ave generally punched 
oY DATA C ARDS vather than being 
Specified in assignment statements. This way, 
the SAME derovaged (or, Ifyou eve fer, 


unlepi dop teran ) Pprogvam Can be used from 


yun to run, with only the DaTA GARDS being 
altered. By branching back and reading in new 
data, you can also have the same proavam solve 
Severa \ problems durin the Same Yun. (That’s 
how the Phone company Sends ovt bills APhey don't 
re write theiv billing preqram for each individual 
customer.) Also When You chan C numbers Using 
data cards instead of by rewriting the actua] 
program, the com puter can save time and 
because it doesn’t need to recheck for evrovs in 
oi program You Say is OK. Arhe other me ythe progvam 
\tself wovla have to be changed whenever a data 
valve changed , so the computer world need to 
recheck the whole prog ram for bugs. hat, mon 
frére, is POOK for 28 geod ana sufficient 


YEAaSONS, 


Aro obtainthis list of reasons, senda stamped, 
self - addressed envelope with $12.36 in GASH to: 


De R.E.KaurmMan 
WASHING TON, D.C.20052 
Te 


Constants punched onthe Data Garps are 


stusted into the appropriate Storage locations by 
a statement of the form 


OfferVoid Where 
Prohibited By Law 


READ(8,20) ALPHA,X, MOM, OIVEY 
ANE ey EEE 


Input List 
Statement Number of 
A gsocia ted ¥ ORM4 > 


‘Wom ber that Jimmy fale tah 
assigned to 


ir he vavial\les whose values are tobe read in 
from cards ave \isted (in the order) they wil] 
Appear on the data cavas. aphe second inte 
in the parentheses on the READ (~, 20) ~~is 
the number of a FORMAT STATEMENT 
which tells WYEre and ‘now the numbers will 
be punched onthe damn cavds. 

Arhe LORMAT statement can 90 (Gnsjuhers 
in the Pregram because it doesn’t do an thin 
Dt just Says how things are to be done. he 


format associated with the Read Statement 
just: given mi ght be ; 


Not floormat, 
you dummy. Lasked 
for the FORMAT! 


Ls Iwas saying, the Pormat might be: 
20 FORMAT (2F10.0, 110/F 10.0) 
———S 


These are Fiela Specifications 


Bach element in the In vt hist has a 


corres ponding TELD eprdlFICAT 1 ON in the 
FORMAT. 


ALPHA, X, MOM, 0OIVEY 


eodub 


Cy field specification that 18 te be used 
severa\ times in succession can be 

indicated by preceding the spec by 
a Repeat Constant, T hus, this format 


is equivalent to 


F10.0,F10.0, 110 / F10.0 


Wn the Format, Field spec} $i cations are 
separated either b Comma. 


s or bu slashes. 
Dna Read Sta ment, list elements) are 
separated by commas, 


Great. But, dah, 
what IS a 
Field Specification 
aad 


LA 1 field is used for 


¥ aes Point Numbers. 


A\n Drea Is used for 


[rteger numbers, 
Den?t Science Simple Wonderful] ? 


Furthermore, OO speciticati on such as 


F 10.4 


means 


"thenext Ten columns are to 
Contain bhis lis € element with four 
pla ces atter the decimal point,” sf 


On of Floating Yoint Numbers, the 
computer is smart enovgh to actually read 
the cards in <=> Correctly (as 
contrasted with “ incorrectly”) even if they 
are punched with a Aiffevent number of 
places atter the decimal point he Aecima\ 
point punched on the data card OVERRIDES 
the Field Specification. DF ne decimal) point is 
punched onthe dake card, however, then the 
Fiela § pec takes contro| and Says, 


Vhis feather brain was 
too lazy or cabbage headed 
to punch a decima) Po} nt, 
or else he Aaidan't have room 
for one on his data card. 1 
Shovla throw him out, but 
instead I'l) kindly assume 
he intended to have one 
fovr columns befove the 
end of the fiela , since 
that's what he saia inthe spec 


F10.4 


Only Fools, Acrobats, and M.1.T. 
Fresh people 


resnmen deliberately leave out the Aecimal 
point Withovt compelling Dve Gavse f 


‘rhere 1S one restr ction on punching Floating 
point numbers on data cards, That 1S) 1F Your 


spec said F920. 6 
then you ve gotta get the number into thease 


20 columns ) even ¢ You don't have _ line 
UP tne decima) point S1X places before the end. 


sagogic epexegetic Comments of a 
para enetic nature now So hy pragoaically 
expressed , We come to 


UNE GEIR FIELD 
GS PEGIEIGAY IONS 
[\ Fiela Specitied bu 
L110 

means “The next ten columns onthe dat card 
will contain an LIN WEEER value, vight- 
justified me For those Who flunked English, 
that means ‘shoved over to the right.” 


You probabl wonder 
where the data cards 
go in your cava deck. 


Bach time You Come toa 


READ (8, ~—— 
you start a new data card. Dats, cards ave 
stacked sequentially in the ovaer they Will he 
read - Arhey o into the cava veader after 
the FORTRA program ana ofter whatever 
Contve| cards are needed to start the 
Program executin : Data numbers can be 
punched in any or al of the 80 columns 
of a cava, since the computer reads whichever 
columns you ask it to vead! 

AA SZ4,977 ina FORMAT means either 
Hell?s Angels carved up your deck, or else 
. okip tothe next data cava.” Arhus, if you 
wanted to skip past two data. Cards, you could 
nave 77 


m your Format. Baz slash brings yov to the 


start of a new card 


D+ Yov want to skip ovey ov \qanove the 
next few columns On aA data Cara, Yyov Can 
use an X field specification, like 

x 
which means “skip past the next Krteen 
columns on this card.” 


Grhos, te make a long stovy longer, 

READ (8,20) ALPHA, X,MOM, olvey 
20:: FORMAT (2F10.0,110/F10.0) 
Means: 

WS When the computer comes to the Read 
statement, it starts a new data card. 

OS The first list variable to be vead in is ALPHA, 
Ww Alpho’s value was punched as a floating 
point constant in the Hirst Ten Golumns. 

S% Kis the next list variable. X?’s value will 
be found in the Ten Golumns following ALPHA, 
he MOM isan Integer whose valve will be found 
Right- Justified inthe next Ten Columns. 
So The Slash means “forget about the restof 
this card and move tothe next data card.” 
& OIVEY”’ s valve is punched inthe first 


Ten Co\umns on this New cava : 


dah... Why Always 
Ten? Why never 
Tree, or Toyteen? 


On INPUT, IT like to standardize on 
fields that take vp a multiple of Ten columns. 
Arhis Is & good habit to acquive, as the nun 
said, because it facilitates kegPinching 
and cheSh'"9 datacards, A glance at the 
cards Shows if the numbers are in the right fields, 

INeypunching 16 Aone ina noisu smoke filled 
room inclose proximity toa lot of rude systems 
programmers who never heard of Arvid Extra 
Dry. Thus, It?s good tobe able te doit quickly. 
Use short shallow breaths (try not to inhale) 
and ask one of the systems people to show you 
how tovsea DRUM CARD. 

arhis drum card fits into the keypunch. t's 
like setting the tabvlator ona typewriter, but 
even move versatile. DLcan make the key punch 
avtomatically duplicate or Skip cevtain columns, 
or even shift from alphabetic to numeric mode. 

ZA FORTRAN drum cava can save you hours 
in punching frograms. irhe one 1 use as L peck 
oway at the key punch avtomati cally shifts to numeric 


Be Sure to downshift for G’s in Column 1 for Corment 
cards. For continuations yskip to 7 and backspace 1. 

[If You follow my advice (learned the hard way) 
and always vse fields of Ten for data , You can make 
a dvum card for DATA as follows: 


RSpace, » SSpace, 
Followed bY ». Falowed by ; 
@ Ampersands : q Ampersands « 


more information on DEVUM CARDS, 
see seque\' to this book, Which 1 plan never to 
write, due to tremendous populav demand. 
({H) owevey, Lam negotiating a contract for the 
TV, and movie rights, with Donny & Gher 
playing the part of the Uncenditional Branch 
and . ames Garner OS Column Six.) 


rs 
Bongos ,onares, Kettles, and Keypunches. 


Perforate dF orms 


How to perform 
Input / Output 


with DELICACY 
and ¢ tyl e 


Note: Originally, the above figure was to have been 
that of a well-known ex - politician, seated ina 
wooden privy. This was censored,so0 as not 


to offend any members of the Birch John Society. 
vy 


Outpt 1S almost like input bot its not 
input, it’s output. Wo correlate the ovtput 
with the input and te veri Fy that the input was 
eut in correctly, IVs a Qood idea to output the 
input along with the out put. 


laving clarified that point, observe that 
one could now output the ineut and the output 
by Saying gus 
WRITE (5,30) ALPHA ,X,MOM, OIVEY ANSWR 
ee 
Ovtput List 


Format that tells how se make it 
look pretty on the page. 


The LOGIGAL UNIT NUMBER 
assigned to ovr Line Lrinter by Qn 
Inebriated Systems Lrogrammer* 

Srhe associated Format might be: 

30: {FORMAT (‘I", 5X, ‘ALPHA b=b’, F104 / 
‘li "bXb=b’ ,F 10.4, "bANDbMOMM Yb=b’, 
2/15 ///"bOIVEY bANDb ANSWERbARED 
37 IX, 2 (F 10.4, 5X)) 

L\nagogic though it mig Wt seem, one quickly 

* Unless Yovr Systems Programmer Was at Duffy's Tavern 

with ovr programmer, you'd better ask him the unit 

numbers Sov your printer and card reader! 


realizes thatthe FORMAT STATEMENT 
has Real ana Integer Field Specifications 
bearing a ON FE to OWE correspondence with 
the elements in the (oul put list} @ne also quickly 
learns that You can't key punch the Format the 
Wa. I wrote it because (haha!) there isn't 
any small b on the key punch! Devious thovah 
i an y T didwt do that just for the fun of 
inflicting misery. Trhal’s Simply an added benefit. 
he reason 1 did it is because small ‘b” is 
Computer jargon for “leave this Space blank.” 
Hverybod y in the Know Knows b stands forblark. 
Formats are one of the few places where blanks make 
adifference ,so when they are im portant (and not 
obvious) indicate them with a‘b” (This js mainly 
important if some one else is goin to key punch 
from your hana written Scraw | ) awen, back 
to the format, 
On output , the F ORMAT Speci Ges 
ERAGTILY how the numbers will appear. /\ 
specification oS F104 wi)\) cause a tloating 
point number to be Leopped OFF to precisely 
four figures after the decimal] point, rhe 
resulting numbey would then be printed , 
right-justified in the ten columns. 


iS ns) Ae the com puter, the entire 
number is still there, even 
thovgh yov only chose to print 
four decimal places. 


LA SSash (///) on OUTTUT, has the 
effect of ski pping to a new line, ( Remember, 
on input it skipped toa new tara.) 

Wndubi tably ) You have noted by Now the 
presence of a new Kind of beastie in the 
Format #30, Dt is, to wit, A LITERAL 
DATA FIELD SPEC,TRA LAV (The 
“tra la" is optional,) A\ny A\ phanumeric 


(Numbers and Letters to you, marplot 
Characters enclosed in a paty of Single Apostrophes 


are printed ovt just as they are written in the 
format! (Phat aoes for blanks too!) 
NOTE! This apostrophe is not a 


*TH15 ble ALPRANUMERICGb DATA’ 
(Phe avostyophes themselves don’t get witten out.) 


Thin gs 9° 
better with 
Keala 


LAs your evealy ptus- eating arborea 
marsupial elder could tell You, that means 


F10.4 5X, F104, 5X 


AM ny Tield specification canbe preceded 
bya repeat byavepeat & vepeak constant asin 
AF 10.4 
which simply MEANS 


F10.4, F10.4, F10.4 


For instance, 13) ('X") would QiVe you a 
whole line of asterisks, 
Of course, the 5X means skip over S spaces, 


YPov can also throw paventhéses around a 
group of field specifications and put a 
repeat constant in front of the whole works. 
iis you really have a thing for asterisks, 

a format like 


GO (1X, 131 ("*")/) 
wov \d give youa whole page full. 


| CONTROL 


oop Vi 
oy 


Wa - S ~ 
RS >. di Ch SLA 


The First Golumn on the \ine printer 


is Non Printing and is used for CARRIAGE 
CONTROL. BPA LWAY Saxe start 
each \ine of printing with a blank ora Geliverate) 


cavy\age contr o| Sy mvo\ oY You CAN yun through a. 


amount of paper damn fast. 


Vf you try to print a I atthe start ofa 
new line, the 1 will be wedpbtesteae and the 
printer wj II skip to the top ‘of a New page. 
NF you want to starta new P29e@ y that’s how 
to get there, [But if, by acci dent, your 
program is in an infinite loop, printing a 
Lat the start of a \ine each time throvgh 
the loop, and never quitting, Smokey the 
Bear will personally teach you how it feels 
to be recycled, 


) 
i> 
id? 
BEGINNING PROGRAMMER —— Y 
STARTED EACH [Z_) ws, | 
ee wiTH A cK a iS VG, 
sly NY, WY why Sad Pn NDEXTRT gM ND 


EN blank as the Carriage contro\. Causes 
the printer to advance normally y one line, 
arhus, YOv generally try to arrange things 
So that the first charactey Yov Print oneach 
line is a noncharactey, that 18, Zilch ,a blank. 

S trietly speaking , I lied to you before 
when 1 said a Slash in a format caused you 
to go bo anew line. Phat’s usually the end 


result, but only because the first character 
printed after the slash in the format is 
usvally a blank space. [It’s the blank that 
Causes You te advance to a new line, not the 
slash. 


So what does the slash really do? It 
Causes the printer to get ready back in 
column | and wait Sora new line of characters. 
Un other words, ee ctively it just Moves You 
back to the left margin of the paper. 
Whatever You print next deter mines whether the 
paper will Move Up one line Or ONE page or 
not at all/ 

LA zevoas the Carriage contro| character 
will aavance the apPpey two lines before 
printing. IN ota Japanese Zero You 
nur. , A ZeYO Zero | his gives the Same effect. 
aS YoU wov la get from an extra pair of slashes 
With a blank cayviare control, Srhat 1S, 

H//"»" is equivalent to /'0" 


Is there no end 


F ina lly, for the co\\ectors of arcane tidbits 
of trivia, 0. plus SIQN..., + ) when used in the 


\Steolumn On the line printer, serves to 


SU 7 ESS the 

ts PAP EZR ADVANGE, 
thereby making it possible to OWERPR MNF 
‘his trick can be useful in certain obscure 


Printer Plotting Programs, or for Making 


DARK HEADINGS _ 
CWO Wi 


ASD 
Sp < e ys Sag) ~ 
~ a Nei i~> 


(by over printin the same information 3 or 4 times), 
or if Youarea MEAN AND SPITEFUL SOB, 
You Can vse overprinting to Saw through the paper, 
the ribbon, and to completely ivvitate and scvew 
vp the Gomputer Staff. Be caresvu\ how You vse 


+ Gayriage Gontrol, OY YoU Can earn an 


unintentional reputation as being mean 


Spite ful, oy just plain STUPI D. 


p) 


Whenever you start a new Write Statement, 
you are automatically at the left edge of the 
paper. Where you ao from theve is up to you. 

Uhus, the otted of 

WRITE (5,30) ALPRA,X,MOM, OIVEY,ANSWR 
1 3 3 
30: {FORMAT (‘I", 5X, "ALPHA b=’, F104 / 
H'bX b=b’,F10.4, "SANDbMOMMYb=b", 


2;15///"b OIVEY bANDLANSWERDARED 
Si 1X, 2 (F 10.4, 5X)) 


i } " PPx ie 
would be bo print a page in vs ‘ 
the layouts / 


l 


Uncidentally, most. line printevs have 132 
columns (131 print positions plus the \* carriage 
Conve | co\umn.) 3 ome have on) \20 columns, 
however, and still other chea pies have only 
80 columns. arhen of course , there’s Mm 
adding machine butits printer is beside the Point, 
Since it doesnt speak FORTRAN. 


Prithee, Stay this 
torture! Ver; ly, 
this rankest folio 
| : of NOK\OvS puns 
Vas & is dreadful 

"4 repetitious 


True, bard of Ayon Galling. Many things 
in computing ar€é repetitious. Io help you bear 
that psycholegical burden, God created the 


AG Statement 


Oe Ge ; 
Guebat pee Programs often have 


Wz alo OP, 


4 
> ve GO for one 
ARS 3. on 


Yar 


Yeason or another, 
mainly another. 


W hat do you mean 


“AL ooping Structure” 
ere 


ed 


Vag 


ee 
Frequently, Programs will ie SF om 

repeat a cer tain qroup of ca\ cU lations 

over and over for some given number oS 8 

of times, AA\ typical looping 

Pi ece ot Potential Frogvam 

might go \i Ke this: 


COUNTER 
D) = (~~ Wrar vais) 


<. 
lO BLOCK OF CALCULATIONS Se 
To BE REPEATEP \op 
TWENTY -FIVE TIMES > la 
x 
T= 1+ WcremenD ( 


-IF(I.LE.25) GOTO 10 


stuff to do next 
TEST) 


ahi $ examole showsa block of 
program which is to be repeated 
with 

[=\ ; T= roe T= TEA respeqzeurs fete, tee) Le 22,1223 =24 1525 
soa G ounting L oop has been 
set up. 


Dn this particular case? 
ff tte Counter» 's not FORMICA _ 
bubis the INTEGER VARIABLE) 


I=| ds The Uncrement) each time through 
\0 ~~ the loop is¥1) in this example. 
——~ —- Fe Whe first time throvah the loop, 
sf the {Initial Value) is als& 1), 
we Wust before the end of the loops 
T=1+\ the counter is Tneremented), 
AF (1.15.2) oy) AN (Test) is made at tne end 
:Q0 To 10 
of the loop. 

SM A (Final Valvey is checked. 0S 
the incremented counter does 
not exceed the Final Value, you 
plow back through the loop. 


OF course, there are many other Ways Yov could 
set up a loop ina pregvam |No matter WOW you set 
vt up stheugh, you would still need each of those 
things with (Clouds) around them: a. CounterVariable 
an \incremenl,, an Initial Value), a vanda Test), 
The partiolar form of loop shown in the preceeding 
example is especially important though, because it shows 
how a DO LOOP works , in essence. 

Dn its most geneva form,a looping proavam looks like Pas: 
ENTE 


Initialize Ry I weREMENT Tesr 
COUNTER NO x" £078 y, eave Wir 


| oe 
Fro avoid writing 
a\\ this a 
te\\- tale technical 
tedium , 

Ay) 0 sak cauetl 


was invented. 
ee A 


Dt combines all the afere mentioned 


elements into statement. Aro wit: 


se 
WW A 
ee aw eo & 
Rog ray “¥ a 
ar <f > 
DO 15 I= | 
te) 
ra 
a5 
a 
3 
SIS 


Trhis has the effect that everything 


throvah éy including State ment #15 
gets _= just as before. Phe Range Limit 
is simoel y the number of the last se ieneah 


in the loop. ‘rhe loop MUST end on an 


execvtab\e statement ~ one like an Assignment 
Statement ora Read | Write. You GAN’T end on 


a Format, Go To, If, or another Do! ¥ 
Wh ve DO me like you do do Dd: 
CE SEA 
Yo help avoid ending Oon-a bummey, CONTINVE is @ 
avs, pen executable statement thatcan 
C3 be vsed as an end to aYOLoop 


or whenever else You wanta 


numbered spot in a program as a place to 
branch to or for some other devious reason. 


Keep your eye on it 
ana \et me run that 
past You once again 


DO 30 K=MSTRT, ND, ING 


INTEGER VAR/ ABLES 
OR CONSTANTS ~- No EXPRESSIONS / 


(78) 


sim ply means; 
LS Set the counter (2 to the integer 
valve MSTRLI 7 


Sa Do everything trom here down throvg h 
ee including slatemenl #350 


EYP Irevement K by the anoont INO) 


Gheck tosee if K exceeds ND) . 


Lf Net, go back iia 80, leave the 
up to the top & loop and 92on tothe 
the loopandy "7 statement after*30. 


rtrw 


Notice that the statement at the Range 
Limit js inside the loop: (Statement # 30 in this 
case) Statement 30 can be any execvtable 
statement , such as 

30 ALPHA = SQRT (BETA) 
ov 30 Y=Kr*4. 
or 30 LAMDA=K+MOM 


The Law ean SAUS statement 30 can’t be 
any sort of a branch statement, a DO, OY OW 
FORMAT. Whese are considered nonexecutable- 


they don’t cause anu actual number evunching. 


h¢the place you want to end 
Yovy loop is smack dab on top 
of a nonexecutable statement, 
then You Can get exactly the 
Same effect by throwing in a 
CONTINUE; 


ae 


25 FORMAT (~~ 


30 GON TINUE i, 
You cantenda here oY here but youcan stophere, 


athe only efect £ the GONTINVE its to give you 

a legal Place to end Yor dO loop so that the IF 
statement will be inside the loop and not outside 
Ot really 18 just a Dummy do-nothin statement 
on which yov can hangaa statement num ber, Eve 
thovah it does nothing , the computer considers it 
an executable statement , So jt doesn't Kick you 


ovt on your 


You can always enda loop ona CONTINUE. 


(Continued onnext page) 


ar hos, ina program, you might see 
Something like this: 


£8) The 
( Be ociety-Promusicg 
€ 


The Sch leptet mD 


C 


Ina FORTRAN program, onthe other hand, 
you mi ght see; 


_. 


DO 50 K=5,9, MILTY 
TOP= BANANA **K 
(. SILY = BERLE 
50 BOTOM = SQRT (BERLE) 
= Pe PAD 
AL the start , K is set to 5 . You od down through 
ti\\ you reach statement 50.4fhat's the bottom 
of the Berle. K is then incremented by 2 which 
brings itto 7. Since 7 does nok exceed 4, Yov go 
hack, UP to the statement ™ TOP= : and meandev 
down throvah Using K=7. Once again you reach 


the ‘BOTOM y (7) dont dave make a pun this time 
or You might Berle me in oil.) i\fter statement 
5O is execvted , K is incremented once again 
to 9. Since % does not exceed the upper \imit 
(which is q), you go back vp to the top. 
Once Yov finish with statement 50, K is 
incremented once again , bringin jt over the 
upper limit, (You can velax- Mere was no corny 
joke this time.) Since the upper limit is exceeded 
you Exit) From the loop and go on to “Joke =” 
and the rest of the program. 

You are probably wondering abovt the JOKE=,” 
Well, the “JOKE” is that when yov ve been 
Uhyovgh the loop with K 4, you're atthe dog gore 


sia ‘You have te handit to Dr. K. 
He has nerve! 


They ll probably never 
let him back in the 
faculty club! 


I hear the painted 
ovthis name on his 
office doov when this 
book was published. 


CC 


dah... What if 


the, vh, increment 


is 1° 


Funn You should ask! 
That?s a specia\ case. 
DE the increment ina DO LOOT is Unity, 
it can be omitted from the DO 


statement anda iwill be 
understood, 


( Understooca ? I don a understand ¢ 


INo, but the computer Aoes. [For instance? 
DO 10 MUSH=30,1QUIT 

means “do everything from here down through 
and including statement #10, Start with 
MUSH = 30> Weep looping through those 
statements , incrementing MU $i by 1 each 
time throvah. top) when you've been 
throvgh the loop with MUSH equal tolQuit, 
Then, merrily proceed onwards to the part of 
the program ‘following statement #]o, 


Don't Trinker with try ing to chanae 

valve of the DO COUNTER, the I ete 
VALVE, the FINAL VALUE, or the - 
INCREMENT while you ave IN a DO Loop! 
You can (ose) these values on the Sp 

side of an assignment statement, 

or in any arithmetic expression, but 
change them! You CaN Diddle with them befove 
ov after the loop, but not during ! 


[A Novma\ Exit (or egr ess) from a DO 
is One in which the sustem bun ws You ovtfor 
reaching the final test valve, On this case, once 
You are out of the loop, the DO COUNTER’, 
Value is UNDEFINED, so don’t try to vse it | 
Pon texpect it to have its s last valve, ov that 
valve plus the increment. Ne kngws what 
value lurks within, Gear piereniae tie 


YT) 


Vie, 
You canalso of a DO loop because 


ofa branch statement within the loo P- Nt 
may be that you perform a test of your own 
inside the loop and decide not to wait for 
the loop to vun to completion. aThis is perfectly 
0.K. even thovgh its called an Abnormal) Bxit, 


~ Normal Exit 
3 sales a Va 
, Oy) and his 


brother, Abnormal 


After an ABNORMAL EXIT, the valve of 
the DO Counter will Still) Keep its last valve, 
SO You Can vefer te that value outside the 
loop. 

‘Only Enter 2 DO \cop at the very beginning y by 
exeavting the DO statement. Once you ave inside 
the loop, Yov can jymp around within the loop to your 
heart's content. But don't jump from ovtside a Do 
loop to inside tne loop Unless Your insu rance Prenivns 
are paid up. IL is Risky Business! 

Whi ch brings US to 


THE NESTED Do: 


Hrhe < Nested DO ( domesticus dears) Sa 
beavtiful thing to behold. Dktis of the form: 


DO 100 NEX T= |,1000 

DO 50 M=ITER,J3 

po’ b¢ LOOP=M,J 
Gus 

DO 50 INT=25,50,2 


(_ an 

aes 

C al vi MUD =1, KICK 
80 


PS § 
a 


Orhis isa Legal Species protected by the 
National FORTRAN, Fish, and Wildlife Naency. 
Qt 1s legal from September 1 through the last day 
of Avaust. The yeason it 1S le a| 1s because Cach 


loop is (Completely > Contained within the next 


larger loop, LAn iNega) species would have 
loops that cut across each othey, asin: 


1 
\ 
V\ Se 
Arhis is Tilegal . be cause these} two loops jum P 

from a point Cutside a loop toa point in the 


ofa DO without Qeing in via the DO 


statement at the top of the loop. 
A 


rf af Dt is O.K. to end several DO loops 
sx >» onthe same statement as in 

DO 50... DO 50... onthe preceding 

page. hose loops ave considered to be nested! 


Dwell onthe 
je following Points ! 


W Phe LVCREMENT, STARTING VALUE 


€f FINAL VALUE of a DO statement 
must all be Positive Integer Variables or 
Gonstants. 
7 You Gan’ have Ich, Teh f 
DO 50 I=710,10 
or DO 30 J= 0,20,2 
Zero isn't quite 
positive enovah. 
or DOGS K=M+*2, 100 
Expressions are 
egal here! 
or DOIOX= -PI, PI 
X and Flare Reals, not integers 


and -FTI is negative 
or DO 15, I=1,10 


LT just wondered ! Gov Were paying 
attention. No comma. belongs here 


Wy Whe STOPPING VALUE’ can 
be LESS than the STARTING VALUR 


You Can’ ty have This cane 


be vegalive 
DOIS LL=9,3,-1_ 7 


This must be 
at least 7 | 
or KkK=2 
DO 40 ITER=5, KK 
LLY can't be fess than 
the starting va lve Which 
oe 


“rhe stopping value Can equal the 
starting valve.Arhis brings up the 
next point! 


¥ W henevev Yov execute a DO statement, you 
ALWAYS Fall throu h the loop art least 


Once ,y Ltegardless of t Final value. 
orhis is because the Test associated 
with the DO takes place atk the 


devriére, ov (if you || pardon the 
term) the End of the loop. Df 


you want to execute the loop 


Zero times under some situations, 
You need to prog Yam jn youy own 
test to skip past the laop. 
(Maybe I szia that wrang eC 
should have said “If You dGon't 
want to execute the loop Zero 
times...” or is it ‘IF yov wayrt 
to unexecute the loop...’ or...) 
Anyway on to the next pint. 


W Yfou can't jump from outside a DO loop to inside a 
DO without going in viathe PO statement. 
Vfovlean't) have 


30 I=1,10 
, f § These loops 
al or cut across 


one anether 


Suppose You had 
5 4 


GOTO 80 
DO.80 LL=1,M 
f Sf 


UPhat’s all right, r\ ght ? 


Statement SO is inside the range 
of the DO, so the GOTO B0is 
an illegal jump into the Do! 
ven 3% it isn't always egal itis 
certainly immoral and proba bly 


fattening. Don’t do it./ 


¥Y [tis STRIGTLY FORBIDDEN to alter 


the imtial value , ina! valve, or incvement 


of a DO while you are in the loop. 
Yfou(Can't) have 


DO 4% I=1,M,ISKIP 


3 
T=1I+| = These are all 
ff f illegal ! 

M =MLAST 
Sf 
ISKTP= 2% 1SKIP. 


44 


tirhose three statements ave perfectly 
legal outside the loop, howevey! 


To Varying degrees, this 
treatise offends my delicate 


sensibilities for a hundred ax 
reasons, but it’s too late eae . 
to get my money back. w3 
WHAT CAN I DO? 


NOS 


rhe answer to the maiden’s prayer is she? 
can stuff her list of reasons into an ORDERED 
ARRAY. LA rrays provide a nice, systematic Waly 
to handle whole lists of numbers, Our friend has 


one hundred reasons for hating this magnum opus: 


R Number : ‘How Bad 
or Subscript : Reason : Scale Coal joe 
(\=Mildly Offensive 
se veee cceusttigese cesessateeterssesesan WO=Potrid) 
| Vulga r q 
3. TooMany Equations: | 
4 ‘Undefi itting a Distinguished: & 


‘(ol ve of Nobel Laureates 
5 Ghawiristicloward Giraffes 2 
*Note: The term She" is used as A Non-sexist personal 


pronoun. Lf You prefer, read ‘He for’ She’ and insert 
a picture fa boy aivaffe. 


Pro ramm in ; (Overflow error) 
99: InsultsMy Intell igence: 4 
100 ‘Didn't Gome With 7 


Grayons 


Dt wovldbe a real Pain & Nuisance to need to 
think vp separate variable names for each reason 
on the list. LA \so, it wovld be hard to mani plate 
such a list when you wanted to perform the 
same operation on all the elements of the list. 
Yfou wov\dn’t want to need to dez\ with each 
item individually. Por instance , You might want 
to dump the entive list onto a magnetic disk, 
$0 you could frisbee it inte a waste lasket/ 

Since the sadists who invented FORTRAN 
felt Yov woulda susfer enovah Ane when You 
Came to mani pul ating Aras of numbers, they 
Kindly made specia\ provisions Soy hand lin arrays. 
You Can give a Single ARRAY NAM to OW 
whole list & numbers and then vefer to the whole 
ball of wax by that one name! You can also do lots 
more, Me\chiov, but \et’s start simple! 

LAL the Very Beginnin oF a program before 


ee, 


the first executable FORIRAN statement , you must 


tell the com puter to set aside enovgh storage 
locations Cor any avvaus you lan to vse. 
his 15 done by means of a DIMENSION 
Statement, Since ovr friend with the twent foot 
larynx hated this literary master piece te }00 
reasons, she could ral the computer to set 
aside space for all those reasons by sayin 


DIMENSION TIT SBAD (100) 
right at the start of her program. INlow, 


She doesn’t need to make up a separate 
variable name for each separate veason, She 
Can refer te JAALL the YEASONS by the single 
Array Name , to wit, “ITSBAD,’ A particular 
*yaison disenchantement,” or Glement) of the 
Orray covla then be specified by means of 
its SUBSCRIPT. 

One covla SAY (Subscript 


ITSBAD (|) = 9 Mabe 
ITSBAD (2) = ¢ 
— ) = | 


IT SBAD(I00) = 7 


hen, if You wanted to ca\cv\ate just how 
bad this book 1S, on the average, you covld 
sim ely say: 


KRUMMY=0 
DOIO I=1, 100 
10 KRUMMY = KRUMMY+1T SBAD(T) 
ROTTEN = KRUMMY 
HOW BAD = ROTTEN 100. 
WRITE (5,20) HOWBAD 
20: ‘FORMAT (WE ALLKNEW IT WAS 
\:bBAD BUT THIS 15 RIDICULOUS’, 
?2:F 10.4) 


Rules for naming arrays ave the same as 
those for naming anu cther variables, Dn the 
example, we haa a One-Dimensiona| Integer 
Avray, since the array name started wi Lh 
an 1,9, KL, Ml, or N and since only me 
subscript was needed jin order te specitey any 
particy lar aval element. (Notice that we 
changea the Tota) Terrvibleness from an inte ger 
valve to a. Fea) valve before Figuring ovt the 
Avera ge Awful ness ) 


\\ ac ° Suppose I wanted to 
Oe | Give more precise valves 
5 to New much I dislike 
iL 


this tome ? 
NO 


[For more subtle , finer gradaxions, one 
could VS5E A Real Away of values by iving a 
Real name to the array. When you covla 
differentiate with nicety between degrees of 
a gvavation like 6.3758 and 6.376025. 

hat mi ght not seem important to you, but jt 
could mean al] the diffevence in the wovld tothe 
Daturday Review of Literate. 

LA\\so, one might want to have data on how 
much each r eC WAS Bisliked. You covla set UP qe 
table like the Following : 


PAGE NUMBER 
| 23 4 ee Oe 
92 44 4.1 83725 99 -.-.. 5.15 4.9 


4 [759 619.99 - + ee ee ee ee 
oe ee re 
QQ). 64 e os & She BoP Pc oe 
100 17-1 65 83 45.-.. 0.2... FBS 


Since there are already w A pages ana 
lOO veasons, one wld need to set aside 
100 X 96 stova ea locations for tne elements 
o£ tis TWO DIMENSIONAL ARRAY 
Not only do ye need to save these 100 locations 
but each one needs tobe big enough for a Real valve. 


Tlowevev, since you can't change the storage 
allocation for arrays while the program is runving, 
you need toe set aside at the stavt” as much 
storage Space as You think you're ever going 
bo need for the array. Since the sneaky Dy. W 
has already slipped another half yoo in on VS 


and since there is no sian ot the end in 4} ht, 


locati ons, 


Don't set aside too 
much extra space, 
ov there may nat 
be enovah Yoormn 
left jin the com puter 
to solve yer problem | 


Once again, space iS veserved by means Fa 
dimension statement at the very start or the 
program. One cov\Aa say! 

DIMENSION AWFULL (100,200) 
AWFULL (1,1) = 4.2 
AWFULL (1,2) = 9.4 
AWFULL(1,3)= i 


° 
AWrUut. (\,96)= 8.9 
AWFULL (2,1) = 6.3 
AWFULL (2,0) =2.6 


AWFULL (2,3)= 6.9 


f 
AWFULL (2, 90)= 5.7 
i 5 § 
AWFULL (l00,))= 7, | 


5 2 ; 
AWFULL (100, 96)= 9.8625 


[Further elaborati n4 en this AWFUT], v 
example y one covla refer te the element jn the 
12 ™Row and the [5% Column as AWFUIL (219 
ana one covla refer to the element in the mh 
row and the 2n-\ ™ column as 


AWFULL (M, 2% N-1D 


Rules for subscripts VA from computer te 
computer but on many systems Svbseripts Can 
@aly) be either an 

INTEGER CONSTANT 
INTEGER VARIABLE 
INTEGER CONSTANT * INTEGER VARIABLE 


or 


ONE of THESE TWO & INTEGER GONSTANT 
Subsevipts Glas) must be positive integers. You 


also can’t have a zero subscript. 
*What an awtol war to spell awful as a variable. 
Bind ovt wnat’s legal ov) your system ‘ 


Tellme, young man.Gould 
I use a subseript like 
65-MOM °? 


Why bless Fae little heart, it depends on 
n 


the com puter! small systems the most 


form for a subscript is usually an 
INTEGER NTEGER\ +¢ INTEGER 
CONSTANT) “ (VARIABLE) = \GONSTANT 


I subscrigt \ike 65-MOM weld get you 
thrown ovt on your corset,since it isn’t an 
integer variable fan integer constant but is an 
integer constant — a variable. When they say that's 
the most gen eval form , they mean it! On many 
other systems however, a valid subscript can be 
almost any sovt & an wibe ger expression, just so 
it is positive when You Gish evalvatin ‘t.@n 
big systems where the FORTRAN has ots of 
’s ana —. *5 65-MOM wold be finel 


While youre around ) Momma, show US hew 
You store AVYAYS in your drawers! 


Arrays ave stored 
by Gelomns. The 
irst column is 
Stored first, then 
the Second column 
and then the third. 
After that comes 
the fourth column, 
and so forth, 


Paraphrasing her 
elucidation, aYYaus are 
Stored in ascending 
Stora ge locations with 
the value of the first 
subsceri pt increasing 
most rapidlyand that 
of the second subscript 
changing less rapidly. 

Uppose you had an arr with 3 rows 
and 2 columns. fhe S\X arra elements would 
be referred to by = Toll owing subscripts 

Colum 


n 
2 


\Per haps this avrar is used ina program to test 
Fiede\baum’s famous theory in nuclear physics. 

( Piedelbaum principle on the Onset of Nuclear 
Destruction w Laboratory E xperi ments. ) arhen F 
using highly technical mnemonics relevant to 
the problem at hand ana Meaninatul te other 
phusicists when you publish You results, you 
woulda hand le this array as Follows on the Computer: 


a DIMENSION FONDLE ee 
he comepvter wovld store the elements the 
arya Gaternalld) 4) the Following order: a Dan 


st (Foe <n 


- 3h 
FONDLE (2,1) Sy a 


FON DLE (3,1) 
nd {rb (1,2) 


Go\lvumn 


ENDING 
ATION 


FONDLE (2,2) = 
FONDLE (3,2) |< 


Golumn oS 
at 


Wotice the computer stored the two-dimensional) 
array inthe form ofa one-dimensional] string of valves. 
Dt does that for You automati cally. You can always 
think oS the array as being two-dimensjona| with the 
\t cUbscript telling the row and the 24 subscvipt 
the co\umn of an array element. But keep in mina 
how thiinas are actually stored! Nlow 
that FOND i is stored, You can massage thedatal 


a Sa ee el 


FAscending location Numbers, that is! 


Mister... 
teach us te 
Read and Write 


Allo read ov write 
GA) the elements of 
an array in the order 
in which it is 
stored, just put 
the name of the 
array inan 
Input / Output list. 
[Remember, it’s the 
job of the associated 
FORMAT to 
determine how 
these valves wi\| 
actval|y look on the page, such exotica as 
how many to a line, what text goes 
between values, where page skips will 
occur, or whether the total effect will be 
O starkly cubist block of numbers ora 
rococo and delicate swiv\ of fiaquves in 


a style reminiscent of the great masters of the 
Igth century, 


75, 
3.7 3.2, 1 I7d 16.4 ye, 
4. a 4, . “L- Portion & matrix 
&O : 
15, 3.1 2 7 printed in thestyle 
5. . aces 
ie Bi 7.14 g,7. 4. of the great 
30, adil Givse ppe Galli Bibiena 
oq ««3- 4G y) 
q 31.9 $35 Cie \'7] 
23, 4.17 74,62 a 
. 28.40 35.19 : er 
. TPES 44.15 Metropol 
4§ 52,14 a Gourtesy ,Metrop 


Lyseum cf Art 


Cphat style of computer printovt reached its 
peak in 1752 with the renowned Flemish 
programmer, Lieder Krantz and the printout of 
the Frenchman , Fromage When Johannes 
Guten berg invented the line printer with 
co\umnay type, however, that style began bo 
decline, thovah it is still practiced today by 
man beginning programmers. 

lence ana forsoo th , one wovld \ikKe to know 
just what a FORMAT will do when given a 
string of numbers to write ovt, Knowledge 
1s power anda well-wr}j tten FORMAT Can bring 
a \|200 line /minvte printer Lo its Knees! 

First , io is comforting to Know that if there 
are more values to be printed in an 1/0 list 
than the Format allows for, the extra numbers 
will not be brutally squashed against 
the closing right parenthesis | 


Suppose you had 
DIMENSION PIMPLE(20) 


WRITE (5,10) PIMPLE 
0 FORMAT (1X, ie 


i 


Thlere ) the format allows for Five vea\ Numbers 
to aline, witha starting blank te get past 


the carriage contro|. So the first line of 
output wovld have the (Values of 


Bmele(i) Pimple2) Pimple(s) Pimple (4) Fimple és) 


Bvt wait | “rhe write statement is trying to 
Send Pimple (6) te the Format, but we're 
smack up against the clesing right parenthesis. 
“Phat would never do, because you should 
never squeeze a pimple, so the computer 
automat cally inserts a slash and goes back 
and reuses the Kormat again! Whenever 
there are more valves to be read or witten 
than OW Format allows Tor, the com putev 
mserts a slash and starts over. thus, that 
FORMAT#lo behaves as if it had been written 
lO. FORMAT (1X, 5 F1I04/1X, 5F10.4/1X, 
ESF 10.4/\X,5F10.4) 


You just need to be sure that if the 1/0 \ist 
1S trying to stuff ovt an Integer , the format 
swt waiting for a Real at that point, or, 
OSs they say here in Washi ngten , ak that point 
in time. 


aan er What, ah, what 
if theve ave parentheses 


inside the format? 


Lm glad you asked. 
That’s an excellent question 


fora PhD candidate ! 


NE there are parentheses (inside) of a 
POKMAT, the anselal doesn’t Gree way 
back to the start of the Format. Instead, 
ue qoes back te the start of the most recent 
inside set of paren theses, For instance, if you 


7 
PAA 20 FORMAT (IX, 215 /(1X,3(F104,52)) 
and you ran ovt of Format 


and still had more inthe last Jside 
1/0 \ist to print ovt the Farenthesis 


Computer wovla inserta slash 

and start ag ain from this point; since itis the 
left parenthesis corresponding to the last right 
parenthesis inside the format, 


+ Wov saia 


DIMENSION VELHOR (49) 

WRITE (5,20) LAMBDA, XNEW, VELHOR 
20: ‘FORMAT (1X 15, 5%, F104 / 

1X, 3 (Feo 2X))) 


and the current valves were 
LAMBDA ....163 VELHOR (5)...72.54 
YNEW.....---25.381659  VELHOR(O..-17.32 
VELUOR (1)....172.36912  VELHOR()...\000 
VELHOR (2). ....19.19958 VHLHOR(s)...-15.19 
VELWOR (3). -.7285-500° =WVELHOR (9)...18.9 
VELWOR. (4).---. 17.0 


then the compoter wovld spew out 


bbbi63 bb bbb bbb25.38 16 
bbbi72. 3 G.bbbbbI4. 17,bb b- 285.50 


bbbbl7. 00.bb bbb72. 59 bbb b- 17.32 
b bb 100 . 00 bb bb- 15. 14 bb bbbIs. 90 


A 


where the smal] “a” marks have been used to 
show the breaks between numbers. OF covrse, 
the b''s wovld be blank spaces on the actual 
printout, but they have been put in here to 
show Now the valves line up. 


IN|} otice how the nine elements of the VELWOR 


array have been written ovt on thvee lines as 
belhor, Ve thor, Velhor, 
Velhors Ve lhor, Velhor, 
Velhor- Velhor g Velhor, 
[Bach element gets printed jn the next available 
Field. 


| dan.... S up pose i 
didwt want the 
whole array y but 
\ust the 4% element? 


Nt Yov want to Yead ov print just On 
parti cv\ar element of an AVY AY y sim ply put. 
that element into your Input /Outpot list, 
\ust like any ether variable. Say 


READ (8,40) X,LL(,'22(3,1) 
WRITE (5, 65) VELHOR (4), XNEW 


When an array name appears in an 1/0 \ist, all 
the elements of the array are Yead in ov 
written ovt In the order they Appear in storage! 
When an array element js inan 1/0 list, just 
that ene element gets transmitted. 


Mo ving on briskly to another example for 
your edifi cation, consider: 


DIMENSION KVETGH (2,3), M (7) 
READ(8,10)M (4), KVETGH,X,Y, Z 
10: FORMAT (15/3 (110, |0X, 1107) 
13 F104) 


Suppose the data cards in the card reader at the 
time this read statement was encountered were ! 


Flow do You think the data valves woula be 
matched up with the variable names ? 


Tl] tell you, 
Ma gets 236 
from co\umns 

I-5 of the First 
data cara. 


Arhis hap ens first because M (4) is the first 
list elem ae first field speci fi cation in 
the Format is Is, which means that the first 
5 columns on the data card. Wi Y contain the 
value of Ma, right justified with the last 
Aigit of the valve (6 in this case ) falling in 
column 5, 

Arhen , the whole KVETGH array gets 
stored, but the way you might quess at 
First | eres how it aets stoved: 
KVETCH,,...I5 KVETCH,,...1B KVETCH,....45 

This was the 2" This a valve came } This was the 4 
data valve of f the 34 card / 
_1O-KVETCH,»...6 KVETCH,,...12 


2,\--* 
Here's the 34S and the 5% 7 ana this was 7%! 


Notice that the KVETCGH elements get read 
in in Golumn order, just as they appear internal 
in the computer's storage. To hus KVETCH (2,I) 
Comes right after K ETGH (11). [But the 

format speciticati ons fov the KVETCH array 


/3(110,10X, 1107) 


which means the |*KVETCH valve Comes 
from columns I-10 of the 2°4 data card 
(because oS the / bet ore the S. Yov go to the 
en cara )5 the er aKevETCH valve isin co\vms 
21-30 on that same card; the slash Says goon 
to the 34 card; and now, because of the 3 in 
front of the parentheses jn the format you 
do the same thing again ana again. (Phe 3 
is & repeat constant y Wlotice as Yyov 90 aown 
acolumn in the actva\ matyi x You ave Yeading 
across a data card! [Rows and columns get 
transposed or interchanged on input and outpt 
because valves get sent to ov fromthe format 
in the ordevy in which the y ave stored internally 
wn the com puters memory. 

[Finall ) Xa Ws ana Z straggling along at 
the end the 1/0 \ist get read in as 

K?...23.68 Y..1%2192 Z...-3.8¢ 

Bince this was Input, the Aecima\ points on the 
cards didnt need to )ine up with the FORMAT! 


But Suppose 1 don't want 


my array transposed ? 


Yo ung Man . You'd better 
stick to Politica| Caricatures. Sy 


Don’t you know The Womens y 5) f 
Movement Has No Go 


Sense of Humor? VF 


No. Hum a few bars Sit. a 


To read ov write array 
elements in an ovder 
other than that in 

which they are stored 
internally, or toveaad 
ov write just part of 
an aryay... 


You Can vse an 
IMPLIED DO 
LOOP! 


Note: 


Dr. K’s wife feels some people won't 
appreciate the satire on the preceding 


page. UF you find that page offensive, 
please disregard it and use the 
following page IN instead.! 


tote: IF you find the possessive term” Dr. K’s 
wife “offensive, disregard preceding note. 


Youna man, 
Do you Know 
This Book Has 


No Redeeming 
Social Virtues? 


No. Huma few bars of it. 


Frequently, situations arise where you 
don’t want your array trans posed ov where 
ov want to write ovtanam in an order 


other than that in which it is Stored. 


To read or write ary 
elements in an order 
other than that in 
which they are storea 
internally or to read 
or write just part ofan 
array... 


You can use an 


Kn IMPLIED DO is sort of like a 
yeavlay DO that was levelled by o. steanveller. 
oF @ wiole loop has been laid ovt Sideways 
inside oF a Single FORTRAN read or write 


statement. 


T’ve been laid out 
sideways For years and 
nobody made a fuss! 


Lhe form of an IMPLIED DO is’ 
(List, Spt = Ms ért, Mend, Incr) 


Lhi tial 
Valve 
Do Govnter 


(Integer Variable) Integ er Variables 


or Gonstants 


This whole works /s a lisd item in an 

input /ovt pot list / 

[-\ s You have indvbitably quessed , the 
parentheses in the Implied Do mark the 
range of the DO. the way this all works is 
very simple, It’s almost exactly like a reqular DO, 


except that the entire thing , range and all, 

1s contained in one statement. Here ave sone 

examples of implied DO loops in T/0 statements: 
READ (8,10) (SGORE(K), K=1,100) 


- ee — ((KVETCH (1,5),3=1,3), 

Ht: T=l,2) | 

i WRITE (5,50) X, ALPHA,M, (SCORE ( 

WIGNT), ICNT=1,M,3), XLAST 

What happens is the DO counter is first 
set to the value & Mstrt. Phen the list 
items in the range of the DO are shipped 
off to the associated format, The counter 
is then incremented by Lner, tested , ana 
again the list is shipped off, but using the 
new value for the DO counter as a subscript, 
Bteetera etcetera etc., ete. 


C Where will it all ena Miss Jones? 
(Asthe doctor asked the nurse in 


the maternity ward. 


Dt all ends when the incremented counter 
ov\a exceea_ Mend! 


Verily, Yyov may ask," Who Cares?” 


1 WANT YOU 


TO CARE 


By being self - contained within a single I/O 
statement, theImpliedA DO allows yeu to contro\ 
the order in which subscripted array elements 
are shipped of to a format for printing, or 
pulled in from cards on input. [Remember, 
whenever you encounter a Read or Write 
statement, you start at the Beginnin oF 
its associatea format. this is trve if 
you have never met the vead/ write before, 
say when bumping into it for the First 
time axt the start of a program. Dtis trve 
if you meet in some secluded rendezvous, 
say at the end of a seldom used branch. 


[nh] owever, it is AUSe as trve even if yov have 
been meeting veqviar\ » SAY in Some 
delight ful litle loop. You can have been 
oing around toaether for some time 
Mone - the - less, Each) time you Yun into the 
read ov write statement, You start 
over at the beginning of its associated 
Formak, 

\ So what ,” you may be tempted to 
ask. LA w Go ahead, ask ! 

is means you cant use a re vilar 


DO loop to didd\e arovnd with oa 


subscripts For stuff to Appear On aA 


line of output, 
SOs AD (ypu ask) 


[Because whenever the DO Changed the 
subscripts of the valves to be printed, jt 
wovld be (outside) the write statement, 
Phoo & Bah Bach time throug h the 
loop, when you hit the write statement 
hoping te print the next valve on the 
line, Yyoo'dA find yoursel+ thrown back te 
the start of the Format again! 


FORMAIR Gg 


WS 


@vppese youhaa 
DO 20 T=1,3 
Ke ool 
WRITE (5,10) K 
lO: : FORMAT (\X,12,"b MOCKING) 
EBIRDS'/\X,12,'b FRENCH bHENS’ 
2/"bANDb', 12,'bPARTRIDGE, *) 
20 GONTINUVE 
Unfortunately, that wovld write 
3 MOCKING BIRDS 


2 MOCKING BIRDS 
| MOCKING BIRDS 


making a mockery of yovr efforts 
Equally futile wovla be 
DIMENSION ere (2,3) 


DO 50 1=1,2 
WO OC 37155 
50 WRITE (5,70) KVETGH (1,3) 
70 FORMAT (1X, 3F10.4) 
hirhese are lain, ordinary gar Aen vaviety 
nested do Veops. Athen sequence KVET GH bbyrows 
instead of columns, but the ele ments unSoytunately 
get printed as one long string of numbers! 


So give them a 
‘For instance . 
of what an 

Implied De can do | 


For instance , SUPpese You had 


DIMENSION BANANA (320) 
READ (8, 10) (BANANA (K),K=1,3204) 
10 roRMAy a 


tthis one vead statement wovld read in 
80 data valves for BANANA (i) BANANA(5), 
BANANA ( Ds and soforth. BANANA 5 would 
be bunched eight toa Cara , so the one yead 
statement wovld read in 10 cards, selectively 
storing successive numbers in every fourth 


location in the BANANA arvau, Ave te the 


increment, 


U$ the increment in an IMPLIED Do j, 
to be a 1, then Yer Can leave it off , just 
as jin the Case of a ve vlar DO. LA\\so, Ov 
can write ovt the valve of a DO countey just 
like that of any other variable. Wake a look 
at the following example: 


WRITE (5, 10) (M,X(M), M=1,20) 
10 FORMAT (2 ('bbX (', 12, ")=",F]04)) 


would write out the first twenty elements of the 
array named X” in the following fashion: 


Aa es), Xa 
X(3)=esee5 x4@)= o> 


[Here , each element was labeled b writin 
ovt its subscript along with Some text , 
just before writing the element itself. Arhis 
0C5 a long way towards makin the printout 
Understandable at a glance, without needing 
to vefev constantly back to the program. 
F\n IMPLIED DO can also have 
other Implied DOs nested within its I/0 list. 
For example, remember ovy ola friend the 
KVET CH array You have been trying so hard 
to forget? Ar here are severa| Ways You Can USe 
implied DO’s to write it ovt without having it 
trans posed. One wat is by Using an implied 
DO to move you across each \ine and a 
requiar DO to move Yov from Yow bo YOw, Af 
YoU are a lutton for punishment, you can do it 


al\'m one WRITE using Nested Implied DO loops ! 


OW Ret gc us bey 1c ae, 
i KVETGH (1, J), 3=1 ,3),1=I 2). XY, 7, 


20} ‘FORMAT (1%,3F 10.4 /2 (5X,"BOWl 
i yL1,'>b’, 3%, 'KVETCH ( y3N,’," yil, 
2:")=" T10, 3X,"KVETGH(",11,","5 ll, 
3=", TO, 3X, 'KVETCH( , IRR TL, 
4")=", T10/) "bX=", F10.4, as SY", 
5IFI0.4, SX, "bZ=’ F104 //) 


esd Gries / 
Be 


orhis Sescenine and formidable Format 
freak’s farctate Foaming would produce 
fecund ana Feracious phonics on the 
line printer, 


W har’s he talking abovt? 


El A Bot Tao Say | Bhat lycidly 
written FORMAT wovla produce the 
Tollowing printea PA9e, taste fv || y annotated’ 


ROW() KVETCH (1,)= 3 KVETCH (U) = KVETCAR: 
ROW(D KVETCH (2,))= eee, KVBTCB (20=) RETA 3 


where ) of Course the S4ul qgles ms wovld have 
whatever numbevs were stoved in those variables at 
the time of printing. 
That’s the wo leoking printout 
/ 


) 
I've ever seen. 


It was done with 
a chain printer. 


T hear the qvy wi th the 


chain broke both the 
printer’s wrists! 


Momma, youve done a 


LOVELY 


job arran ging your drawers. 


‘You don't mind if I suggest 
you put some variables 

beginning with I,J, K... 
in big drawers pst to make 
your chayming house a 


bit less dull... 


Integers fit in sma\\l arawers and 
Floating point numbers only fitin big drawers, 
hrivring B. Octal, the famous Computer 
science theoretician, said it best in his 
of L-quotea “Principle of Gonservation o$ 
Drawer Svace” in which he said 
\Entropiano al numero integerto coda alla 
Vaccinara minestrone Ai Fai ol, itsa. nota 
gonna fit alla teenetzia drawer!” 


which rovahly translates as “You can’t cram 
a bia number into the same drawey with the 
socks,“ inlowe vey, it (is) possible to ive a big 
dvawer a name that world normally be 
associated with a small drawer, 
you tell the computer that You are 
deliberately switching conventions for that 
particvlay variable only. 


AS long as 


dah?? 
NU 
Un other words, Momma set UP Asa 
purely arbitrary convention that variable 
names starting with IJ, ,L,M, or IN 
would be understood to be integer names. 
Following that IMPLIED PE convention 
where the variable Lupe is implied by the 
First letter in its name, helps momma Keep 
her drawers neat. She’s less likely to make 
a mistake like trying to put a. Yea\ number 
in a sma\l drawer or a girdle in with 
her dentures, 
TH] owever, sometimes for mnemonic YEASNS , 
oY even platonic veasons, You might want te 
override the implied type and vse a variable 


name like GOUNT as an Integer oy LAMDA 
as a Real. 


Yfov can do this,so long as YoU tell the 
computer explicitly that that is what you 
plan to do, a, the implied 
convention, you start your program, with 
on EXPLICIT TYPE STATEMENT 
in which You hist those variables which don't 


follow the implied convention. 
For instance You covla Sau 


REAL LAMDA,K4D0G,MAX , LAST 


‘| INTEGER COUNT, DIME (10), DOLAR(I20), 

‘I! PENY, ZETA, Z1G,AND, ZAG 
Erhen. throvahout the program. the computer 
socth, kao nal eae the usual me 
Conventions, DIME was an Integer arva 4 
name with ten elements, If an array name 
appears in an Hxplicit T pe Statement, then 
all the elements of the avra get that ein 
Yfov can either a\ve the seh the 
array right thenanda there w the type 
Statement , or else you can ee the 
Aimensions in a following D TEN SION 


statement. (Gotcha!) 
Me 


Uncidentally , just be cause Dy. K Vses 
droll yvaviable names like K4 DOG doesn't 
mMéan You shovld. Dr. K can’t stand it in 
other people's programs, Hle’\l senda note 
home to your mother! It’s a terr} ble 
disg usting habit to aet into. Youv program 
wil} get filled with fleas while You try to 
vemembey what K9DOG stands for. Stick 
to nice dv\l meaningful names like PI and 
XxX MIN. four giv friend Or boysriend's NAMNC 
may be meaningful, but leave it ovt of Yovr 
pro qram too , especiall if Yov ave mavV ied . 
When Drv.K writes actual Computer programs 
they ave duiler than the Vladivostok 
tele phone Airecto v4) but they wovk! 
Gomputer time costs at SO vy to make 
aeuNe [eeqrans as clear and legible AS 4ov 
CAN $0 You wont blow a \ot of runs | be thrown 
ovt os school ger $s ver) {ror WYouv job, get 
bitten by your dog, have your car repossessed, 
Your mort gage cance\ ea, You” wife leave 4ou for 
oO trave| ing calli ope salesman andyin afit of 
Aespair , throw yourself under the print hammers 
On a page printer and be pummeled to death. 
That's what happens to people who 


write cute variable names! ! 


124. 


That sounds 
sort of Logical | 


oe 


Speaking of Legical, there ave cther Kinds 
of variables besides those that stand for 
integers and veals. In particular, You can have 
Logical Variables which would be drawers 
that would either contain the Constant valve 
STRUE. ov else the constant FALSE... 

tre tell the com puter that certain 
variables are lo ical ana will only be set to 
TRUE. or FALSE. , vse an explicit type 
statement like 

LOGICAL BOOL, DEBUG 


Athen, In Yovy program, You can assign 
logical valves te these Variables by 
statements such as 

DEBUG =.VRUE. 
or BooL =X.GT.Y 
In these assignment statements the expression 


on the Wight of the _. sign is a 
Logical Expression. Bither X js bigger than 


Y or else it isn’t. UF it is. then BOOL is 
set to .TRUE. Otherwise, BOOL is .FALSE,. 


Hlardlu a big dea), but sort 
of a middle-sized deal, Later 
on inthe pro ram, You covla test the 
contents of BOOL oy DEBUG and use them 
to help decide what to do next. For instance, 
you could say: 


IF (DEBUG) WRITE (5,IOXY,Z 


and if DEBUG was .TRUE., then X,Y, and Z 
would get written ovt according to format#]0, 
Otherwise, that particular WRITE wovld be 
skipped. While you were devel opin Yovr 
program, you covlad have a let these 
write statements print ovt intermediate 
results, Once the program worked , by 
Changin one cava and setting DEBUG & 
FALSE, , you could suppress all the exta 
printout. DE Yov ever wanted it Again theugh, You 
covla get it back. Don't letthe A. 1.A pet DEBUG 


in your pro gram: 


This sounds 
awfu) GOMPLEX ! 


Dt really isn’t, it’s 
simple, but speakin 
G omplex you Can also tell the computer to 
use complex arithmetic in its ca\cvlations 
if Yov want to. If along with the other type 
state ments at the start of Yyovr pregram, 
You throw in one that Sans something like 

COMPLEX VEG,Z2,DELTA,WVEC 
then that tells the Cornputev that VEG and 
Ze are not plain ola variables from Plains, 
GA, but are complex. Jf DELTA is to be 
OO Complex quantity ) then put that variable 
name inte the GOMPLEX type statement Athen, 
the computer is ready when You are. 5 OW) 
| vp never heard of a complex Number 
haloes, you’ A better skip right past the 
next few P2465 » ?cavse is apt to be 
con fusing though amusing at a.casval perusing. 


[Because complex numbers have two parts, 
they make a dandy way to mani pulate certain 
Kinds of problems in which you have two 
humbersto Keep track of. For instance, it’s 
hard te stuff a road map intog computer, 
since the machine seldom drives into an 


Exxon station. 

Fill her up,ana 
check the cara 
bin! 


[Nlone- the-\ess, it’s 
easy for you to handle geometry 
on a Computer if WOv Use Com plex av ithmetic. 
One part of Yovv complex number can represent 
Yovr xX coordinate and the othey part, your 
Y coordinate. Woila | What could be easier! 


Instead of needinato Keep track of yovr 
Jatitvde and longi We separately, you can 
represent Your positi on \o ene com plex valve, 
‘There Aare \ots of other Good things Yov CAY) 
Ao with complex numbers , SO if Ov don’t 
Know much abovt them , L commend them ts you 
for further serutiny, 

Since a complex number consists of an 
ordered Pair of rea | Numbers | the Com puter 


sets aside two storaae locations for each 
complex number, [Nlow, a5 everyone knows, a 
complex number consists of a real part andan 
Imaginary part Where the Wmaginar part 1S 
as réa\ as the real part if Yov Can imagine 
that. What may seem unrea) is that the 
real meaning of the wovd ‘veal" doesn't 
mean real.” Fukesd. 4 means “real.” het me try 
that once again. ‘Toa mathematician, the word 
“real” isused to talk about all the rational 
ond iwration a\ numbers, IN ow we a\\ know that 
FORTRAN is ivrationa\ 5 ana Bolzano’s 
bambino would have wet his diaper 1f he 
Knew that computerniks A > 
would Someday vse the 9) 
word “veal” to mean be 
“numbers with decimal 
points,” e., F\ oatin , 
point numbers. Aarhus, the wova ‘veal isa rea} 
bummey becavse it has so manu meanings it 
can make YoU veel. Usually what is meant 
is é- perfect clear Srom the context, however, 

RF\ COMPLEX CONSTAN T ic written 
in FORTRAN as 


(3. 625, - 25.4) 


oFhat 1S, itis written as a pair of real WB, 
(floating point ) constants Separated by a 
Comma and suvvounded bu a set of 
parentheses, What exam a shows the comple 
number 3.625 ~ 5$25,4 OY, if Yov prefer, it 
shows the complex number 3,625-125,4. {t's 
the Same number — which notati on Yov vse 
just indicates whether Yov AYeC A mechanical 
engineer, electrical engineer, mathemati clan, or 
an itinerant Peanut farmer, Wn an case, the 


lL ov | simply represents V7". hus, in FORTRAN 
(19.16, 4.85) 


means “the complex num ber whose vea\| ayt 
is 19.16 and whose imaginavy part is 4.35,’ 
+ You wrote an assianment state ment 
Usin complex Variables and constants, it 
wovld \ook pretty much \iKe any other 
assignment statement! 
VEG=((0.,1.)¥Z2)/WVEC4DELTA 
Arhis means,‘ multiply the complex constant 
(0.,1.) by the complex variable Z2 and then 
divide that by the complex variable WVEG, 
[Final jada on the complex variable 
DELTA and then store the result in the 
complex vaviable VEG. Use com plex 
avithmetic Lov all the calevlations, since 


the numbers are complex.” 
0 read or write complex variables, 
simply allow two real fields for each complex 


variable in the I/O list. For instance, 
WRITE (5,10) VEG 


10: ;FORMAT ("b THE REAL FART OF, 
:VEG 5b", F10.4,",AND THE IMAG 
:@: INARY PART 156’, F10.4) 


Generally, you shovidn’t mix comPlex 
numbers with reals or integers unless You 
Know it is O.K. on UY ovr com puter. IF ov are 
still reading at this point y You must ve 
interested in Computin y 20 it is worth our 
while to qe out and buy ao FORTRAN 
lang vage manua | published by the Fine Folks 
who designed your particular computev.Afheir 
manval will tell just how thinas ave supposed 
to be done on theiv system. Who knows. It 
might even work like they Say, 


WishSul thinking ! 
We have a KLUDG 
7000 


oro he\p you use complex arithmetic in 
Oo program jalon with conventional avithmelic, 
theve are some specia| hlorary Concti ons 
that ave nice to know about. You aan lok 
vp all the available functions in the manua| 
yov just bought duving the last gen 


but here ave some of the most vseful ones: 
CMPLX (X,Y) 


Sa handy danay little Cunckien that 

takes twe Floating point variables and 
makes them inte the real and ima inary 
parts of a tom plex num bey. You cou\ ad 


VEG= DELTA + GMPLX (PIBETA) 


ro find out the real” part of a complex vaviable 
You could Sa 


REAL (VEC) 


Hrhis shovla not be confused with the REAL 
explicit type statement at the Start oF Yovy 
program! “Ihe REAL (VBC) is atotally 
aifferent thing. This isa floating point 
function that takes a complex avaument 
and Gives back the ‘veal? pvt as the valve of 
the fonction. 

AIMAG (VEG) 
1S similar, but gives back the “imaginary! part ; 


Not by a long shat! 


There's \ots of 
impor tant FORTRAN 
yet to come! 


i+ You repent , you've yourse\+ to blame, 
Yov shovla never have let me begin: 

[Bot speaking of the end , that is enovg h 
abovt complex nombers. Nlow You know they 
exist in FORTRAN ‘ LF You want to know move 
you can read about them in your computer 
manva\, 

IL et me te\| Ov just how near the end 1s. 
rhe very last cava in any FORTRAN program 


must be END 


Strangely enovah that’s called an END 
statement, End of discussion. 


Tell us ab ovt 


SUBPROGRAM 9! 
seh 
{ V7, A\\ my workas a 
cy Gleep requires the 
GRISELDA (ZIZZER PHO) 
function. 


UF, a Yovr professional CAYECY ADA Gleep 
You neea a \ess prosaic Cunction than the sines 
ond Square roots vsed 4 the viff - rat, then 
it iS aad to write Yovr oun special femetjon. 
Arhe function will bea complete is ea by 
itsels, usea alona with Uovv main progyvayn. 
Since it is not a standara program already in 
the computer, you'll need to add it to the end 
of Uovr Main Whee Wass but ahead of your 
data caras. You pre bab|ly need some special 
contro\ cavas in Yoo deck aS we || ) but they 
Vay Grom Computer to con putev So T can't 
tell Wov abovt them. 


(Rleve’s how You a write the Griselda fonction’ 


WEIRD FUNCTION GRISLD(ZIZR,PH) 
DUMB= SQRT(ZIZR)/(2.* PH) 
STUPID= PH-ZIZR 
IF (DUMB.LT.ABS(STUPID)) GOTO 10 
GRISLD = STUPID 
RETURN 
lO WRITE (5,20) DUMB 
20 FORMAT ('b DUMBb=b’, F10.4/) 
GRISLD= DUMB 
RETURN 
END 


IN otice that this fonction iS a complete 


Program by ibselt (Look! It'sa compete 
Ut ends with the (program by itself! 
dyaditional FORTRAN ~——— 


stake ment , “END " “his is required . Yfov 
would only be wasting time trying to ena with 


ARIVER DIRTY 


ov 
AU REPERVOIR 


“he computer only accepts ‘EN D", 
Wlotice also that you can tell from the first 
card that this isn’t an ordinary program , but 


is a weird function. When the computer SEES 


that Cox i, it knows just what it 1s dealing 
with , right off the bat. 

Sto use this Function, you would add 
ie bo Your cavd deck just aftev the END 
card of Your main \ine pro ram, Any hata 
cards wovld come later, along with any 
contro \ caras You miant Yequive for You 
particular com puter, hen, in ovr main line 
ged Nae, Yov woula Vse this just like any 
other vunction, Your maiy proavam might Say, 


SILLY= $5 R (GRISLD (X 2.x FOOL)) 
PQLY A R= GRISLD (SILLY, DILLY) 


Notice thet in this statement, SILLY 
ana DILLY would neeAto be REAL numbers 
With valves that weve defined earlier in 
the program. These Walves are the Arguments 
o-the GRISLD fonction. In the phunction 
icselt , G1ZK and PY are merely Dummy 
Ayguments or place holders.4Phe Valve of 
SILLY is used whenever ZIZR Appears in 
the function andthe valve of DILLY is 
used wherever PH appears. Whe variable 


names ave matched vp between the actval- 

variables in the program that is vsing the 
function and the Aum mu arguments in the 
sub program. 


PGLYAR=GRISLD( SILLY, DILLY) 


Svb Proar 


WEIRD FUNCTION GRISLD(ZIZR, PD 


Usina the values of SILTY and DILIX in 
place Pe ZIZR and PH, the steps cf the sub - 


ee a ave then cary ied out ana A valve iS 


assigned to the Lunetion ) etther by the 


statement 
GRISLD = STUPID 
or by the statement 
GRISLD = DUMB 
LAt this point , the subpreavam 5 logically 
(ov in this case, Megically) ovev,, so Control Is 


returned te the main pregram which asked 
Co tne Function to be inveked.Athis is done 


by the RETURN 
statement, 


dah... wheve’s 


the answer ? 


What was the questi on? 
Gh yes the ANSWER comes back as the 
valve assigned to the Name of the function, 
Nt works just like the cther Sunctions you've 
seen — ones like SIN & SQRT. Arhe whole 


works , 
GRISLD( SILLY, DILLY) 
ets replaced witha ‘Sjnal aum bey which 
is the pone of i eS 
IN] otice that this function happened to 

afer Crom One like SQRT, in that it had 
two ot pia Ponct; ons Can have ads Manly 
Aaravments a5 You like , but they can only velm 
a single final result. 


Whats all the vacket 


in the computer room? 


Its just two fonctions 
having an avqument ! 


ae. in ae twat when & main Procram 
a Sunction sub vain. the list of 
arguments ok the facies "is used asaway 
co get valves fom the main PrYoarawm iMtte 
the function. Urhe names of the arguments 
in the main progvam Can be entirely Aifferent 
+rom the Names oF the covvespendina dvommy 
arguments in the eub program. The we lists 
of avauments covld vse the same Names it 
makes no. difference to the Computer, ENN the 
compute CAVES 15 that the actva\ arguinentts 
agree with the dummy avquments in 


% Number £324 am 


x Type GOoABE 
* Dimension BIS 


In other words if the dum ag ava vments Consist 
o+ two veals an integer ayvray with ten elements, 
ank another real ) then the actual arguments 
must correspond, 

Because the MAIN yegrars 19 a completely 
sepavate eae fyom the Sub program ) Aou Can 
have duplication of statement numbers between 
the two progvams and not need to worry that 


the tom puter wi \l ger into trouble, For 
instance , the main cola have a FORMAT 
#10, but when the Sub progvam hit “GO Ta jo’ 
ie wovlan + bry to jure into the Wain. 
Statement num bers only need to be Un| que 
within the parti culay pre gram in whi ch they 
oCCUY. : 

Since, yn this case the function name 
GRISLD isa REAL Name the Value of the 
function wovld normally be returned asa 
REAL valve. Revles Cor namin functions ave. 
the same as those for naming any variable 
lad we Wanted GRISLD +, be ix ee 
function, (despite) its nam2, we. coulda have 
T'yped it explicitly by Saying 

INTEGER, FUNCTION GRISLD(Z\2R, PA) 

Type 
DS the function TYPE 1s implied by the name , 
then it can be omitted. Puncbions can be REAL, 
INTEGER COMPLEX, o LOGICAL. WEIRD 
5 aon appropriate , though iMegal + pe 
speci si cation. LA part from that sthovgh , this 
program wovld yun ' althovah it dves some 
weird things, 


My workas a 


Pure Mathem ati Cian 
ana an 


Appliea Bambino 
1s highly technical! 


é ae 


Lneed Severa\ answers 


back from muy subprograms! 


Well Now) we cant let Bolzano’s 
Bambing burst his bonnet , can wes Dn 
that case, write a SUBROUTINE 
Subprogram instead of a function. AA svbrovtine 
returns valves to ‘ts calling preg 

Yoo Hee! by storing New valves 

Subrovtine! into some of the names 

acl in its argument list. 
ithe name of the svbrovtine is 


vsedas a way of getting back the answers! 


Subroutines are of the form: 


: ‘SUBROUTINE NUTS (WALNUT, 
i: GASHEW) 
DIMENSION WALNUT(5) 
DO1O NUTY=\,5 
lO :: WALNUT (NUT Y)=CASHEW ¥ 
‘]; WALNUT (NUTY) 
RETURN 
END 


Arhis subroutine subprogram miaht be used 
by a main\ine program of the follows ng ilK: 


DIMENSION PEGAN(), HAZEL (3) 
READ(8,10) PEGAN, HAZEL 
10 FORMAT (5F104 /3F10.4) 
WRITE(5,20) PECAN, HAZEL 
20 : :FORMAT ('bPEGANS ARED', 5F104 / 
TE"bHAZELS AREW, 3F10.4) 
CALL NUTS (PECAN, HAZEL (2)) 
WRI1TE(5,30) PECAN HAZEL 
30: :FORMAT (">NEW PECANS AREU’, 
NiSF10.4/ ">HAZELS STILL AREL’ 
21, 3F 10.4 ) 
CALL EXIT 
END 


Note that the kK, are 


7 
invoked b Y CALL Statements. Where are twe 
svbrovtine subprograms being vsed b this 
particvlay mainline proqvam. First Jet's talk 
abovt the Second.trhe second we're throvah 
we'\l discuss the first. 

EXIT isa program which the Cow pute v 
\ockeys who Aid the S$ stem peg? Cov 
Yovr com outer probably Wwéve goed enough to 
provide, It sits in yor torn puter's ey manent bag 
of goonies y alona with SQR Tl and SIN, Whenever 
Program Yov write is all throvah doing its 
thing and ready to call it qvits , don't. call quits, 
GALL BRIT .~CKIT then escorts Your Program 
off the computer, calls the cleaning lady to bid 
vp Mommas Avawers , and ushers the next 
progvam in. On some systems Yov say 

OY 


insteca of GALL EXIT. What's the difference 
between GALL BAIT and END, yov inqvive © 
WN tel you. END 15 always the Very last Card 
" evev\ program. What it Aoes 1S te the 
Computer \this 16 the VeryLast Cava in this 
program. GALL BAIT yon the other hand \ can 


o anuwheve jn the ero Yay thet ov Want 
Jt. You can have ceva | GA Nols EXITS. 
one at the end of each branch of your 
flowchart that comes bb a dead end. 
nevey youv program decides it isall 
Yhrovah , it tells the Com puter- ej thev by 
saying “STOP” or else “CALL EXIT", 


Wow , \et’s get back tothe sample main 
program and the GALL NUTS (PEGAN, HAZEL(z)). 
Drhis ca\l statement invokes ovr homemade 
Syb program which we will add te oOvY deck 
alon WA th an other functions OV Subrovtines 


iaht want, ' 
“es oo We dent want any! 


a 


the ar uments ina Subrovtine Sub pregram 
work like a twoe- es street .AMhey are vsed 
(Goth) to get Valves into the subroutine ana 
also to yet valves back ovt! his is 
different froma function ASheve the arguments 
only lead a you Gan'b) get things out threvgh 


the argument list of a function, 

Observe , sil vous plaik,y that theve is a 
one-to-one corres pondence between the array 
ot PECANS vsea by the main \ine proqvam 
and the array of WALNUTS used by the 
subyovtine . LA\\so whatever the subrovtine 
does with the CASHEW is really being done 
with the HAZEL (2) in the main.4Ahvs, the 
Names in the subrovtine ave Dummy Names 
which veally stand for the same storage 
locations vsed by the covres pending variables 
back in the main \ine Ppreavram., 

On cther wevds ; the (same) 
bureav Arawer Aces by bath 
namés , by one name jn the 
main and by the other in 
the subroutine subproqram, 


Notice thot the on \ thin the two eae 
have in Common js the actual bureau Arawers 
that they share throvgh the araument list. 
Bither Prcamm Can pot things inte or teke 
then pvt of the common stovaye locatien ) but 
neithey proavam knows or cares what name 
the other program uses tov its vaviables. Not 
only are variable names Totally independent 
between the poo ee but also so ave statement 


numbers. Inhe subrovtine has no idea that the 
main \ine has a statement # lO, and vice vevsa. 
Dt was pyre ly coincidental that the state merit 
numbers were the samé, 

hus , in this case) it wovld seem we 
hada randomly mixea bakch of nuts. 
What's almest bot not quite tyrve. VF the 


Main program had an INTEGER like 
INTEGER PEANUT (5) 


— — (PEANUT, KOKNUT) 
° 


then the KOKNUT wovlda get cyack.ea when 
it rolleA aveund in place of the real GASHEW. 
[Pur ther, you would end up with peanvt butter 
back inthe main program , because the 
subroutine wovld cram in valves thinking 
it was Stovina reals! 

Or hus ) CVE V\ though You Can change ON, 
integer toa real (ov vice versa ) ACVOSS 
an equal siqn,as in 

Y= 
Yov Can't Change types \n qeing threvgh AY\ 


arg ument list, 


What about the 
Name of the Subrovtine? 


Wovu can vse \eaa\ variable name 
a5 a svbrovtine name, because the answers 
are passed back and Sorth through the 
Ce om ent list, Not threv h The prog Yam 
name. Whus the fact that NUTS Was an 
integer name was totally \rvelevant , since 


no valve was ever stored in that name! 


Psst | 
Wanna Know 
a pecret ° 


UNG 


Sa 


IM) ost computers Aow + bother te check 
array AiImensi ons Lor thinas 
argument list, 4fhe program that 
subprogram has already set aside 
storage lecations For any arrays 
the Way momma always arranges 
When AO Main proqvam Uses ae subprogray , it 
really tells the sub progvam where to find 


the drawers Containing the valves jn tts 
argument list, 


Dc, 
When aN Prvray Name appears inan 


argument list of a sub program reference, 
the Main prea ram 16 actva 4 te\lin the 
Subpreqvam where to tind the first element 
of the array / Since the sub program Knows 
how Momma always arranges hey Avawers, 
a\\ it needs to Know to Find a pay ticular 
ara element is 


Poy Where the first element & the 


Array 15 stored 


and 
KF Which element it is looking for, 
Ht then finds that element by conting 


drawers, startina with the one the calling 
program said was the first. 

Hence and forsocth , mu dear Watson, 
caren this impeccable logic to its inevitable 
Conclusion, we start to see the possibilities Gor 
incredible skulduggeru, devious intrigue, 
conniving, underhanded trickery , and, if I 


(3 


may be so bola, SRROR? Your objective, 
snovla vou chocse tw acce pt this Mission, is to 


learn to avoid the 
KOK 
ERRD 


whi le bending the computer to Your superior 
i\\ b lpet v ons, 
Wi 4 emp one all ne? ¥ wiles and weap ns 


PART 1 


Avoiding the E rrors 


les arrays are Aimensionea the same way in 
the Main PYeavam aS in the sub program , You 
ave safest. [Anyone who helps Yov debug Your 
program ana throws your cava deck at Your 
head (Provably), is Syustvabed over some other 
stupid thing Yov ’ve Aone, ins corres ponding 
AYVAUS are dimensioned. the same aiid, in 
both programs, at least the computer will 
follow the same rules for finding a particular 


avvaAY element in both \aces, 
Pov Instance aA 


Por instance, suppose you have an array 
dimensioned 
DIMENSION P2NIA(5,7) 
D+ You asked the computer for 
P2NIA (3,2) 


the Com puter wov la find it Seven locations 


away from P2NIA (1,\). 


dah , where did ya 
et this Seven from 
all of a sudden? 


Seven, my friend 15 one \ess than ejaht. 

IN| ow e\ Wes, historically, is twice four and 
\s the cv be root of 512 which is known not 
to bea prime numbey, [nl owever, more 
pertinent to the case at hand is the fact that 
the P2ENIA ONT AM has five vows to it. Since 
OXYAUS AYE stoved b columns, we. need te 

o down five elements to ger past the first 
column, and then another three to get to 
the thiva element inthe second Column .By 
Lobachevsky’s third theorevn on iwational 
Spaces, this 15 Seven Awa fyom the First 


element. Q.E.D 


IM] ore generally , if You have a two- 
aAimensional OVV AL dimensioned Nyows by 
SOme number of Columns and Yov wanted to 
find the (1,5) element of the arvay, it 
wovld be Souna in the 


I+ (3-1) ¥ Nrow 


th storage location , counting from where the 
Orvau begins. IF or instance, PENTA(,1) 


isin location 


I+ (1-1) #5 
Lea. location |) ana P2NIA (4,5) 1S in 
location 

sae oO ase 


(i.e. ) location 24 ) 23 away from P2NIA(IN), 


How about 
P2NIA (4,7) 
2 


et’s see, Ferhat would be the 
A+ (7-1l) #5 


th element. FPhat?s 3d, Wold it! 
Were past the PENIAs and inte the 
briar patch, 


PAKT 2 


Bending the Gom ute ry Lo Your Will 


Dn the Cinest 
domestic. computers, 
one reserves ample 
storage S pace for 
AVvas4S in advance, 
prey te their actva\ vse, No Bostonian woulda 
have it otherwise. In cther words, when the 
maitre A? com pevtevr asks at the start of 
the pvogram “and how many array elements 
shall we expect will be coming P" Wovr main 
progvam shovia sa 


DIMENSION P2NIA (5,7) 


where the 5 and 7 ave integer corstants 
that reserve all the stovage that the 


P2NIA array conceivably might need. Dn other 
words, the first subscript of PENIA is never 
expected to exceed J and the second subscript 
is never expected to get bigger than 7 

Since Stora ge space for these 35 P2NIA’s 
has now been set aside, O subprogram that 
Wanted to vse that ayrauy would not(itself) 
need to allocate storage for the aVrau. 
Rather, it would just need to know where to 
go to find the first PENIA, and how many 
to pluck. Srhus, the subprogram Could use 
variables in its dimension statements , where 
the actua\ values of these dimensions were 
passed down to the subproqram throvah the 
argument list from the calling program, lou 
could have a subroutine like this: 


SUBROUTINE FLOWR(BUD,M,N, ROSE) 
DIMENSION BUD(M,N) 


ROSE=O0. 
DO 10 I=\,M 
DO 10 J=\.N 


10 ROSE = ROSE+BUD(I,5) 
RETURN 
END 


This routine might be invoked by 
GALL FLOWR( P2NIA,5,7, THORN) 


Un this subprogram, BUD isa Variably 
Dimensioned Avray. Trhe subrovtine Knows 
the exact size of the array because the 
calling Peqvem Pass es the FL OWR is ba 
the dimensions of PeNIA just as they are 


given in the main program. 


MAIN SUBROV TINE 
CALL FLOWR (P2miA, 5,7, THORN) 


SS_—_- 


SUBROUTINE FLOWR (BUD,M,N, ROSE) 
DIMENSION BUD(M,W) 


In the subroutine , BUD isa dummy array 
name, used to stand for the batch of drawers 
in which momma has stored the P2NIA array. 


arhe main saus 
(eek FKL) the es gram Says 


Bo the main Says “PENA, dimensioned 5X7 ° 
and the ore Says, ‘Hi sr2NIA mind if 


; 


ca\\ You BUD?” and the main says 


O.K. , but none of that 
mixed- mode stuff you hear ° 


ss se 82 & 


sé 
95 847 74 «+178 80 


Some of you probably would like to get even 
more deeply immersed in the subject of arrays. 
Crhere ave lots of neat tricks that computer 
jockeys can play by fiddling around with 
subprogram aAimension statements. Sometimes 
they can pv) the wool over momma’s eyes by 
clever\y switching array dimensions from 
program Lo program. Rnowing how Momma. 
aVTANAeS her drawers ) you can sometimes do 
sseful things b breaking the “Rules” and focling 
the computer. tHIwever, you can also get into 
BIG TROVBLE that way , if you don’t know 
exactly what you ave doing. Frill you are 
an expert , you shoulda stick to the 
rules and not try faking out 
the computer, hike they say on 
the hair tonic bottles in the 
barber shops, tnat stu ff ‘ca 


Tl = 


Well, | may not be 
a senius but I’m no 
foo! either. 1’ \\ stick 
to the straight-forward 
stuff for now! 


Say - speaking of simple stuff, is there an 
easy way to write simple functions ¢ 


FAs a matter of fact there IS, Murgatroyd, 
AY STATEMENT FUNCTION is a diminutive 
little function defined by a Single FORTRAN 
Dtatement at the beginning of the program that 
needs the function. Now what cov\da be easier 
than that , 1 ask rhetorically ? 

Gottaminit ? VV show You 3 or 4 of 
these Statement Functions In their entirety, 
complete with expository comment cards thrown 


in: 
C FUNCTION TO AVERAGE YOUR NUMBERS 
AVRAGS4 (X1,X2,X3,X4)=(X1+X24 X34 X4)/4. 
G FUNCTION To FIND AREA OF CIRCLE 
C GIVEW THE RADIUS 
AREA (RAD) = 3.14159 * RAD * RAD 


C SQUARING THE CLIRGLE-FUNCTION TO 
C FIND SIDE OF SQUARE HAVING SAME AREA 
C AS GIRGLE OF KADIVUS RYO 
SQUARE (RYO) = SQRT (AREA (RHO)) 
GC POLYNOMIAL AX&*34+BK 4x2 +CX4+D 
: POLY (X) = AX (AX 4X3) + Bx (XXX) 
+ GxX +D 


Srhese Statement Function Definitions 
neea to qo right ak the beginnin of the 
program that wil] use the functions. They 
go \before any FDS an os 
statements but (after an e or dimension 
statements. You ee os Betaet little 
power hovses almost the same way you use 
conventiona\ function subprograms, For 
instance, you might Sau 


\ AVE = (AVRAG4 (P, @,R,S)+AVRAGC4 
“hCT,U,V,W)) 72. 


"CIRCLE = AREA (DIAM /2.) 


however, since each statement function 
is aefined € com pletely contained Wuithiny 
the program that uses it, there are some 
important ai fferences between statement 
functions and function Sub programs, 


Ferst: oa statement function can ONLY be 
used within the program that contains 
the function. Other programs don't have 
the foggiest idea the function exists! 


Second: 1f a statement function uses a 
statement function itself (like the 
SQUARE function which used the KREA 
statement function ), then the statement 
function that is USED needs to have been 
defined ‘in the program , 


Chitd : the dummy arguments used ina 
statement function definition must be 
non- subscripted Variables. 


Fouthls: if the TYPE (real, integer, ete.) 
is not implied by the first charactey in the 
Name of the statement function, then it 
must be specified by a Type Statement 
above it in the program. 


Lastiybutnokleas tly: the definition 
of a statement function contain 
variables which are not arguments, These 
are called Yarameters of the function. 


For instance, back in the example, POLY(X) 


was a Statement function 1 made up that 
had X as a Dummy Avgumen t, but which 
vsed A,B,C, ¢D as YParameters. These 
parameters are actual variable names 
vsed elsewhere in the program. Whenever 
the statement function is used, these 
parameters will have whatever values those 
variables were given as of that time in the 
preare™: 
Example: 
C STATEMENT FONCTION WHIFF 
 SWHIFF (AROMA)=AROMA* (LILY + NOXI0Q) 
‘i+ 180. ¥ AROMA 


f } 
Mine 
NOX10S = 12 
GAS = 83.2 


f 

30 ODOR=WHIFF(GAS) 
Since the parameters LILY g NOXI0S Were 
already assigned valves, the WHIFF function 
Knows their valves even before statement 30 
is reached. AROMA is aitfevent , however. Ak 
isn't a variable name used elsewhere in the program. 
Jt is a dummy argument or place holder. When 
statement 30 is exeevted the actva\ argument is 
passed to khe statement function. 


Ugh! Give exam ple of 
Ty thagor ean Lheorem 
Statement Function! 


ORK6 Savaw In picture sitting on 
hippopotamus skin, gave birth to twin 
PApooses. ( Papeese?) Squaw who sits on 
alligator hide had one papoose, Squaw 
who sits on bearskin had one papoose. 
iheorem Says,” whe Squaw of the 
hippopotamus is equa\ to the sum of the 
Ssquaws of the other two hides.” 

Crummy joke shovida not be dignified by 
illustrating it by means of Statement 
Function. Better Lo quit ond go on to 
subject of GOMMON STATEMENT S. 


One more joke 
MiKe that and 
Dr Kavfman 

divested of scalp! 


Joke wasn t $0 bad. 
At \east he didn’t 
bell the one abovt 
Indian who went inte 
bank and asked to see 
the Loan Arranger! 


his kind of everyday chit-chat is not 


what is meant b 
COMMON STATEMENT 5S 


In FORTRAN : COMMON refers to anarea 
ot storage twat Momma has set aside for 
COMMON ~ use by anu programs that want 
to COMMON gerit. Normall ; each program, 
from the dinkiest ttle function to the biagest 
mainline, runs its own set of storage locations 
Bach program has its own private Wttle batch 
of \ocations in which it kee ps its own constants, 
variables, AYYAUS, and $0 Forth. No other 
program Can barge mM and put something into 
another program's drawers, unless the other 
proaram has aiven it Access. 


< Dont even ask! 
CL? Dont even ask! 


Yfou wanna Know how one program gets 
access to the locations controlled by 
another proaram , right ¢ Wrong ? Tt doesn't 
matter- Vm going to tell you anyway! 


Controlled 
by 
Subprogram 


Drawers 
2 | Gontrolled 
SS 
s Cex : ro va 
ASP ain ols 
ee 7 i 


“rhe argument lists of subprograms 
are \i Ke doorways through which values get 
passed from program to program. When Q Main 


program vses a subprogram , Say by a 
statement like 


CALL SHNAPZ (VODKA) 


it really is telling the gubprearam where to 
find the drawer in which it keeps the VODKA. 


The subroutine probably starts with a 
statement like 


CUBROUTINE SANAPZ (GLASS) 
where GLASS is a Dummy Arqument he 


subprogram doesn’t need to set aside one ot 
its own drawers for the SANAPZ GLASS 


because it simply uses the same VODKA 
drawer used by the other program. 


Nie Hy hic, thash clear | 
AP Ma 
= 


When you have a very long list ot argument 
or. when you have a group of programs that 
vse the same arguments over and over again 
‘+ can be a real TAIN & AGGRAVATION 
Keypunching up the arqument list each time, 
especially i§ most of the arguments don't 
change from one time to the next. Also, the 
computer wastes some of its valuable time an. 
drawer space shuffling a\\ the arguments bac} 
& forth from program to program. 


0? Get to the point. 
Tm starting to doze off ! 


So each variable in GOMMON (doesn’t) 
need to be included in an araument ist. 
Every program that vses the “COMMON 
storage locations Knows where Lo find all 
the GOMMON Variables, once and for all. 

Bach program that wants bo use the 
COMMON area has a non-executable 
statement in it called —wovld Yyov believe— 
ao GOMMON statement! Erhis GOMMON 
statement goes up at the top of the program, 
just after any Type statements theve might 
be (REAL, INTEGER, ete.) and after the 
DIMENSION statements , but before any 
Statement Functions or plain old executable 
FORTRAN, 

The commonest type of GOMMON 
statement consists of the word * COMMON’ 
followed by a list of all the variable or 
array Names which that program uses ana 
which are to be Kept in the GOM MON AYeA. 
i\ typical COMMON statement might be: 


: CONMON SENSE BOND,DEER, PIN, ALONG, 
STOCK ,TATOR, WEALTH, COLD, GETIT, 
‘2°PLACE, THRU, DANT, SURATE, TALYVU, 
3! NOUN, MAN, NOLIDG , LAW, NAIL 


Storage locations in COMMON are 
assigned consecutively , in the same order 
that variable or array Names appear in 
the COMMON statement. Thus, in this 
example, the variable PIN is in the 
fourth real storage location in COMMON 
and NOUN is storea as an integer, right 
after the veal TALYVU. 

Gince aifferent programs may vse aiffevent 
names for the same variables, it is important 
that their COMMON blocks line vp with 
one another, Suppose a second program 
only needed the contents of the COMMON PIN 
location and the COMMON NOLIDG. To insure 
that its COMMON locations lined Up with the 
other program's COMMON, it wovld need to 
start off with three real variables, ist 
as throw away space fillers. Achen , if it 
used the variable name ‘TACK’ for the 
quantity which the fivst program called “PINS 
this secona proaram might Say 


: {DIMENSION GARBAS(3), DUMY(IN, 
JUNK (2) 

| ;COMMON GARBAJ, TACK, DUMY, JUNK, 
‘INFO 


TAis lines Up wi th "PIN" 


This lines ae Mo one cares what happens from 
with "NOLIDG INFO on, at least jn this prearam| 


Due to the Space fillers , TACK in this 
program wovla be stered in the fourth 
COMMOW location - the same place the other 
program refers to by the name ‘PIN.” Iso, 
“INFO in this program corresponds to 
NOLIDG" in the other. Whe different 
‘drawer ” \enaths requirea for integers and 
reals have all been care fully compensated 
for, so that numbers don't hang ovt like 
dirty socks. 

ry hel p the computer in keeping things 
aligned »it 1s qood to oraanize COMMON 
variables so that COMPLE Variables come 
first in COMMON , then come RE ao. then 
INTE GERS , ana finally , LOGICALS. Of 
course, 1f Your program doesn’t have an 
COMPLEX variables , Yov shovian't qo inventing 
them just for the sake of adding Vari ety to 
your COMMON statement! 


She wrote part ofa 
program and | wrote 


the other part.... 


“and she used different variable names 
than I did for the same thinas. She 
called 3.14 Pi in her par t oP the program 
and I called it PYE in muy part. Wow we 
want Lo Combine ovr eH orts inte one 
biq program, Inlow can we tell the computer 
that her (PL and m FYE; shov\a both 
refer to the same storage location in 


the computer eo 
MOFHENG CPaLo BE EAMERT 
Simply vse an EQUIVALENGE statement 
and Sau 
EQUIVALENCE (PI, PYE, PIE) 


and the computer. will Know that these 
Yhree variable names shovld all stand 
for the same thing in storage. If You 


want to te\}| the computer that her array 
named “GLURP; dimensioned (\O X3), and your 
ae Named “BLUBER" Are reall the Same 
and are to be put in the same A rage 
locations, Yov sim ply need +o EQUIVALENCE 
ann one. element of BLUBER to the Same 
element cS GLUBP, Because of the way arrays 
are stored ‘in adjacent stora e locations 
once you have lined Up one Clement of 
GLURP & BLUBER , all the other elements 
will avtomatica\| fall inte place, a BLUBER 
for every GLURP. aThous Yov covla SAY 


‘ VEQUIVALENCE (P1,PYE,PIE), (N,LL3), 
I C@LURP (1,1), BLUBER(I,1)), (VALUE, XX) 
2: ant a ) (UISE@).1Y); CISTI) 
ne 1Z 
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ach pair oS parentheses encloses a 
list of variables or ONVAy elements which are 
to occupy the storage location. No 
unsubseripted ONva4 Names Can appear in 
an equivalence statement. Instead, you 


specify what is to shave storage with a 
sails element of an AVA and the vest of 
the OV AY avtomat cally aAvogs ‘into place 
because & the ioe MOMMA always ar ranaes 
her Arawers, 


Enovgh on YPQUIVALENCE | 
You’re making 


On 
mountain ovt of a.molehil] | 


With EQUIVALENCE that’s easy todo! 
Yfou can also make tin inte gold, Sows ears 
into silk purses, work into play, ants into 
elephants , night into day and other 
wondrous conversions. Just say 


EQUIVALENCE (MY TAIN, MOLHIL), 
i (TIN, GOLD), (SOW, SILK), (RAR, PURSE), 
2: (WORK PLAY), (ANTS, ELEFN'T), (FARRA 


>) 


3) BELAY, (FAWCE’T, ABZUG), (ANITAB, TANG) 


The compoter must 
need to know which 
variables ave equi valenced 
before it comes to 

them in a program, 


doesn’t it? 


You guessed it | TYPE, DIMENSION, COMMON, 
ond EQUIVALENCE statements al\ need to 
Qo up at the top oF Your pro ram. ihe Y¥eason 
is that these are SPECIFICATION statements 
that ell momma how to aAVWVANAE her drawers, 
Till she Knows what drawers shes going to 
need ; she cant start working on the vest of 
the program, 


x 


dah... it?s ttin 
pretty crowded a 
at the top! 


“Phe (order) of statements in a 
program shoulda be 

FUNCTION oy SUBROUTINE statement, if any 
TYPE statements (REAL, INTEGER LOGICAL ete.) 
EXTERNAL stachements (We hauer't seen them yet!) 
DIMENSION statements 

GOMMON statements 

EQUIVALENCE statements 
\DATA statements (We haven't reached them either!) 
Statement Functions 

ft this point don’t fovaet to put in 

The Program lise £ 


fe\lowed by an END card! 
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Since you mentioned it, 
jst wh at is an 
EXTERNAL 


statement e 


Subprograms often will themselves use ather 
subprograms. 4Phat’s a perfectly normal course 
of events and it doesn’t requive any EXTERVAL 
statement 44 ana 44 /\00 5 7% of the time. 
Once in a blue moon (hike the one on the deor 
above ) you \I find that Yov want to write a 
subproavram which uses one subprogram the first 
hime it is called and uses a aifferent subproqram 
the next time it is called. Mlaybe at the time 
you are writing this subprogram, you might not 
even Know what svb progvam you are going to 
want it to use ! Gon fusing yeh wot ¢ 

Actually it is very simple. Quppese for 
instance You want write a subroutine to find 
the roots of a function. Df possible, you 
would \ike your subrovtine tobe versatile and 
general purpose. when you can Save it as part 
of your own personal program library. One 
Lime you might want Your root finding program 
to Find the roots of polynomial , another 


time you might use it to solve a 
trigonometric equation, ana so forth. Dn 
other words, Yov want to write the root 
Finding program in such a way that the 
Main proavam can tell the sub program jest 
what ae Sees function to use while the 


proqram is running! 


How can it do that? 


By passim the Name of the particular 
Sub program to be vsea down throvah the 


araum ent list : 


dah... could ya give US 
an eEXam ple? 
a 

Sv ppose Uovr root finding sub program 1s 
named “ROOT” and it knows how to find 
voots of Some function namea FUNK’ The 
root Sinder might. be something like 

| SUBROUTINE ROOT (FUNK, A|B,ANS) 
Deviovs 


Intrigue TE ( ABS (FUNK (R).LT. TOL) ANS=% 


; 
RETURN 
END 


We will assume that this routine 3$ Filled 
vith some DEVIOUS INTRIGUE by means £ 
which jt finds the root cf FUNK. TH 
returns the root throvah the dummy argument 
‘ANS? 

The ROOT sub program Knows thal FUNK 
is © Function name. ow does it know, 4ov 
ask® Tt Knows becavse it realizes that 
FUNKE covidwt be an array name because 
4 isn’t dimensioned anywhere . Porther, since 
the name "FUNK" appears in. the argument list 
of the subroutine, the subprogram knows that 
FUNE js a dummy Function name ana that the 
actual subpre ram name will be passed throveh 
the argument hist when the svborovtine Is 
actually called. 

irhe MAIN proqvam that vses this ROOT 
Subprogram might SAY something hke 


CALL ROOT (CUBIC, START ENNDRESLT) 


when it wanted KOOT to Vse the program 
named CUBIC” in place of the dumme name 
"FUONR’ and it might say 
” $ 
CALL, ROOT (TRIG,ALMA BETA,VALVE) 
when it wanted ROO* to use the “TRIG 


program ‘in place of FUNK. Yn other 
words , FUNK’ was vsed asa Avm mi sub program 
name which SCEYVES A285 A place holder inside 
ROOT. One time ROOT vses the CUBIC 
function and the next time iL vses the TRIG 
Function. 

Unfortunate) J the Main ae has 
no way of Knowing that “CUBIC” and ‘TRIG’ 
ave subprogram names and not just. plain ola 
variable names. 4Phe main program might he 
tempted to lead you astray — but wait! 
ere’s the TEPER, se en! to 
clear Up all the confusion bb explaining to 
the main that ‘CUBIC “and TRA" ave Subprogram 


AMeS and not Variable names! 
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In order that our mild-mannered PATERNAL, 
stakement not arrive tao late to save the 
situation, it has to be right vp at the top, 
just after the ty pe state ments. IMere’s how 
it’s done 

REAL MONEY 
EXTERNAL “— , CUBIC, ALEVIL 


CALL a a ce ine 


CALL pre (TRIG, ALRUD BETA , VALVE) 
CALL foot eerahanneers MONEY) 


rhe EXTERNAL statement tells the 
program that contains it that CVBIC, TRIG, & 
ALEVIL. are subprogqram names and not 
variable names. Phe only time you need an 
EX TERWAL statement is when Yov are Na 
program that 1S trying to pass a subprogvam 
Name to another subprogram. The externa] 
statement 02s after the type statements 
and before the DIMENSIO statements in 
We program that is Acing the passing. hus jin 
the Kample, ALEVIL is “a function name bei 
passe 4p ROOT and MONEY is the ROOT of ALEVIL. 


Like the farmer in the 
corn field said, 


L hope this book 
is about over. Data 
\S coming ovt of 
my ears! 


Don't get cornfused, Cornelius! ihere 
are other ways to store your data! VF you 
are constantly Storing constant constants, 
you'a do better to Use a DATA statement. 
bet me explain. 

So far we have studied two Ways te put 


constant valves inte storage. One way 1S by 
o READ statement like 


READ(8,10)X,M, (DEL(D, 1=1, 30) 


Bach time the read statement is executed, 
new valves get vead in for X,M, and the 
30 elements of the DEL AYYay . oro change 


these valves, all You have to do is change 
the constants on the data. cords. You don't 
need to change the program itself, 

Whe other Wau we studied for putting 
constants ‘mto storage locations was by means 
of Assignment Statements hke 

PL= 3.14159 
Srhat’s a perfectly 000. way to put the value 
3.14154 into the location named Fl, but it 
\sn’t necessarily the best Way. Suppose you 
had that statement inside of a subprogram 
that vsed the variable named PI. Bach 
time the computer comes to that assignment 
statement, it takes the constant valve 
3.14159 and stores it in the drawer nameaA 
PI. Of yov use the subprogram 50000 
times in a row, iL does that 50000 times. 
[Pretty dumb, huh? Jt mau not take long 
to do, but it does take some time and 
little bits of wasted time can add UP) 
especially inside of a loop. Why waste time 
doing something like this overand over again? 
Aster all ) some constants are more constant 


than others. trake the valve of Dy ue 


for CXaM ple. {E ver since Christoph ev Golumbs 
dropped the apple on \ssac Newton's head 
proving the earth was the center of the 
universe , most scientists and the majovit 
of leais\ators have been willing to Keep 
PT =3.149159, ven Bert Lance Uses P1L= 
3.14154 although once or twice he was caught 
Using Pi= 4. 

So yov ask yourself, \ Self. (low can I put 
dota valves into storage before the program 
starts to run by a non-executable statement 
that won't get executed over and over while 
the program 1s running?” And the answer 


You dive yourse\f is, ‘Darned if 1 Know.” 


Since you don’t Know 
the answer, Vil answer 


the question! 


‘rhe answer is Yov USE Q@ 


DATA STATEMENT. 


stora e 5e fore » the YOQYam starts executin 
on Raa SS Brome Yow on , the dala J 
statement 15 ignored , $O it Acesn t take 
any more time. Wlere’s what a data 
statement might look like: 

DATA PL LAST, ELOG/3.\4154,56,2. 71828 / 
“This sim pl means witialize PI to 3.14154, 
LAST to 5b and HLOG to 2.71828. Notice, 
theve is a list of variables to be initialized, 
followed by a \ist of imtial values for these 
variables.aphe list of values is suvvounded by 
a& pair of slashes to set it off from the 
list of variables. Is You would expect , the 
First val ve gets assigned to the First variable, 
ana so forth. 

[Avra elements Or AYTA4 Names Can also 
be initialized. Valves get tucked away in 
aa mM the Same oraer that the QAYYQ 
elements appear in storaae , that 18) in co\umn 
order. When severa\ consecutive elements 1h 
ao data statement OC t, be initialized ty 
the Same constant value , the constant can 
be preceded by OW repeat constant as in 


30% 5.6,-18, 10% 0. 


Srhis means there are 30 consecutive values 


of 5.6, then the integer cons tant-l8, and 
Finally ten real constant valves of 0.0. 


irhe genera Form of a DATA statement is 


{DATA list \/ valve list 1, list 2/value 
Mist 2A \ist 3/ value list 3/, etcetera 


lAlere , the \ist is list ey and so forth aye 
simply lists oS variable names ) array elements, 
or array names. Urhese are the names of the 
variables You want to initialize. (Incidentally , 
a data statement ima subproqram can't be 
used to initialize dummy arauments of the 
Subpro ram since storage for these ar uments 
is under the Rightful “dSurisdichion of some 
other program.) Bach \ist of names is folloued 
by a s\ash.Arhen comes the list of values in 
One to One Covrespondence te the elements in 
the name hist. Finally comes a Closing slash. Vs 
there is to be another list & names} separate 
it from the preceding slash by means Fa 
Comma. F\\so USE A libera\ amount of COMMAS 
to separate list elements. (Dont be too liberal. 
One comma between Cach paiv of elements is 
just the right number! ) 

hast but not least, stick the DATA statemait 
provaly UP at the top oS the program ,aSter 


the, Bquivalence Statements (which means 
aSter the Common Statements (which means 
ofter the Dimension Statements (which means 
after the External Statements (which means 
after the Ty pe Statements (which means 
ofter the Function or Subprogram Statements 
»)))), 14 any, That way Momma will tuck 
OWO the Anitial Valves when She is 
settin up her dvawers and before she 
actually starts to execute Yovr program. 
Arhus, a proaram might start with: 
DIMENSION MARAY (2,3), ARRY (25) 
: {DATA PLRAD,DEG /3.14154, 00174533, 
1!57.29578/, ARAY /25 XO, /, 
235,X,Y/50, |4.2,-3.85/, 
3: MARAY /3,2,-1,7, 6,5/ 
Chis MEANS inti alize Pl to 3.14159, RAD to 
0,0174533, and DEG to 57,29578. While 
you've at it, store 25 zeros as startin 
values for the array named ARAY. Don't git 
ut give JJa starting va\ve of 50, give 
X a starting value of 14.2, and give Ya 
starting value of —3.85. One thing move - 
initialize MNARAY SO that 
MARAY (1,1) = 3 MARAY (12) =-/ MARAY (13) = 6 
MARAY (2,1) =2 MARAY(2,2)= 7 MARAY(2,3)=5 


There must be 
a Zillion things 
to remember 


about FORTRAN ! 


Ut?s not really that bad Most of it will 
come naturally Lo you once you get the hang 
of it, 

But if you want to use numbers \ike a 
zillion, Yov really shovld Know about & 
format and the exponentia|\ form for rea\ 
constants. 1 didwt tell you abovt it before 
because you had plenty of othey things to fret 
obovt, 

i\s yov indubitably know, very large ana 
ver 4 sma\\ numbers ave easiev to manage when 
yov use scientific nobation. Instead of 
writing | 

8875190000000 0000 C00 CoC Coco, 
it’s easier to write o7 
387519 X 10 
ond instead of 
0,000 00 cccc cc Co0000 Gee 00 00000642 
it’s easier to write 


692 x10" 


TPhat’s how you ho it in algebra, tow do you 
to it on the computer © 

@n the computer you can either write a 
real constant the wa yov learned back at 
the start of this learneaG scholarly work 
or else Yov Cay) write ik in exponential fom, 


In E form, the constants You just saw wovla be 
3.B87519E27 ond 6.92 K-24 


Of Covv se, Uov coulda also write them aS 


387.519%25 and 0.692 7-28 


In this notation, You have a real constant 
followed by the letter Ky and a one or two 
digit integer exponent.The exponent Can have 
a. plus or minus sign in front AFhe exponent tells 
what power of ten 1S multiplying the rea\ 
constant. 

Hleve are some more valid rea\ constants: 


In fact since the power of ten exponent 
allows you to effectively "float" the decimal 
point back and forth inside the number, 
FORTRAN also allows Yov to write real} 
floating point constants in the form of an 
integer constant followed by an Kanda 
one or two di it exponent .4rhus, the real 
constant 3.)4159 could be written in 
FORTRAN as 

3.14159 or 3I4I59E-5 or O314/59E| 
or 3.14159 20 or 31.4154 E-01 or 314159 E-05 
or 314159.E-05 or any of a zillion cther ways, 
These are all considered valid Real Constants 
even when written without a decimal point as 
in 314159 E-5".4Phe veason is that the E 


exponent implies the presence o& the decimal 
point even thovgh it doesn't appear explicitly. 


doh.... low do 


Ya write ovt 
numbers in... dah... 


BE format ? 


BriMiant, mu boy brilliant | AFhat’s yust. 
how you write them ovt. Yo v vse EK format 
and you SpeciTY it by means of an 
K field specification code, as in 

8OFORMAT ("SDISTANCE TO TWINKLING 
iE STAR ISb',E10.3,bLIGHT YEARS.’) 
Ihlere , Ine” \O"means “the next | O columns 
are going to contain this number.” She *.3" 
means “write out 3 Significant digits.” (If 
the number being written ovt was 
184395200000 


it wovld come ovt on the printer as 
DISTANCE TO TWINKLING STAR ISb b0,164E bI2 bL 16... 


Prhe number ~483759% lo* written in 820.5 
form woulda come out 
bobbbbbb-0,18376E-06 
tt ee et ttt 
a eld le el 


7 More Columns! 
20 Columns Total 
Wlotice you need to allow at least Gevemextra 
columns, in addition to the number cf significant 
digits you want to print. Fhese are needed for the 
sign, exponent, decima\ point and leading zero, 


And in conclusion, 


Td like to say... 


there’s lots more FORTRAN we covld 
discuss, but Now YoU Know plenty Loy getting 
started . l\s you ge along Yov Can easily 
pick vp the rest. 

FA\\ 1 can do is tel You some of the 
rules & principles of FORTRAN. When you 
first see them all in a heap seem 


mind - Yfou think to 


LA\ ctva!| J it’s not a\\ that bad. Dont 
try to \earn FORTRAN by sitting down 
and memorizina this book, Someone will 
have you Commi tted. Yfou ll be babbling 
incoherently within the hour, LA\ mazingly 
enough , the same Thin g is trve + You take 


that approach with a formal fORTRAN 


mM anva\ I 


So if we arent supposed to 
plow throvah this book what 
aye we supposed to do with 
ick? It?s too late to get ovr 
money back. 


orhe right Wau bo vead a FORTRAN 
book 1s bu a series of skimming passes. 
Bach time through you pick up a bit move of 
the nitty gritty derail. When You pick Upan 
IBM. manual ) for instance, first f lip through 
it locking for the jokes. Where arent ANY , 
So 90 back and flip throvah again qetti 
Familiar with the overall idea. ‘When flip 
throvas fvom back to {vont , side to side, 
ond top to bottom. If theres a center fold, 
yov ve m the wyon publication! Bach time 
Through , you \\ be looking for specific details 
in veaard to what uov are doing at the 
time on the computer. f\fter a while, the whole 
works will have seeped throvgh your skin by 
Osmosis, 

Didn't they take Osmosis 
of the market ? 


frwo thinas will motivate You to \earn 
al\\ the details of FORTYRAN. 


ll One. A\ sincere desire to 


pass this course, 


“= Two. Whe fact that Programming 
acomputer can be a hell of 


a \ot of fun. 


Yov learn how Lo program not b 
reading about it or \istenina to lectures on 
programming. You learn to proaram by 
writing programs, Yfou \earn how (hot. 
to make mistakes by makin mistakes. 
Each Lime Yov make a grammatical boner in 
EOSIN, the com puter will, with infinite 
patience, point it ovt bo You, At wil even 
go $0 Far as to tell You the Breed of the 
Blunder. [Being of sound mind and body but 
having O Finite Vondness for Futile Frustration 
ca maKe O mental note not to Gommit that 

in st since 7 Was totally lacking iy”) 
Compensatory pleasures. 

n the process of correcting Your mistakes 
your fingers will do so much walXjn’ through 
the FORTRAN manuals that the pages will 
turn yellow, Before You know it though, you'll 


you've pubconscious\y — a|\ 7 
idd\in re ett altr articulars 0 
ORTEAN Soh c Asevie Sl wake up and 
find yov ve Aone from Tonoramus to 
Intelligentsia and are enjoying it to boot! 
Once You reach that stage, it will all be 
second nature to you, [t's almost like riding 
a bike. Once You've mastered it, UOU Can 
come back to it years \ater and it will 
still be etched in Your memory. Vou mau 
think you ve Sorqotten it but aftera 
brief review it will all come rushing back 
to you, assuming You really learned it in 
the first place? 

Srhe worst mistake ad can make mn 
learning FORTRAN is to be ofvai nd to bry 
Some thing on the com puter for feav of 
making a mistake. Sure, You shov\a try to 
make Yovy programs as acevrate as possitle, 
but don’t be paranoja abovt it. Mistakes 
are @ natoral part of \earnin and iF you 
Wwory too mucn abovt makin mistakes Yoav 
won't ever stretch your nevrons by trying 
any thing new, 

Im going to quit now, since my violin 


string just broke. 


maiple Problems 


EXerCISeS 
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Wome -~Ganned Pickles, Preserves, 
And SQUARE ROOTS 


A Problem You’l\ Relish! 
Suppose, st to live life the hard way, You 


wanted to write your own computer progvam 
to calculate Square roots, Yfou might start by 
observing that if 


y= VX, then y= XK or yoo 


yu | . // 
This last equation gives you a. “brilliant ~ idea, 
Wyo find the square root of x, you simply need 
to divide x by its sqvave root. Phe result 
is Nx 7 Joyfully you picture this in your 


mind: 
) oars J 


LA \\ you need is a machine that wil\ divide 
a number xX by a number u . he number 
that comes ovt wi \\ be the Square root & x! 


Wait a mnvte! Where will 
you get the number ye 


Ovt of the machine | We?) ca\\ that Y P 


EXAMPLE I-1 


Svdden\y You are chicken. You need to 
divide x by UY in order to get the answer, 
But y is the answer, so you need the answer 
St to divide into x, but you are eqged on 
bu wh WY comes ovt before ib went in’? 
Wh , oh why, O Ys 


Bolzano’s bambino tells You 


You are usin Y 
in its own oP niton, 
Mathematicians 
call it an 


"Tmpli at Function: 


DWicit you covla relate to, but implicit yy, 
doesn't help you at a\\. But you etan idealys 


4 GUESS 


What will come out ? (Shredded wheat? 


No! If the quess 1S \ess than VX, then 
what comes ovt willbe more ! And if the 
qvess 1s too big , wnat comes out wil\ be too small. 


EXAMPLE |-2 


Qo? So it’s a\\ sewn up! Whatever l 
gvess, 1 Know the trve Vx lies 

between my quess and whatevev comes 
ovt! Othese two numbers are Upper ana 
lower bounds on the trve Vx. 
So if I average those two numbers ; 
T\\ have an even better quess. 1 can 
again and again till Im as close as Ll want 


to get! 


Wait a min ute! 
It would take you 
forever to find 
the exact root 


And the computer only 
Keeps abovt Six figures 
of accvracy in its 

calculations, so even if 
ov \et it run forever, 
it still might not find 


the exact answer. 


So who needs to be that precise ¢ 1 \) 
stop when I’m within a reasonable tolerance 
of the true answey, Hleve’s what Wl do: 


EXAMPLE !-3 


~H Bs acl 
> y ss 
Read Number = < € = 
of X?’s a ae 
3 ‘3 > -~ 
> ti. 3 
Pia Se 2 
Read xX € tol = Su 3 
ES ae 
“Se ol 
a 
er s & 
Sl Tt yl 
"ia Perel 2 
0 
Initialize = ui GAD 
+ Ay © -< 
guess for \ S es 


Divide gvess v 
into xX ” 
No A\\ 
Average to Done 
get new quess y 
No 
Yes 


Calevlate 


EXAMPLE I-4 


Z\nd here’s my actva\ program . 


11 JOB DR. K These control cards will be 
// YOR different on your computer! 
C SQUARE ROOT SAMPLE PROGRAM 
C READ NUMBER OF CARDS WITH ROOTS TO FIND 
READ(B,10N 
lO FORMAT (15) 
C ECHO PRINT THE INPUT DATA 
WRITE(5,20) N 
20 FORMAT C'ITHE NUMBER OF ROOTS IS’, 1I5//) 
C READ NUMBER WHOSE ROOT 15 DESIRED 
C AND TOLERANCE, REQUIRED 
30 READ(8,40) X, TOL 
40 FORMAT (2FI0. 4) 
C ECHO PRINT TO BE SATE 
WRITE (5,50) X , TOL 
20: ‘FORMAT ("bX=b! ,F10.4,5X, 
‘1: "bTOLERANCE ISv FIO. 4) 
C Grr Y GUESS 
YG=xX 
60 TEMP=X/YG 
CALCULATE NEW GUESS FOR ROOT 
YG6=0.5*(YG+4 TEMP) 
CALCULATE THE ERROR- CHECK AGAINST XIN 
CASE OUR THEORY ITSELF HAS A BUG 
ERROR = ABS (X- YG X YG) 


(796) EXAMPLE 1-5 


C LET'S WATCH THE ALGORITHM CONVERGE 
WRITE (5,70) ERROR 

70 FORMAT (15X,"ERROR 1$b’, El2.5) 

C LOOP BACK IF NOT WITHIN TOLERANCE 
IF (ERROR .GT. TOL) GO TO 60 

CG PRINT FINAL ANSWER 
WRITE (5,80) X, YG 

80 |:FORMAT (/'bTHE SQUARE ROOTOFD’, 

‘EFIO.4, "bISb’, F15.6//) 
C ARE THERE MORE ROOTS TO FIND ? 


N=N-1 
IF(N.GT.0) GO TO 30 
WRITE (5,40) 
40 ::FORMAT( //"bTHTHTHAT"S ALL , b 
FOLKS") 
CALL EXIT 
END 
MoXEQ 
bbbb4 
2, 0,001 
7.62 : 0.0005 
25. ; 0.0001 | 
8.36 0.001 : 
/END 3 


Columns 1/-20 
These cards Mag di theron Ovr 


computer system, so Check: 


EXAMPLE I-6 


*Tn Sooth, 1 Know Not 

Why lAm So Gad. 

*Tis Ascribed To All This 
Homework ! 

A Pox On The Cad” 
The Merchant of Venice 
Volgate Watergate Edition 
Antonio Lé.1 


Nicomach US, Known to his friends bach in 
\S* century Athens as “Nicomachus,” wrotea 
book (using ao cold chisel), In his 
“Tntroductio Arithmetica” he asked the 


musical question, 


How can the cubes be 
represented in terms of 
the natural numbers? 


LPs Greek to me! 


His answer was that 


Cubical numbers are always equa to the 
Sums of successive odd numbers, % 
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For instance 
:* \=| 
2° = B8=345 
3° =27=7494\ 
4° =64 = 13415 417419 


d 


; : a : 
n= P= 3 bp 


where m, is given by 


and the other m’s are the following 
oda inte ers, That 15) Me=M, +2, MZ= Mare, 
and so forth, 

Note that the sum of the first n-1 
integers is given by 

> iL ._ nn (n-D 

L=I e 

Write a flowchart and program using 
integer orithmetic to print ovt a table of 
numbers and their cubes, as caleviated b 
Nicomachus’ scheme. apey to print out the 
table for valves of n from \to 33, Something 
funny may happen at the end. See if you can 
Figure out what! (Then too, nothing may happen!) 


PROBLEM }-2 


Speakin Q ot 

VAN O players 
is minds me os 
another qveat 
mathematician, 


Fibonacci! 


a, 


Everyone Knows Liberace, but did 
You know thak in 1202 Fibonacc) wrote 
the Liber abaci’ or Book of the Abacus?” 
(Actually Fibonacci was sort of his 
nickname. le Was also known As Leonardo 
of Pisa.’ Scholars arent svre what Fibonaci 
means. Some translate it as “Son of aood 
nature,” others as “ Prosperity,” but the 
most popular interpretation 15 “Son of an 
Ass.” TD think that last translation is about 
half wi ant » 

[Be that as it AM little Leonardo once 
put his fertile mind to work theorizing how 
rabbits mu\ti plied. Hie figured ib had 
nothing to do with theiy fertile minds. Whe 
way he Tiavved it, it takes rabbits about 


PROBLEM 2-1 


ao month to mature From then on, he 
conjectured , O pair o& rabbits cov\d brin 
forth a \itter everu month ‘ro keep things 
simple, he assumed oath lither woulda have 
exactly two rabbits, one female and the other 
the ot her: no rabbits ever passed away 
since they Wwe re havin such fon 5 and once 
rabbits were mature, they wovld do what 
came natovrall every month, naturally. 
Suppose we \et T,, be the number & 
pairs of rabbits alive at the end of the 
kt month. hen, the way Fibonace} 
Figured or sly ttle devil that he Was, 
15 yov started with one pair of newborn 
rabbits, Yov made a bj mistake already. 
Por at the end of the fivst month, you'd 
still have a pair of rabbits, but now they'd 
have a gleam in their eyes. At the end ef 
the $e cond month , you A have anice new 
Wtker of newborns. But you'd also have 
OW pair o& matvve advits makina whoopee | 
Wext month , they’ d have a second \itter of 
newborns, but their first litter would be 
Yeady for a little orgy of thei own, 
Bach month , all the vabbits who weve 
axound last month are still theve. Boutin 
addition , a\\ the rabbits who were avound 
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two months earlier have babies | Suppose 
each picture & represents the number of pairs 
alive ot the end & the month. Phen ovr 
month Lo month “Planned “Parenthood chart 
wovld have the following hare-brained setup: 
Month F Rabbit Pairs Alive 


EGU 


YUM WZ UM 2X5 EY 


34 Kelp | Rabbits are 
29 dropping ofS the edge 
2 ah of the page. |’m 


So DON DN AWN — 
DOU wr ——/| 


drowning in droppings! 
Notice that F, =F =| . For the k® 


month, we have 
¥, 7 Fe t Fez 


Orhis “Recursion Relation allows us to 
calculate next month's rabbit population, 
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Knowing the current month's population 
and \ast month’s total. 

OF course, most of you find this deadly 
dv\\, since You don’t plan to raise 
rabbits ina box under Yovy bed .ihere’A 
be a terrble odor and uou aren't svre the 
rabbits covld ger used to it, Nonetheless, 
Fibonacci’s Aiscovery of how rabbits 
procreate laid the aroundwork Foy 
much of the subject Known today as 
in finite series, So it’s valuable to know 
even 1f You arent a zooKee per. 

With that background , write a 
computer program to generate a table 
of the first 100 Pibonace numbers, It 
shovld start ot as \,\,2,3,5,8,).... 

Be sure bo number the numbers when You 
write them out, so you'll Know that 8 was 
the 6” Fibonacci number. 


What’s the 


\2 th 
Liberace number? I think it’s 
‘Chopin by 
Candlelight” 
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Ph. \D), Wiceiex 
A led Higher & Deeper 


Uf one inseribes a polygon of n sidesina 
circle, the perimeter of the polygon will be 
less than that of the circle.Ifone circumseribes 
a. polygon of n sides about the circle, the 
perimeter of the polygon will be greater 
than that of the circle. 


How do you cireumsize a polygon? 


The werd is “cireumscribe, 
Tdiot! 


aro 1G decimal places, the perimeter of 
a unit circle (one whose diameter 16 One 
and the same as one) is 3,141592653589743 , 
o number somewhat close to TY 

Develop & flowchart & program an algorithm 
bo aetermine upper ana lower bounds on 
the trve valve of Tr Lo within any 
specified accuracy. Do this by first 
mseribing ana cireumseribing an equilateral] 
triangle (n=3), Galev\ate G print their 
perimeters. Phen check whether their 
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difference 1S less than the allowed ervvor, If 
ib is not, then double the number of 
siaes and repeat the procedure, 

Include yovy derivations ana fi qures 


with the algorithm Slowchart. 


When shovla we quit? 


Quit if You have more than 4b sides, 
oy \f You have reached a tolerance of 


better than 0.0005 . Or hen sit back and 
study your ovt put. 


FFE Took Some a a 
mean 
Pi got closer ana closer for 
a while and then it started 
getting worse! And the 


perimeter of my inscribed 
polygon got biager than 
that of my cwcle! 

How deo You explain that? 


Wow do you © 


————e 
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One for the Birds... 
Z Heron’s Formula 
Wp Hy) a AN AN AAAYSY 


A~_A_AAHANAAYWY 


ne L\ convex, non- requiay 

—E Sey YD polygon with M ver tices 

ee is specified by Giving the 
X,3 coordinates of each 
vertex: 


Floweh art , write, & run a Program to 
read in the vertices of a sequence of 
polygons, print the vertices, and calculate 
and print the total area of each polygon. 
Use the following method: 

Subdivide each polygon into M-2 
triangles, all radiating from the first 


PROBLEM 4-1 


vertex. rhe area of the polygon is 
A= 2. ae 
\= 


where A, is the avea of the es triangle, 
he area of On triangle 15 easily 
ca\cv\ated by means of Heron’s formula: 


A,= NS(8-ay(8-b)(5-0) 


where S=y (@+b+C) and where a,b,c 
ave the \engths of the sides of the jth 
triangle. 

Meron : of Course, 1% best Known for his 
vnvention of the Heronbone suit, Inlis friend, 
Mero of Alexandria,invented the Submarine Sandwich, 

Bince this is all Greek to you, vse the 
Pythagorean Theorem to find the distance 
between vertices? 

(Xz 542) 


a 
(X,,41) 


A=" (X,- XY + C427 yy 


Wlotice the first triangle is bounded by 
vey tices Gyo. the second by vertices 
Osun the kth by vertices | Jet, det2 5 and 
the last by 1, M-1,™M. 
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Read in one vertex per data card. Phe 
last vertex is recognized because itis 
the Same as the first, Pest Your program 
on the following Figures ; 

Figure i Figured Figure 3 Figure 4 
X KX 4h ¥*%)Y X 


2 Re 
QOMLROL 


Figure 5 Fiaqure ¢ Fiavre 7 

So Sh eee 3 

Yov wnvent! Likewise! X 3 
-10,28) 16.15 
+10.28)-16.15 
—10.281 16.15 


Did your program handle *3 and 
#7 O.K.? TF so, gooa for you! If not, 
let that be a lesson to you, Be ready 
for anything and be sure to rogram 
sensible provisions fov properly handling 
SPec\ al cases that m ght arise ! 
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AU, 8. Govt. Inspected Prime Quelity IProblen 


EN prime number is a number Aivisible 
without yvemainder only by the whole numbers 
unity and the prime number itself, Hlere’s 
oa flowchart for Oo prime number algorithm. 
Program it vp and Cnd all primes between 
lana 200. 


Print: PRIMES FROM 


Sneaky very Is 
sneaky ~ 


Next 
exactly divisible 

by ~Itest? 
Next / Itest xItest 


First 
Prime 


, 


Try next # 
(Odas only) 
Next=Next+2 


PROBLEM 5-!] 


‘When Sorrows come, they 
come not Single spies, 


But in battalions. 
Hamlet WV, 5 


InJaving bu\\dozed throvah one program 
to calcviake prime numbers, try it again 
with Dtule & Grace or perhaps Mildred. 

nyway y notice that You coulda write a much 
more efficient Prime Number progvam by Using 
AYYays, You covla store the primes you find 
aS you 40 along, by putting them in @ 
one -dimensiona\ aArvau. irhen ‘ to seeifa 
positive integer 1S prime, you only need to 
divide it by the primes you have already 
calcvlated. It is Ridievlous to waste Lime 
dividing by numbers \ike 6, 1,12, or 4, once 
You know they arent primes! 

Write a program to find the first 75 
primes Using ONY AMS, iSemember: \ 1S Not a Prime, 
Yov May assume 2 is prime, but figure the rest ovt 

Pot a nice heading on Yovy printout ana then 
write all 75 Primes USing at Most 5 additiona) lines! 
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Xssac Newton’s Very Own Method for Solving, 


Wonlinear Equations 
Suppose You had a function Y= (Xx) and 


ov wanted to find its roots, 1.e. places 
where §(X)=O. On the picture, XH is the 


unknown root we are looking For, 


4 


P(x, , x) 


If ovv Intrepid Root Finder starts by 
guessing that the root lies at a point X,, 
he will find that he isn’t at the root at 
aM, but is standing On a point P way 
Up on the curve. Tlowever, notice that the 
Tangent to the curve at P cuts the Xaxis 
at a point Xe . Note also, that Xz is much 
closer to & than was ovv startin Qvess Xp. 

IS we can ealewlate the Slope of the 
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Tangent ax PB, then we can calevilate X,, 
We cov\ah then repeat the Process , USin Xo 
in place of X, Etcetera <iisiabeeasiinian 
After a while, if we ave on a nice, 
beginners Wi) and not a bump expert's 
slalom Course, we would find we had zoomed 
in pretty close to &. When we saw that 
OvY £ (x) Was sufficient clese to Zero, we 
cov\d quit, vashl y assuming we had Cound Xn, 
Since the tangent ak FP is 


y-F(X)= £7 (KX) [ X-x, ] 


it 18 obviovs (as they Say whenever anything 
Is obscure) that 
(X,) 


Xo Xx f(x) 

¢ You vse this as Your basic iteration 
scheme, you'll zoom jn \iKe gang busters, if 
the curve doesn’t have an 4 Hips and bumps 
between you and the root. It the correction 
term has n zeros after the de cima 
point, then the resvit is generally good 
to abovt 2@n decima\ Places, 

Incidentally, a 900A feature of the 
Newton method 1s that minor errors wi\\| 
corveck Xhemselves. his is nice if yov are 
vsing we by hand, 


(2) PROBLEM 7-2 


1 

Nothing mn education 1s SO astonishing 
as the amovnt of janorance it 
accumv\ates in the form of 


inert facts | Henry Adams 


Arve, Wherefore try the following 
home work. (Snicker, snicker ! ) 7 

Write o FORTRAN proaqram VSing 
Newton's method to find the recots ff a 
function in the Vicinity of some Given 
starting Qvess. Your program shovld use tup 
Function subprograms, one called FUNK and 
the other called DERIV. You will also be 
wr} ting these two function subprograms, 
One will contain the Function whose Yoot is 
desired and the other will contain the 
derivative of hat funckion. 

Got that ? OK. Now vse Newton’s method 
to find the three roots of the polynomial 

¥-8X-4=0 

One root 16 near —-2.% one is near O, : and 
the last is between 2 and4.5. Find the 
roots to within four decimal places of 


ACLU Yacy e 


Verity each root that Y ov find by 
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checking to see if it really does satisfy 
the polynomial, Inave the computer calculate 
R= X-8X -4 
I? shovld equal zero fer each root, but 
there wi\\ actually be some residve dve 
to small vround-off evrors, truncati on CYvors, 
and the fact that You quit before the 
compute v died oS exhaustion. Whe smaller 
the valve of K , the better the voct. 
Gheck all three voots and print ovt 
both the root and the residve for each 


one, 
YW 


Want to try some thing else ¢ Write Qa. 
program Using Newton’s method to 
determine square voots ana cube roots of 
numbers, Prove Your program on a list of 
at \eask 1O numbers , including 64, 10, 0.3, 
and 3.)41549 . Gompare your results with 
those given by SQRT and ¥%* (1./3.) 
Print ovt atable showing the number: 
ibs Square vook by Newton!s methed } SQRT'S 
Squave root} the Aiference (ink Format); 
the cube Yoot by Newton 5 the computers 
cube root 5 ond their Aifferen ce in Ly Sorat, 
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Yow does Newton’s Method 
have any thing to ao with 


square roots 0? 


Hint: het 
F(xy= A X° 
where A is the number whose root you 
want, 


ChE 
\ddea Little Brain Weaser . 


orhe fiest sample program we saw 
(Pickles — reserve $s, and Square Roots on Page |12) 
also calcv\ated square roots, For $64,000.00 
(oy maybe some extra credit), is there 
any connection between that method for 
Tindina SQuare roots and Newton's 
method for Finding Squave voots? 
What is the connection 2 Ws it French or 
via Peoria? Gan uov get there non stop ¢ 

[nei dentally that otner method for 

Tinding Sqvare roots \S Known as Mechanics 
Rule. Oo SUp pose ve was named after Dlusees 
S. Mechanic , the Samovs actor who Coy Years 
played the role of RouRogers horse Trigger. 


PROBLEM 7-5 


How can I find the area 
Under this curve ? 


E\n easy way te calcviate the area under 
& curve Y= £(X) , from ak XS b, 


is to subdivide the area under the curve into 
vertical strips of equal width. The top of 
each strip can be approximated bya straight 
line as shown below: 


J 


QO ath b X 
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arhus, the area under the curve is 
approximately the (Gum) of the areas of N 
trapezoids of width h. 
nh=b-a 
‘The area under the Lrapezoid 
Ss U— Ff Gas(edh) 


~ Lg (atkh) 


+ (k-i)h Laikh 
1s given bu: 


RU ¢(a+Ceph) + F(a+ kh) | 


rhe total area under the curve from atob 
is then approximately 


A=h poe +f (ash)+f (arzh)+ flash), .4F (a+(n-h)| 


Notice that at No Time aid the word Integration’ 
pass my lips! Yrou don’t need to know calevlus to 
find OY COS tis way. he above procedure 1S 
obviously extreme\y easy to understand and to 
proqram. 

Some people ca\\ the preceding formula the 
Trapezoidal Bule for Integration. Of course, by 


eee § Diese eee 0 «6 EER 0 | ESE 


this name,it is much more difficult to grasp. 
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B\nother method for area estimation which 
is still easy but slightly more accurate than the 
Trapezoidal Rule involves replacing the 
straight line tops with parabolas: 


a ath arzh hath 
4 
areh ats carraieh b 


Yfou pick sets of 3 points off of the curve 
and then find the parabola that passes through 
them. (‘Irhere is only one!) A\fter hairy 
algebraic tedium, youfind that the area 
under these parabolas Is Given exactly by a 
simple formula. DF the points are even| 
Spaced a distance h apart, the area under 
the parabola from atoar2zh is given by 


kis (a)+ 4 (ath) + F(aseh) | 


Similarly, the area vundey the parabola from 
ar2h to arth is: 


RY F(as2h) +4 (ar3h) +24 (ar4h)] 
and so forth, 
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orhus, in toto, the area under the pavabolas 
approximating the curve from a to b is! 


AREA=2[£(a)4+4f (ath) +2 F(as2h)+4 f(a+3h) 
+ 2£(as4h)t...4 26 (axn-2yh) 
444 (a+ (neiyh) 4 £() | 


if there are n divisions. (Be svre nis even!) 


ea pn mmaemeamatind eed paneer as 


Rule Soy integration. 


(This integration stuff is easy! 
ee 


It really 16, $0 try a few problems. 


Pirst , write a numerical integration function 
suv program based on the Trapezoida\ Rule. 
Once that works , write a second function 
Subprogram based on Simpsons Rule. Whese 

wo programs shovla each be able to integrate 
onother function sub program named ‘FUNK, 
Using any desired number of divisions and any 
oven starting and ending points aGh, 
Next write a mainline program, ‘I his 
Main program is goin to compare the results 
of the trapezoidal rule with those from Simpsons 
rule, [It also is going to see what effect n 


PROBLEM 8-4 


has on the aceuracy of the results. Here's 
how your main \ine is to work? 


Integrate by 
Trapezoidal 


Integrate by 


Simpson 


Store answers 
im an array 


D 
oo 
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For each h, calevlate 
the ai fference between 


the two answers 


Starting with 
n=4, cad levlate 
the change mn 
each answer 
from the 

pre ceding one, 
That is, Compare 
the trapezcidal 
results for n=l 
with those Loy 
N=8, andso 
forth )for both 
methods 


Print everythin 
with clear 
héeadinas 


Now, test ovt your programs lou seeing if 
yov Can integrate the Following functions 
between the aiven limits. You can either 
write a separate FUNCT progvam for 
each one and-yun them one at a time, ov 
else you can try using an External statement 
if you know what that is. I suggest running 
them one at Qa tine Cor now, 

ry these: 

Y= ay from XZ O to X=2 


(Srhe answer shovld be 0.69315 ) 
| 
= rec from X= O to X= | 
(GFhe answer shovld be T ) 
“x2 
y= e€ e from X= O to X= 04 


(1 havent the Foagiest idea what the 
right answer is but L wovldnt be 
sur prised if it was 0.389 something!) 


y=(X"- cos X) e* From x=~\ to X= 


ae 
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Hlere’s a nominal *nomial Tov Poly: 
x3 14.275 x° 451) x ~54,3¢8 
Any polynomial of deqree 3 can be written as 
a,x? + a,x° t a,X% + a4 
Whe a’s stand for constant coefficients. Por 


instance, in the nominal > nomial] we have a,=l, 


O,= -14.275, and so on. 
Whips he ae ek) 
the thiva degree ¢ 
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Trhe valve of a polynomial of anu 
degree N mau be computed fora particu lar 
X Using the fo\lowing algorithm: 


Set polynomia| valve 


to zero 


Set a counter, k, 
to zero 
Add \ to ke 


Replace the polynomial) valve 
by Ap plus X times the 


current polynomial valve 
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“A sadder aed a 


wiser man 
He rose the 
Mu 

morrow morn, 


~ Colerid ‘Ge 
The Ancien t Mariner, 


Write a function subprogram 

3 based on the preceding algorithm. 
Your function shovla be able to — 
valve of o\unomia\l of dearee less than 
100. AS? “6 the function should be 
the value of X, the degvee N, ana the array 
of coefficients, 


O.K., Chief! 
TVve done that part, 


Now write a main program to vead 
in the polynomials and print out the results, 
Jt shovld read in the degree of the 
polynomial ) followed by the coefficients Arhen 
iL shovld wate out the polynomial, After 
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that it shovld read in a valve for the 
variable, then evaluate the pol nomial 
and print the valve of the variable together 
with the valve of the polynomial, Keep 
looping back ana readin in new valves 
for the vaviable til Yyov read a valve of 
-494.914 Phat 1S a Tag" which SAYS 

‘Quit evaluating this polynomial 

and go on tp “the next one.” 

Pe 5 
499.99 
\nlere are the polynomials to try, 

X2-5X46 ot XE-2,-15.7,0,2.18,-3,4 


I4 GX BX°-3.2Xt1Z1 at X= 1.6,2,-2,18 
XO 7 325 X2412%63 at KE-27.25,-9.3, 18.6 
2XP43N RH 4.24 4XU186 at X= 0,86, 1.2, 4.15 


Xt- 526K 44 ot X5-3.6,-2.1-1.5,-1, 0, 
| 552.1, 3.6 
On Covyse, You Can check some of these in 
your head .4fhe First one Comneés ovt 20 at 
X= gar the last one is 4 at X= 0. | 
Make VP SOME test polynomials yourse) § | 
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‘And nothin 4 to look backward to with pride, 
And nothing to \ook forward to with hope 


-Robert frost, The Death of the Hired Man 


WN ow, suppose you wanted to find the reots 
of one of those polynomials you were playing 
with (Yes, I know you never wanted to see 
them AQAin , but just SUPpose You wanted. to 
Sind the rects of one of them!) 

het’s assume you Know that the polynomial 
has Oo root somewhere in the intevva| 


X,4X<X, 
How wovld 1 Know that? 


For example, 1f you evaluated the polynomial 
at X, and at X, and found there was a 
change of s1qn, then yov"A Know there 


e a root between %, and Xz since 
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the polynomial 16 Continuous, 


aga oe Me. How can qT 
get from plus to minus ¢ 


1 reckon You have 


to pass through 
Zero +o get there. 


lynomia/ Va lve | 
at X) 


Cl 


AA quick and “ier Lasy se wae to Find 
Xx* is to evaluate the polynomial ata 
point Xz, midway between X, and Xo. 

_ X + Xo 
37 2 

A\ssuming the root, X* | 1S not exactly at Xs 

(and itwovld be pretty amazing if we 
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were that lucky ) then we have two 
possi bilities for where the voot is \ocated 


y-Cositive) 


. “rhe voot X™ lies between X3 and Xp if 


the valve oS the polynomial has the 


n ak Ky and at X3. Example: 


in the case shown above the Sign is 
positive at both X, and Xz , so the root 
lies between X, and X,. lad both 
Signs been negative , the voot wovla 
still be between X, and X3. 

He ‘rhe root X*™ \ies between K, ana Xz 
if the valve of the pelynomia) has a 
different sian at Ay from the $\4n at 
X3. 


Same $14 


see 
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When we started, we didn’t Know wheve 
the root was, but we knew it was somewheve 
between X, and Xz. That span \s called 
on “INTERVAL of UNCERTAINTY.” Now, we 
Still don’t Know exactly where the rootis, 
but ovr interval of uncertainty is only 
half the size it started out. We can 
repeat the same procedure over and over 
until the interval of uncertainty is as 
small as we like! 
hat’s move, this same 
technique can be used 
to find roots of Jang) 
continuous function, 
not just polynomials! 
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Urhis method is genera\ly called the 
HALF-INTERVAL, SEARCH technique. 


I don’t want to 
ae immodest, 
but you might 

mention this is ako 
Known as the 

Bolzano Bisection 
Method 


true. A\ny Wau , eschewing sustained 
the interrogatory and 

beating arouna the bush, 

based on this method. 
which will take in 

ending valves For an interval 
, together with a tolerance 
valve o§ . Program 

or 

the tolerance. 


that Valve of X aS 


it, have the subroutine 
1 endpoints on the 
along wi th the 
function at those endpoints. 
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CiBoth ends should generally be olmost Lhe 
Same point and Khe function valve should be 
nearly zero at both ends of the final 
interval. Use E format to write out these 


small numbers.) 


Hint: Lf Wov want to 
Know if two numbers 
have the same sign or 
different signs, multiply 
them together. L€ the 


result is positive , they iw 
pave the same signee / "Us \P 


Yfour subroutine should be able to find 
roots of (any) continvous unction, given 
vpper and \ower bounds on the woot, Now, 
to test ovt your subroutine , write a 
main line which will vead in polynomials 
together with tolerances and bounds on the 
locakion of OW voot , Using the polynomial 
program from Problem is find one OY Move 
roors Sor each polynomial. \s a test case, 


find a root of 
X7— 14.275 X2+51X- 54,368 


to within |Ixlot. PX voot lies somewhere 
between X=7 and X= 40. 
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JAIM BIN OS 


Flowch art and write a computer 
reoqre™ Lo calculate the answer to the 


fo\\ owing quod libet: 


How many bambinos did 


i, oes have ¢ 


Print out your vesvits ina matrix, S\X 
bambinos to a \ine. 


Hah, hah, the joke’s on you! 
The answer is‘ None ko speak of," 
Bernhard Bolzano may have been 
a Father, but he waswrta Dada. 
€ wag an Avstyrian 
Catholic Priest who 
lived from 1761-1848, 
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The WH omework Gets Cheesier and Cheesier? 


EX field mouse, & church Cor synagogue) 
movse, and a doormouse find a nice 
triangular wedge of Swiss cheese. Such a 
Govda find cheers them Up because they 
were Feeling kind of blev. (But 1 digress.) 

Gamembert to the subject at hand, each 
mouse starts nibbling at a different vertex. 
The field mouse tunnels towards the 
door mouse, tye doormouse nibbles towards 
the crunching sound of the churchmouse, 
and the churchmovse heads towards the 
fie\d mouse, Y wovld have included a 
tik mouse \ut was afraid to offend my 
Secretary, 

Be that as it MAY , write a proavam 
using complex arithmetic to caleviate and 
print ovt the loci of the three mice 
(and the ovtline of the cheese). UF Yovy- 
computer has a plotter, use it to plot evt 
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the paths the mice Fo\low., IF You don’t have 
a plotter, print ovt the (X, 4) coordinates 
of each step and plot the paths on gvaph 
paper by hand. 

E\\so caleviate and print ovt the 
— ae a three mice travel 
assumina equal nibbling speeds for each), 
Assume ci cheese is bedialecs, one 
mouse measure ona Side. 

Ts the meece’s locus sensitive to the 
nibbling speed! Size of nibble ? Brand 
of cheese? 

Remember: at each instant, each movse 
always heads straight towards the other 
rodent. Since the other mouse is also 
moving , the mice end up tracing ovt 
“Pursuit Curves,” Put tick marks on the 
curves after equal nibble intervals, (Say, 
every ten nibbles.) Athen yov can check the 
accvvacy OF the caleviations by laying OW 
straight edge on the plot, tangent ata 
tick mark. Tf the program is correct, 
the tangent to the door Mouse's path at 
the 7TH” tick mark shoulda point straight 
to the 7™ tick mark on the path of the 
churchmouse.Whis kind of check is the 
prooS of the pudding, ov the cheddar in the 
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cheese. Ghecking the progvam is as 


important Os writing the program, 
Pa 1 Q 


, © 
a” =e 


3 y 
“er 


Tangent at q 
4th tickmark 
points straight to 4° mark on other heus! 

For variety , you might try 4 mice ona 
squave, 9 Generals on a Pentagon, etc. Unless 
Your program makes USE of symmetry , You 
neea not VseE a requiar Figure, althovgh 
some mMovse might Come w \ast. 


,;mm (Mouse Measure ) = 
10 straightened mouse tails 
in length, by the way, 


Use Your compu ter system's re ference 
manuals for information abovt plotting and 
also abovt com plex arithmetic. Be suve to 
Type a\\ home-made functions that “Return a 
complex value, D§ uovr sustem has it, You 
might try using the IMPLICIT COMPLE X 
statement , just to try it out. 
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/ Whar? Did you say 
‘NT haven't any iden 
how to begin ad 


Rts really very simple, Neep 
track of the mice using a 
complex coordinate system: 


CMOUSE-DMOUSE 


X 


Erhese are cheese- fixed coordinates . Suppose 
the doormouse’s position is given ,atan time, 
by the complex vector DMOUSE. A unit vector 
from the door mouse towards the churchmovse 


wee Oe : }CDUNIT =(CMOUSE-DMOUSE) / 
CABS ( CMOUSE -DMOUSE) 

UF his nibbling speed \s ‘BITE’ then his 

next position will be 


DMOUSK= DMOUSEYBITEXCDUNIT 
See how easy it is 0 7 
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Wore’s perience Gauss Against lhe Grain? 
(Seidel vp Lo it slowly! ) 
Suppose You had the equations 


eX, X= Il 
X,+3Xz = 18 
YYou covla rewrite the equations O§$ 
- {| _ Xo 

X= a @ 
~o eee x : 
Ia-X%, : 

Ke > aa a 7 


~ noe wwwenwew &* % 


Le sim ple iterative way te solve these 
equations for XK, and As is as follows: 


Guess oa value Sov K, 


Try itin equation(A) 
Take the resviting XK, 
and try it in equation ® 


Take the resulting 
Ko and eee 
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Bach time thy ough the loop, the valves 
would get closev to the trve roots. When 
they stop changing by much , they are 
close to the exact answers. 

One Wau to picture the sitvation is 


Vf this doesn’t work, then You are going 
the wrong way around the loop In other 
words , if instead of converging to the reots 
things start to blow up, then Yov should 
be solving equation (A) for Xz and ®) for 
X)- 


This whole thing 1s Stupid! 
I can see by inspection the 
answers are X,=Iand A,= 5 
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— 
AnadTI can see by 
inspection you are 
an argumentative 
trouble maker! 


Avctvally , it’s trve. Por two equations 
and two mknouns this is a stupid and 
cumbersome way to so\ve the equations. 
You can easily solve the problem in closed. - 
form with no need for iteration. 


So what's the peint ¢ 


orhe point 1S) this Game methed covld 


be vsed if you haa man equations "n many 
unknowns. In a case like that, it would be 
hard to so\ve Lov the answéy's by aN 
Xexact” closed-form mathematical procedure 
What's more. the ‘exact’ answers that came 
out wouldn't veally be exact, dve te 
aceumviation of Yound-osf ana truncation 
evvors. This iterative scheme keeps 
improving the accuracy & the answers 
each time around the \oop, so its 
“approximate” answers m actually be 


better than the “exact” answers fomad by 
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solving the equations directly. What's move, 
it’s very straight - Sovward to program | 


Wow! What do you call this 


~ methed , anyway * 


] ca.\\ it ‘Dy. K’s Very Own Supev Duper 
A®F\ Patented Equation Solvev,” Avhe vest 
of the wovld calls it the “GAUSS- SEIDEL’ 
method, Just because Gauss was considered 
the “Prince of Mathematicians” and figured 
aM this stuff out o century and a half ago, 
I don't see why L shouldn't get some of the 
credit, 

EB schewing further obfuscation ; the 
denovement at issue in the present perquisition 
is... Cor, to put it ancther way, ) 

Here's what ya. gotta. do? 

Plowchart and write Oe preqram te solve 
two equations of the form 

aX, tbxX,=C 

AX, +X tf 
using Gauss~ Seidel iteration, where a,b, 
cde, and £ ave veal numbers to be read 
in Off of cards. Stop itevating and print 
out the answers when X, and Xo change 
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by less than 1X10 fvom their previous 
valves, ov if no solvtion has been found 
after loo itevations. (Print ovt an 
appropriate message in that case, 
ry at \east \0 sets ot data yinclvain 
the sample example from page 237 

or extra credit, genera \ze your 

ogvam to handle N equations in 


unknowns... 
en does this method work? 


hrhis method works when the coe Kident 
of a different unknown in each equation is 
lar e compared with the remaining elements, 
Solve each equation for the unknown 
has the biggest coeKicient (in abselute vale) 
hus, in the example, we solved the first 
equation Sor X, and the second for X, ,since 
2 is bigger than | and 3 is also bigger than I. 
eve’s a problem to try the generalized. 
program on: 
ISxX+ L6yYt+ 72Z- WE llo.ts 
l2X—174- 3 ZS+bl\we 72.1 
X+ Yt BzZt2we= -l2e 
X= 4HY=- 4 = €43.2w= 1982 


One of the answers 1s 4, another is -3.6! 
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Razors pain you; 

Kivers ave damp; 

Acids stain YOU 5 

And drugs cause cramp, 

Guns aren’t lawful; 

Nooses give; 

Gas smells awSul; 

You might as well live, 
-~Doro thy farker 


Suppose Yov had three towns, A,B, and 
C, connected by roads as shown on the 


Following map ° 


Distances between these towns covld be 
represented by Lhe following avYay a num bers? 


To Town 
A BC 
AjO 3 tS The street 
Ton J Bi3 08 Ese 
C 00 2 O is 2 miles long! 
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Arhat’s infinity ; dodo, andis because of 
the one-way street. You’ can’t take a single 
street from town C to town A, SO effectively 
the lenath ok the diveck youte from C te 
Aw infinite, 

Arhis 1s calleA a CONNECTION 
MATRIX. Le Row of the matrix corvesponds 
to a town at the Start of the roads, 


IOOEOOEOIIIDE GD 


of the roads. 

~Urfhe 4 in the matvix means there 
isa 4 mile road from A to C. The 
ZECYOS On the diagonal n the maxcrix mean 
there 16 no distance from Qa. town to itself, 
Pnfini ty in the matrix means ‘You can't 
get there from here ” by a aivect route, 
even though you miglrt be able to transfer 
in Peoria and take the mlK train from 
there, 
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INlow, suppose one wanted to find the 
hortest) distance from ‘B te A.|SPhe 


e@e00 6 OMHO* 


distance from B te any other place .Athe 
First column shows the distance from 
anyplace to A. 


Ar hus , the shortest rovte From BToN 
qo by any combination of paths 1S the 


2+ 0O , OF3, or 2+ W 
3 fe) CO 


Trhat is ; Sr 


Simi larly, the shortest vovte from Ctoh 
1s given by looking at the third row ana Fist colunn 
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pee 0 00 + O 
aa 1° minimum | 2 +3 )=5 
*~-{00 O+0/ = 
arhus we cov\a Gorm QO new matyix 
showing al\\ the One oY two voaa distances 


from anu town Lo any other, OAhis 
matrix wovld be 


B 
Two Road 
A 3 4 Connection 
Bs 3 02 Matrix 
C15 20 


Nf we wanted to find the shortest 
eae or 3 voaa path from any place 
to anu other, it would be possible te )ook 
fov the shortest 2 road path to any 
intermediate spot followed by the shortest 
| voad path from neve on te the destination. 
IN} oti Ce. however, that since there ave only 
three towns, you wovld never need to take 
Move than two roaas to get to any other 
place. 5N 
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Page Jand the xy IO) rotten expletive 
deleLeA@\ hasn't even toldus the question! 


Suppose we call the element in the ith 
yow and im column of a connection matrix 
Ciy. hen it’s pretty easy to calculate 
these elements. Suppese jin general, we 
n towns and we wanted to find the 
shortest 4 voad connections. Ifhen 


(4) ne 3 () 
Ci; = een (0, + Cp ij 
ke=\ 
Do it fer.all Shortest 3 read Shertest 1 
ossible Ss revte from road. revte 
tron |ten Lto Rk from k toj 


Pry will then be the shertest 4 yoad path 
from town L to town }. 


Very interesting ; but 
what do Yov want 


us to do? 


Write QW Function subprogram which, 
qiven alist & n valves stored ina 
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one — dimensional Arran, will Find the minimum 
value and ge that back as the valve cf the 
function. 

Trhen write a Subroutine which, given 
an m-~-road connection matrix anda 
one- path connection matrix will find 
the m+ | path shortest connection matrix, 
Trhis subroutine shovla work for any 
number of towns, rhe actual num bev of towns 
wil\ Vary , So pass it te the subrevtine through 
its argument list. 

Finally J write a M ain\ine which w} )) vead 
in First the number cf towns (1 to10); 
then the \- vead distance matrix Sor 
those towns. ead \ row of the matrix 
ola time , using an implied Ao. Since there 
ave \O towns maximum, use a. Fela S eclumns 
wide. IP. 8., the distances mat be Fractions 
like 4.5 miles. Cirhey miant even be in 
Kilometers , for that mattey, or they miaht 
be Trave\ Times instead of distances! ). 
(hen you Aa be finding a Quickest Path, ov 
Brachistochrone, instead of a shortest path!) 

+ there are YW Towns , UY ovv mainline shovla. 
vse the subrevtine m~| times to Sind the final 
distance matrix, Print everythin 4 out nicely and 
4? back and vead ma new Aata. set. Make UP Some 
good test data ! 
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‘As The Worm Turng” 
When last we saw ovr family oS 


simple - minded quadrile worms, they were 
livmg on a piece o% screen, a vectanqular 
mesh dipped in honey, A\ typical worm 
wovld eat the food along the wire and 
turn Right when it came to a junction, 
assuvmina food was available on the new 
wire, le the wire to the right had 
already been nibbled, the worm would go 
Gtraight at the junction. If the worm 
can’t go right or straight, he tumslé&st. 
F noSood is \eft on any branch, he 
expires, poor chap. 


TP 
afl 


Yfou are to Program & Plot the erratic 
course followed by each worm. Galevlate 
and print the total path length each 


worm travels, prioy to his untimely demise, 
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Incidental , to Keep this from being a 
cold, calcula 4, study of worms writhing 
in hunger, it would be @ nice, human touch 
to Name each worm, orhe Names oS the 
worms, together with their starting 
coordinates ave to be punchea on data 
caras, When a worm dies, it wovld be 
thoughtful to write his or her name 
beside the spot where he or she Kicked 
the bucket. F\\so , when yov write out 
how tar the worm travelled , Yov can 
refer to it by name. 

o read or write variable character 
infor mation, Yov need to vse A format. 
Since details of how Computers handle 
character data tendto be machine- 
dependent , Vil let Yov vead abovt it in 


yovr computer manval for your oun system. 


What? 
Gace “dependent? 


It’s one that the Computey 
lists on its W-2 Sorm. 
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Yfour progvay shov\d be able to handle 
from \tolO worms. Multiple worms should 
start from nice, symmetrical positions if 
you want pretty patterns. Instead A. 
reading na data cava telling Yow many 
WOYYNS to vse avrin OW particular wer ming, 
Keep reading data cards with the names 
of worms till you come te a cara that 
SAYS ‘GOWORM? SFhen plot thak batch of 
Worms meanderings ana start over. 

single worm woulda do this: 


Notice, he dies where he started, 
pavvre pe (re. 


Try the _- \\ owln starting patterns ) 
plus any others that “strike yoor fancy, 


rhese are all tested worm names. (Andy, 
Vicki s Larry , Jom, and Erland weve m long - 
suffering hard- working Teaching Assistants at MI T) 
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Let’s Get To The 
ROOT Of It All 


US. Patent 
#328,899 


Fr reel 


In 1885 My. Andrew J. Morrison of 
Buffalo, New York received U.S. Fatent 
No. 328,899 for an aerial railway run by 
gravity, fA\s shown in the illustration, the cars 
coasted along beneath A. continvovs cable 
suspended from balloons. arhe balloons covld 
Qo up ov down to achieve the necessary 
downhill angle. 

Und ovbtedly both My. Movrison ana hig 
potent have expired. IF not, he covld 
easily sell his idea to the Department of 
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Transportation. But since the patent has run 
out, we ave free Lo use the idea to Find 
Roots of Nonlinear Equations. (Ors COUYSE, 
the idea wi\\ vequive some minor modification 
for ovv pur poses. 

Suppose yov wank to find the root x* 
of the function Y=F (xX). 


£(xX,) £(>*)=0 


<U$ (x2) Ye Olde 
Root 


Xx 


We cov\a start by making OW pair of Gvesses 
ak the location of the root. Let’s cal\ these 
quesses X, ana Xz When we eva\vate £(xX) 
and 4(%,) we find we weven’t anywhere 
near the YoOt , but being too stubbern te 
quit we plow on. 

Following the suggestion of My, Waldo 
B. Secant of Krotz Pprings, Lovisiana, 
we make an Kdveated Guess that the 


function can be closely approximated bya 
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straight \ine in the vicinity of (X,) and 
(x2). We Pass a stvaight line throvgh £(x,) 
and £(X,) and slide down it toa point 
X, which is closer to X* than was 
either of ovr two starting guesses, 


Encouraged , we do this again ana again, 
each time passing the line throv gh the 
Ae-\ * guess and the foth qvess to predict 
the k+1* guess, 
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As you see, it doesn’t take us very many 
guesses to slide right down to the answer. 
We probably covla have gotten there even 
fastev it some jokey hadn't wrapped 
barbed wive around ovr styraiaht \ine. 
ZAnyway , when YUdV $ee that lf(x)| Is 
within some reasonable tolerance from Zero, 
quit and call your last xX the answev. 

In honey cf Waldo Secant , this is 
called the SECAN'T METHOD. Grr 
using it to find the roots & the following 
equations * 

eo * - Sin (BS) = 0 
(A root lies between Oand \) 


2Sinx= 3 


Cor OF 5-7 2SinX ) 


Enough, already, 


A\M right ; 
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The Distillation of Human Knowledge 


<——~O) 


CEE 


A g 
G00056000 
y eote 4 
one ‘ 
we NED 


LTA cA h— 
Figure 1 
An experimental catalytic convertey being 

developed in the hills of Kentucky in a top- 
secret research project converts sour mash 
into a lead-Svee Sve\ which can put new lite 
into Your Lived old carburetor. his work 
1s being conducted under a jeint vesearch 
rant from the MaSia andthe Cen bral 
Thalapaia Agency. 

Exact details of the device ave, of corse, 
classified. llowever, Oo Key component 
the system iS a long Square pipe with a 
second Square duct yunning down the 
middle. (See Figure 2.) trhe space between 
the two pipes is packed with fermented rye 


mash, A critical questi on is 
“Whar is the exact temperate 
distribu tion n the material ? af 
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Linear 

Tn crease in 

Temperature 

tiie 

32° 
oe 


Q Q 
a) 
Constant 


Tem perature 
Bath ?) S T = 32° 
32°F gi! 


Figure 2 
Drawn by: Moon bean MeSwine 5 In ventor 
Signed 2’ X-—+ZHis mark 


We need to use the computer to calculate 
the temperature distribution in the rue 
packed between the two SQUAVE ducts. rhe 
outer pipe is 4 inches square and the inner 
pipe is 3 inches square and runs down the 
middle. 

Es part o§ the Likker Condensing 
Dystem, Lhe \owev part of the ovter pipe 


1S submerged nm @ Constant Temperature 
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Bath, consisting of ice cubes floating in 
Kickapoo Jou Jvice, this keeps the 
whole ovter wall of the lower part at 
exactly 32° ¥, 

L\s part of the Likker Kggstrakter 
System, the innev pipe is filed with 
steam which keeps its walls ata constant 
2l2°F. 

rhe top surface is part of the 
A\Ky Regenerator & Revenoo Offiser 
Roaster, Arhat surface is Kept ok a 
constant 00° by the hickory five. 

“rhe sides os the outer pipe vary in 
temperature , from 100° atthe top down 
to 32° at the height of the ice bath. the 
bath comes 4" vp the side of the pipe, 
as shown. irhus, the bettem 4" are at 
32° and then the pipe wall increases linearly 
in temperature up to 100? 

hus ) we know the Boundary ‘Temperature 
all around the walls in contact with the rye 
mash .4rhe Catcher in the Kye 15 
How can we Tiqure the 
temperature in the materia\ 
itself, knowing the boundary 
temperature s ° 
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? \\ tell You how ? Nts not ver hard ak 
a\l, Just Keep an eye out Sor the J eno 
Offiser while L tell You. 

Suppose we weve to subdivide the Space 
between the pipes into a lot o& small 


Square regions aS shown below ¢ 


First Column Twentieth Column 
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Notice added anextra row of Squaves 
inside the inner pipe and ovtside the ovter 
pipe. WN tell uov abovt those Phantom 


Savares later. 

: Now, if the SQUAYES aYe Gmalenough then 
it i6 reasonable to assume that the 
temperature is uniform within each square. 
In othey words, we covla say the Square in 
the 7% vow and 12% column is at 53° 
ana we wovldant need bo Say how the 
Lemperatuve varied acvoss the square. With 
sm al| squares , there wouldn't be much variation. 

AM\so, common sense tells us that the 
temperature in each squave wovla be the 

Averag of the temperatures at the 
fovr neighboring Squares. 


AN typical Spoare 


. . Im. in the 
Ra and 
| fBens| iCal 


Thy = z ( Tei + Tiat3 +1i + Ti ia 
MAGIC FORMULA 
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Looking at the figure on page 258, 
we see we covla set Up a big 20 bu2o 
matrix in the com puter ana use it te 
store the temperatures of all the little 


Sq Vvayes. 


But we dov’t Know the 
temperatures of al\ these squares! 


Arve but we Do know the te m peratures 
o§ all the Phantom Squares avound the 
boundaries. We can fill hose valves into 
ovr matyix, since they arent going te 
change. 


hat abovt a\\ the Squares 
between the twe ducts? 


hrPhey’re the problem, aren’t they? Just 
to ger started, lets arbitrarily set 
them all to some reasonable valve 
between 32° and 212° SAY 90° or so, 


You call that Science? 
It sounds pretty random ty me! 
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Mt is pretty random because it 
Doesw't Really Matter. Aphis is just OW 
starting quess. Now we’l| itevatively apply 
our Magic Formula and use it ovev ana 
over to improve ovr estimates of each 
temperature. When we’ve throvah , we'll 
Know Very closely just what the actva!) 
temperature of each squave really is. 

Start with the first imside square and 


iteratively step from squave to Square, 
applying the (Magic Porro) to each 
inside square. Stove the new temperature 
Sor the square ‘in place & the ola valve. 


What do we de when we 


reach the last squave? 


Bo right back te the first Square 
and ao it all over again, A\tter a number 
oS passes , the temperatures of the sqvares 
won?t change much from iteration to 
iteration, so you'll Know you have pretty 
much found the final temperature 
distribution Yov weve looking Cor! 


Dan that again ts 
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O.TKY. Hlere’s the process. You cycle vom 
square bo Square applying the Magic Yormula 
OS you go. Rac AME y befove recording the 
new temperature for a square, see how 
much oF a change thatis Srom what it 
Was be S ove. Don gi bother recording this 
change unless itis a bigger change than 
ONY other Square made Avring that Pass. 
Ashvs , during the pass, you keep track of 
the biggest (abselute valve) change in 


temperature of any square. 


That’s easy en ough 


When you finish with the last interior 
square, check what the biggest Change was, 
14, during that Pass, Some square Changed 
tem peratvre by more than 2°, o back 
ana do it all “over again. [AStev a while, 
the valves will settle down and no 
sqvave will change by move than |/2° 
during a pass, 

EKt that point , quit ana write ovt 
the final tem perature matrix. 
be —— z —_ time, $$, ana space in 
Computer making use cf the syummet 
about the a De a fi the duets, el 
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Hint #2: Make svre you don’t try 
applying the Magic Formula. Lo one of the 
boundary Phanton Squares! orheiv 
temperatures never change! 

Uncidentally, You have just so\ved 
Laplace’s equation, Vru=o . US you 
Wan ted O- MOVE ACCUVALE ANsSwey , You Cov 
vse smaller SQUAVES ana a tighter tolerance, 


Solved it? Tdon’teven 
Know Mr. Laplace! 


Hint #3: Watch ovt for the Revenoe 
OFffisers! 

Gracker Jim sez he was ovt shi n’ on the 
lake t other evenin when two Revenoors 
rowed by him with a still they had just 
contiscated, JAn ole man and a bo 
passea them rowin’? in tother divekshun. 
They rowed on quite a WAS withovt 
sayin’ OW word, but Tinally the boy spoke. 
Gracker Sim Sez they WUZ purty fay AWA, 
but voices Carry over a lake in the still & 
the evenin. He covla make ovt the boy 


ae Paw. ls thet ovr ?n? 


PROBLEM /7~9 


Iiey punching FORTRAN? 


Statements... 
Columns |-5 Statement Mumbers 
Column 6 Continvation Marks 
Columns 7-72 Statement Body 
Gomments,.. 
‘Cin column |, vse rest of card for 
the actual comment. 
Data Gards eee 
Punch data anuwhere in columns 
\- 80 as specified in FORMAT 
C909 2909 & 09409 49 OF HO 00 WNW WO OW 


Gonstants 


Real or Floating Point... 
2.653 or - 37.29 or \.2E-lé 
(Limits on magnitude are usvally 
from abovt 10% to 107°) 
Inteaer.., 
165 or —623 
Depending on computer, Maximum 
Magnitude mau be 32767 or else 
2147483647 as typical limits.) 


Complex... 
(0.,1.) or (2.1495, 8.623) 
Logical... 


sLRUE. or .¥ALSE. 


NUTSHELL 


Variable Names? 
\ to G characters long 
jst character must be alphabetic 
Reals.. 
ALPHA ,X, BIW, DELTA, THETA 
A \st Ae eter can tbe T, J KL MorN 


unless variable name 1s typed “BEAL 


erin SLASTX, K,M80, I, LAMDA 
a een en alek vel, 5, K LM, orN 
unless variable name is typed INTESER 
Complex ms cae Dovble Precision Set Cas 
Any | ega\ vee Bic name a iv Auk 
be typed 
CIM OD 09 09 WAT DOD WM MDW MWMWM MWD WwLH 


iE x pressi ons° 


eal... 
PL (2.6 4SQRT(X-Y)) XBl2/(C-X)+ 2.4 Y 


Integ er, 


(4 T}KKI5)X MTOP /NTER 


Complex 


et 6, -2. 8) x CMPLX (XY) /(VEC+(1.5,0.)) 


Lo eal.. 


X.GT Y-2. -AND.Z.NE.ZBASE 
Expressions ave Combinations of Variables, constants, 
fun ctions, operators, & paventheses 


NUTSHELL. 


Operators: 
Avithmetic..,. 


~AND. , .OR. , .NOT. 
Relational... 

EQ. ,.NE.,.GE. , BP re <i are 
Hierarchy... 


Evaluation of functions; exponentiation; 
multiplication and divisi on; addition 
and subtraction : relationals (.EQ. ; 
NE.,GE.,.LE.,.6t.,.LT.); NOT. ; 
AND. 5 .OR.. Leftto right, 
parentheses override, 


09 08 09 YW WO CO OY 90 TM OC 0 209 09 AO 0 OC OOO 


[\ssignment Statements? 
Arithmetic... 
ALPHA =(PI+SIN(THETA)) /2., 


Arithmetic Expression 
Means " replace by “ 
Single variable name 


Logical. ; 
BOOL = TOL.LT.VALUE 


Logical Expression 
Logical Variable 


(68) VUTSHELL 


Gontrol Statements 3 
Uncondi tional Branch.. 
GO TO 50 


Arithmetic \f. 
Lk (X- SQRT(Y+2Z)) 20,10, 30 
Arithmetic Expression 
[f negative 
Lf zero 
If positive 
hagiee. Tf. 
TF (X. GE. TOL. AND. Y- SQRT(X) 
1 .NE.0.) “READ(G,60)A,B,¢,4 
Gets done if expression 
is TAVE. 
: ‘IF (BOOL .EQ. .TRUE..OR.K4.EQ. 
. 2*K8) G0 TO 80 


Comp vted GOTO. 
ko TO (10,30, 30, 10, 70,50, 10), K4 


Means” if KF a A 5, bone! 
to statement number 70% since 
that is the 4% statement number 
in the list, 

GO TO (50,20,20,10 20), M 
Means ge to statement #50 jf 
M equals 1; to 20 if Mis 2,3, or5; 
or to /O if L0 equa/s 4 


NUTSHELL 


DO LOOPS... 
DO 30 K=1,M,3 


RR: 
to ing Valve 
ae om bonen Starting Valve 
DO GO I=1,10 Positive Lategers 


DO 60,3 =5,150,5 
DO 40 M2=11,3 ested 


$ 4 § Loops 
40 { fF 


60 X(1)=X(1+2) +DELTA 
Increment & test occUrs at end of 
DO loops, DO Counter can be Used 
but not Abused within loop. 


Coc 798 0D tO fC! Ow OO 029 MH 19 WM WT] Dw CL CT 
ZArrays ° 
ubscripts.. , 
Moe. ek 2xMAT +3 _ahert , ete. 
Legal subscripts : Must be one of these: 
Integer Constant 
Integer Variable 
Integer Constant Integer Variable 
or one of these two *Iateger Cons tant 
Subseri pt values must be positive 
inte gers, lSSubscript $ Dimension Valve, 


(270) VUT SHELL 


Storage cas 
By columns, Wi th jst subscript varying 
most rapidly, Storage is reserved by 
Dimension Statement at start of program 


which specifies Maximum $1Ze of subscripts, 


xample o ee 
DIMENSION X (50), Y000,2), MMXB2,5) 
| Integer Constant -Lavgest expected 
M 


iv 


valve possible for svbscript, 
X(1,),K-2) =2x MMX (J,2, K-1) 
Cae 
SWbhscripts must be positive integers. 


CMC WW WW WVDWMNMNDOVDwOW WT NY WT 7 ~S 


Japot /Ovtput Statements 2 
READ(8,10) T3,2Z, (ARAY (1),1=1,10),M. 
of ut Input List: variable or 
array names, arrar CLEMEN 5 
OY format) or implied DO’s 
FORKAT 2FIO4/2 1X,5F104/),15) 
Field Specitications 
WRITE(3, COALPHA 13, LL, MTX3 
but it Unit Ovtpvt List 
LOcFORMAT (TIALPHA 1Sb!, F10.4/ 
HOX,'T3 ISb’ EB 15.6,"bAND LL ISb’ 
2,116/"bMATRIX ISb'//4(3X,15, 1OX 
3.15/)// "THAT'S IT") 


NUTSHELL 


Field Specification Godes. 
End of record (\ine ov card) 
CX DKie 6 columns 
110 Integer field ClO columns wide) 
F\0.4 Real Field (4 figures aSter point ) 
EI5.G ‘B notation, 6 Significant figures 
Ym! Literal field 
SHTEXT Litera\ field 
A2 Character dota 
L\0 Logical field Not inText 
DI0.4 Double precision 
G\5.6 Generalized field 
Speci Fication Statements3 
REAL LAMDA, NUM (50) 
LOGICAL DEBUG, BOOL 
4ype \ INTERGER ZVAL TOTAL(2,5), SUM 
COMPLEX VEC,DELTA 
EXTERNAL POLRUT , FUNK : SUBPRO 
DIMENSION X (10), SUM (2,2,10) 
COMMON X, Pl RAD,DEG 
EQUIVALENCE (TL, TONE), (MAX, LIN) 
DATA DEBUG, X /.FALSE., 10 0.0/ 


After these come the Statement Functions, then 


the body of the Program ana finally the 
END statement. iti, . 


NUTSHELL 


Subprograms 
Statement Function 
SUMPRO(A,B,C C)= “(M4 B4C)/(AXBX C) 


a Lx pression wi thevt 


A P -g v ments iatlons | elemen ts 
Fone tion definition most. ge at Lop 
of ac tva/ Pregram that will Use it. 

“ on Sub pr. 


=2,% SARTCX- § SIN (ALPHA) *X ) 
Las tion References» 


alae ic ale a LAST (X\,%2,ND 


RETURN 
END 


S ubrovtine Subpro yvams, 


CALL EVALUS(MTRIK, INVES 50) 
CALL SWITCH(P.Q) 


—_—"s resuv/ts are returned throvgh 


vment list t, not rovtine name, 
SUB OUTINE SWITCH (A.B) 
TEMP=A 


B=TEM?P 
RETURN 
END 


NUTSHELL 


ONOD)Px¢ 


The tollowing topics NOW oY may not be 
found in this” edition. If present, they 
should. or shovld not be taken seviovsly, 
because they are or are not basedon reality. 


opic Kage 


A Format... Aa) 4 
Adi jv stable Dimensions. eeeeunigeeceneceees a 154 
Alaorithms... ineatimionsaotesieace ree 
Alka Seltzer... shapeiendboteiasenaenancenll 

A\po.... pea puarintadenesnes 

Ampersands (6)... oceee ce ceee cessecse 4 

Areas of Polygons... .. 206 

Arguments... 18, 136, 145, 158, 162 4273 
Arithmetic Assignment Glatements........6 268 
Arithmetic Expressions.... een NT 
Arithmetic Lf... 200.02 37,269 
Arithmetic shinai eT eee: 
AYPOUS 0... sens nner M27 
Array Storag ae apd Ua 147 
Avtichokes (Bec Chi Na, Whe. om 23 
Assignment Statements... 6 125, 268 
Bananas... AZ 


Bananas ois Cee ayes Uncle Milty).. ae 87 


INDEX | 


at hae ee rae OT 


Bolzano’s Bambino... eal 41, 193, 230, 232 
Bostonians aii bon aes 52 
Branches... eieighageeuioe: 34, ft 
Bin caw Dugwers 

Butterflies... bocce cee cececcecses semen, 649 
Cahrtah Timmy Lec eccceceeccacee cscs ane OF 
Cal| Exit. vec cevtetteceee se  143 
Call Statement... emnwenewseden O50 73 
Carriage Control... ves cessecececeseee 04 
Characters (Literal)... crovenernn BAG 2TE 
Characters (Li berary).... a |-273 
Characters er M2. 34, 51; 4 
Chebyshev, siaaideie 

Child, Julia. - 
Com bor tin 4 Words £, or When en Your RogranDies 
Comments... ..214,266 
Common Sl lewments... we "160 
Complex Arithmetic.. 127, 235, 267 
Complex Constants... sessnrne ned OF, 266 
Computed GO TO.. nr 4 | 
Conditional Beanthes... ese OT Tice 
Connection Makvix ....0........c ccc cccec cece 243 
Constants 0.0... cece ee eee 646,129,266 
Continvakions...... ccc ccecce cre cesereees 30,2 
Conlinve ... 2... cece cccccece cu ccees cess sees 78,270 


INDEX 2 


Contro\ Statements............34,3741,74,267 


Cbbiaal Nombers .....c.seaisroGisnee.. alae 
DaS$od)\s.. causa snbtewhesiintee’ss dt eee 
Dahan Cards.. ; rai rie per ic 
Data Statement. tas teil eric 
De Morgan... lee ae oA 
Dimensions... oie 149, 152, aa 
Dimension ee ceck. 7 94. lid, 27 
Dimensions, Variable... ccc. 9458 
Distilling Cie. Sy as 258 
Do Loops .. Awes 270 
Dow’t Cry You" Warp Ye Nouri ie: 32768 
Drivel.. = het 
Drum es niet 
Dummy Arguments ( 6Sed e Arguments & Tiffs) 
E Format... 182 Alon 272 
End at ae 133, 135, 143 
Kaquivalence.... aaa eee a ee 167 
Krrors aes 11,16, 3) 
Bar\bives ee eee we eS 
eat Nebstion histyeea tev 2. Pie 
Exponentiation... Sa ee eee csawchigd Godt 
Expressions. atesssctonve coon | ghey 
External Statement eine. 2emyeel 
Format........ necesinsos centsicqeehhoeaee 


INDEX 3 


Fibonace).. vc ceceeweeussessnee.2OO 


Field Bees Fen ions: Loc sessesene 12 96,272 
ivema ee oy 
Flamboyant Bin swanks.. aetiireseseatesvonies aos 119 
Fleating Point Numbers... sev cveveese by @b6 
| shes ogg Oe Bb 
Floormats... eee veet tite bettie BS 
Flowcharts... iperiaertonceeene 41 10 
Foolish Mistakes... oie deas anes auesiden 1 
Formats.. oon $8,272 
FORTRAN iv Na Se hell. sssies REOOS OVS 
Functions... py S15 156,273 
TP Scimenka\: ere ey 58 


Funnu Answers (See' lites e ne hiidvs e Ov’) 
si a 4 
Garbage In , Garbage Ovt (See * Miso Answers” 2) 


Gauss ~ Sebiie| Method... 7 exe 
Gert Lost Statement (Sea e “Systems Frogramners) 
Girdles. 122 
Gol Te a Barry. Wdeitioer iene 
Go lo Bee iake.. View and Pm nk 
H Format... cu he dela anietal ae 
Hale, Wilhan.. oon oho dex 
Eis - Ribenin\. ee VM in idaaneeet 
Hamlet . seneitetionunate 
IN oleh oecssneriiy cknpaneieee 
Hhevon’s Formula.n...c.cconnnndad ean 207 


INDEX 4 


Homilies... occ ccc cceceeccseceseeccececetevsresessesees [BG 
Hypertension and the Computer usecase 189 
] Bovmat oo... cccccccccsececeeteecerneseuesee ee 9 6,272 
Tf Statements (Avithmetic)............37,269 
If Statements (Legical) enna 4 | 269 
If Statements (TWooical yo... cece SM 
Implied Do A YY) Se 271 
Implied NPC ier suzaeceicactonencmiammnsl Moen 
Indiaression . 20... ..cccc cess cece seas eee eet esos 5 
meey, OUR BUR cons atecentcenaentets §2,102,131,271 
Input / Out put OF Arrays... cece 102,271 
Lg a res: | 266 
Integer Field Speci Fi CAL ONS cccccc cece eeeee OF N4(4 
Integroh yper polydiffere nti a\ Projective Functions e: 
Integration by Simpson’s Rule.................218 


Junk You Should Know........... ec eeeeee. 1273 
Kavfman, My Sen the Dv..........14,54,82 124 
OY PUNCWING « -. oes eects cree terete reteeeenees 2 266 
Koala Bears... cccccccc cece cece eee creer, €By 134 
L Fiela Specifications... cece eee OTC 
Laplace Py QUAL Oiiescsccnra easscaumcassencsss C20 
Learning Por bvan oc cece cccccccc ee ccnteeee caves IBGE 
DOG re cacecetoni sic ecueeteept epee panietaenetven, 
Liberace .. oo. cce cece ccc ceseccee eeteenese nese. 200 


INDEX 5 


Library cemregtte: — tT 


Line Printers... seas 73, 120 
Lint Belly button... 4 
Titeca! Data .. _ 67, 249, 272 
Little - Known ee 17| 
Lobachevsky, Nickolai Ianovie 150 
Logical Assignment sande 125, 268 
logical Constants... 46, 266 
Legical Hxpressions... ain "268 
Logical Lf Statements. .41,269 
Logical Operators... _ oo 46 268 
Logical Variables ....... 0 ccc cee cece I25 268 


Long Time, No oiaila ned "Loops, Insite 
Loop - The - Loops... 


Mechanic’s Rule... 192, “A 
Mice, Field Church Atheist Door, 6 Other... "233 
Mini- Sum Avithmetic ... .. 242 
sc MILTIB 
Mixed- Modes... eecaeeeeese _ 
Mnemonics & Minemoricas.. 

Modes, salmon 7 64 in 
Momma....... ec eeeae ey \2| 


INDEX 6 


Monroe , Marilyn... cbeeeeueas 16 


Morgan, Henvy.. 44 
Names im Vase aivec & Svbprogramsz 47.267 
Nested Do Loops.. . repr 270 
Nested Lolied Op Leops pean bode oak oe H3. 114 

Newton’s Method .. ne. 

Nichomacus.. Kadai d-hant tuwe-adeeaeesniueeades 1 IO 
Nixon Richard... 34, 39 
No No’s (See Also Nixon ny Richard). q IG °3 I 

Nutshe|l\.. 265 - 273 
Octal ‘Touring B. .. 124 
Oiling the Computer (See Maintenance Handbook) 
Operators , Avithmetic .. we, 268 
Operators Big T ime (Sel Lance Bert). paces 178 
Operators , Logical ... 1 46,268 
Operators Relational... Leen case ceee 43,268 
YoSnutij. hs) 


Painful Humor ‘(See the IBM Manval #0286515) 
Panic Statements (See “Hardware nnn 


Farameters......... , — .. 158 
Parentheses ....0 cece cece ec cree cee ewes 13, 15 
Yorker Dorothy ... - 242 


Partridges, Prench Hens, & Mocking Birds .. WE 
Fause Bi le weat (We Don't aie Time For 1h) 

Feev\iar FORTRAN Featuves... 1-273 
Per Fovated Bons ...... 0... cece cece cc een eee eees 64 


INDEX 7 


Terforated Vieers (See “What do you mean the 
Susteyn will be down Sor two weeks 2) 


Ph.D Candidates... Leecccecesccesseersee, LOS 
Pimples... srcrnnnnnnnerenenes OF 
PI Program. ee a ee ee ee are C04 
Polynomials... st aie ate doeeeedtesacens COC 
Prime Weavers. er eee rn 207 
Printers... _ 13,120 124 
Pythagorean Theorem... aeees “160 
Programs... - 20. 19¢ 


Programming the MINSK 90: in ‘\@yOBBOVIAN: ’ 


a non-executable, uninterpretatable lanquage. 


ProneS...cccccccccececcccees 0 teavauceecusueeeeceues 45 
Punch Caves: sansduntinden eeaitiutdasnueeses Ot 
Pursuit Coeves.. seueioe wsetence Licecep ea cecsqerces OS4 
Queer Statements......... ccc cece cece eecesecers 135 
Questionable FORTRAN... ssammasnnicesing TOO: 
Rabbits.. esta biceacee deinen and COU 
Range of De... wa 

Ranges, 6.E. ae Hotpoint See" "Hove wares) 
Keada Bi tenieale x. a ar cy Ae ea 
Real Numbers........0.ccccccceceeeeeeeceee | bb 
Recherché Broken sneer Jo 
Recursion Relation... cccccsccseceee es 202 
Relational Onecdioce.. 43 268 
Relations. Inbred (cee “Rabbits” eee ...c0O 


Whu do th 4 have a ‘Syn Tax’ in FORTRAN ? 


INDEX 8 


Repeat Constants ........0..cccecces cece cece cee 68 


Repugnant Slotements.........cccee ccc ccesee wees lll 
Resign Statement ......0...... 0 cece ceeeee ees 135 
Return Statement... 137,273 
Root Finding... _ ell, 226, 23 /, "25 
Dalesman, Ga Wane: Texvelisae. vee 124 
Sample Problems & — see eaesee: 19 263 
Sample Program. shia aiden Aunaseen AWG 
Secant Met BA: co vctetceeteseeescer 1:25 
Secant, Mr Waldo B... icceuusmedeaneeee ee 
Shales peare: Mr. Wm. . jemenee cone es 44 
Simultaneous Kquations... Seva 281 
Sin (As in’Sine”, which isnt any onl. 
Slashes in TP emate.. 60, 2p 
Slashes in Tires (See “Hell °s Angels!) 
Society Pro Musica.....cccccccccccec cceusees S| 
Specifizati on vinoeshsalaaian ..-. 170,272 
Spiro Whatsisname.. yeqevsepeaceahetn OA 
Sart.. wpnedncinseaioces 10 
Rate Rost Algorithm. . aeaeeeincse Ie 
Statement Functions... ee “eves 156,273 
Statement Numbers... .........0006- 30,35,206 
Statements, Order OF... cee. 170,272 
Stop i ee oe 143 
St torage.. 199, 21. 145, 152, 162, 27| 
Struct wed. Programming ... me : 
SUbprograms..... eevee ce eee cea seeees "134,273 


INDEX 9 


Gubroutines..... cee cece ccc c cee eee cee ces [4),2 73 
Subscripts... Lec cesese cence ee 4 F598 270 
Sundry DL akemenks ede de eee seenestees |- 273 
au" bollic Variable Names......2,/24 ee re 


System P yroarammers.. 


T. Aw’s (Andy Vicki jhareyiTom, Sack . Heland).. "250 


Tea Cosy... rr 0) 
Tem seratove Distributic on ‘Peblew. geaeeee 255 
Things You Learn The Hard Way........+- 124,189 
(See also “Structured het aati 

Trying To Cope With identi . 188 
Tupe (Function)... lene .. 140 
Type Sia cements... sugeeeetevesa ones 125, 127, 272 
Unary & Ornery Operators... er 
Uncle Sam.........ccccccceceeecccececceccesseseas 14 
Unconditional Brancheés...........2006 34,261 
Variable Names 2 124, 26] 
Variable Names, Cute a THIS 1... - 
Vaviably Dim eacioned Ayvays.........eee eee 


Viokvostek Telephone Direchaeites. saan — 
Who to do till the Systems ko henien comes, 
What to do once he baa there.. . Vol.2 
“Worms , Quadri\e.. oes 248 
Wri te Statements... . 5, 27l 
Wrong!.. sien eae ea ences Al 


INDEX 10 


X Field Specification... 0.0.06: 60,272 
You Can’t Get There From Here Statement 


Zat’s Ze End, Folks! 


GALL EXIT 


END | 

* LINE #0015 ERROR #0047 INVALI 

* LINE #0023 ERROR # 0049 UNRECO 
// XEQ 

EOF 

UNDEFINED SUBROUTINES: 12 

X¥**% FORMAT ERROR #0034 ILLEGAL 
ERK 0004 27CA OVERFLOW IN $Q 
ERR FD AT OB7¥ FATAL, STUPID MIS 
EXECUTION SUPPRESSED 


INDEX [1 
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MONEY-BACK GUARANTEE:’ 


Tf you do all the programming 
problems at the end of this book 
and don’t learn FORTRAN in 
the process, you deserve your 
money back. 


* Offer VOID within or without U.S.A. 


